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FOREWORD 


This  report  describes  the  inhouse  effort  conducted  in  the  Mechanics 
and  Surface  Interactions  Branch  (MLBM),  Nonmetallic  Materials  Division  (MLB), 
Materials  Laboratory,  Air  Force  Wright  Aeronautical  Laboratories,  Wright- 
Patterson  Air  Force  Base,  Ohio,  under  the  contract  #F33615-83-C-5056  with 
the  University  of  Dayton  Research  Institute. 

The  work  reported  herein  was  performed  during  the  period  1  July  1982  to 
31  June  1983.  Dr.  Stephen  W.  Tsai  (AFWAL/MLBM)  was  the  Project  Engineer. 

This  report  is  a  modified  version  of  AFWAL-TR-81-4073. 

The  author  wishes  to  express  his  deep  sense  of  gratitude  to  Dr.  S.  W. 
Tsai  for  his  extremely  fruitful  guidance  in  the  course  of  this  work.  Figures 
1  -  2  and  Tables  1  -  7  have  been  taken  from  the  book,  "Introduction  to 
Composite  Materials,"  by  S.  W.  Tsai  and  H.  T.  Hahn  with  their  permission. 
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SECTION  I 
INTRODUCTION 


This  report  presents  a  FORTRAN  computer  code  for  the  solution  of 
composite  materials  problems. +  The  computer  program  can  conduct  the  point 
stress  analysis  of  general  laminates  including  hybrid  laminates.  The  effect 
of  hygrothermal  and  mechanical  loads  on  the  strength  of  composite  laminates 
can  be  studied  using  this  program.  The  analytical  formulas  based  on  the 
lamination  theory  have  been  used.  These  formulas  are  available  in  Reference  1. 
For  the  sake  of  completeness,  the  relevant  relations  are  presented  in  this 
report.  All  the  notations  used  in  References  1  and  2  are  followed.  Plotting 
capabilities  to  investigate  the  in-plane  strength  of  layered  composites 
based  upon  six  commonly  used  failure  theories  have  been  incorporated.  A 
large  number  of  options  to  obtain  failure  surfaces  are  included.  CALCOMP 
plotter  is  used. 

This  manual  has  been  written  to  provide  the  users  results  for  their 
problems  with  a  minimum  of  effort.  The  input  instructions  are  explained  in 
detail  and  are  supported  with  examples.  Material  properties  of  five  well- 
known  composite  materials  are  stored  in  the  program  and  can  be  used  by  giving 
the  relevant  material  name  for  pure  laminates  and  by  giving  the  material 
property  identification  number  for  hybrid  laminates.  The  material  properties 
which  are  not  included  in  the  program  can  be  used  through  the  input  data. 

In  such  a  case,  new  material  properties  will  take  the  storage  space  provided 
for  already-stored  material  property  data;  and,  therefore,  the  corresponding 
material  names  or  material  property  identification  numbers  are  to  be  used  for 
further  calculations  in  that  computer  run.  The  new  materials  are  to  be 
supplied  in  SI  units,  according  to  the  units  used  in  Table  1,  while  the  output 
can  be  obtained  in  English  units.  For  obtaining  results  in  English  units, 
a  command  at  appropriate  place  has  to  be  given.  For  pure  materials,  this 
command  is  ENGLISH  at  the  11th  column  of  the  material  card  and  for  hybrid 
laminates,  IUNIT  =  2  in  the  $LAYERS  card. 

On  the  basis  of  six  commonly  used  failure  theories,  the  in-plane 
strength  characteristics  of  multidirectional  composite  laminates  can  be 
investigated.  In  the  case  of  maximum  stress  and  maximum  strain  failure  theories, 

+Modified  version  of  AFWAL-TR-81 -4073  [12] 
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TABLE  1 

MATERIAL  PROPERTIES  GIVEN  IN  LAMDATA 


LMPI 

1 

2 

3 

4 

5 

6 

7 

MATERIAL 

T300/S208 

B4/5505 

AS/3501 

SCOTCHPLY 

1002 

KEVLAR49 

CORE 

ALUMINUM 

PROPERTY 

GRAPH/EP 

BORON/EP 

ORAPH/EP 

OLASS/EP 

ARAMID/EP 

Ex(OPa ) 

181.0 

204.0 

138.0 

38.6 

76.0 

0.0 

69.0 

Ey(OPa) 

10.3 

18.5 

8.96 

8.27 

5.5 

■a 

69.0 

Vx 

0.28 

0.23 

0.3 

0.26 

0.34 

0.0 

■a 

E  (OP a) 
s 

7.17 

5.59 

7. 1 

4. 14 

2.3 

■a 

26.5 

«  frm/mj/k 
X 

0.02 

6. 1 

-0.3 

8.6 

-4.0 

0.01 

- 

o^Cum/mj/k 

22.5 

30.3 

28. 1 

22.  t 

79.0 

12.5 

- 

P  Cm/m) 

X 

0.0 

0.0 

0.0 

■a 

■a 

0.0 

- 

P  Cm/m) 

Y 

■a 

0.6 

0.44 

■a 

■a 

0.09 

- 

X  (MPa) 

1500.0 

1260.0 

1447.0 

1026.0 

1400.0 

0.09 

400.0 

X' ( MPa  ) 

1500.0 

2500.0 

1447.0 

610.0 

235.0 

0.09 

400.0 

Y  (MPa) 

40.0 

61.0 

51.7 

31.0 

12.0 

0.09 

400.0 

Y7 (MPa) 

246.0 

202.0 

206.0 

118.0 

53.0 

0.09 

400.0 

S  (MPa) 

68.0 

67.0 

93.0 

72.0 

34.0 

0.09 

230.0 

hfl(M) 

125E-6 

125E-6 

125E-6 

125E-6 

125E-6 

0.001 

1.0 

FOR  A  SANDWICH  CORE  IN  fl  LAMINATE  LMPI=6.  IF  NECESSARY  NONZERO  CORE  PROPERT¬ 
IES  CRN  BE  GIVEN.  THESE  VALUES  ARE  TAKEN  FROM  ‘INTRODUCTION  TO  COMPOSITE 
MATERIALS'  BY  S.W.  TSAI  AND  H.T.  HAHN. 
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strength  prediction  of  hygrothermal  loads  are  not  included.  However,  one 
can  calculate  that  on  the  basis  of  mechanical  and  nonmechanical  stress  and 
strain  components  calculated  for  other  failure  theories,  or  for  these 
failure  theories. 

A  list  of  other  computer  codes  for  composite  laminate  analysis  developed 
at  the  Air  Force  Wright  Aeronautical  Laboratories,  Materials  Laboratory  is 
given  in  the  Appendix. 
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SECTION  II 


STIFFNESS  PROPERTIES 

1 .  UNIDIRECTIONAL  LAMINATE  STRESS  STRAIN  RELATIONS 

The  stiffness  of  unidirectional  composites  can  be  defined  by 
appropriate  stress  strain  relations.  These  relations  can  be  expressed 
in  terms  of  engineering  constants,  compliance  components  or  modulus 
components.  A  detailed  treatise  on  mechanics  of  composites  is  given  in 
Reference  1.  For  completeness  the  relevant  relations  are  given  in  the 
present  report.  All  the  notations  used  in  Reference  1  are  followed  in  this 
research. 

The  two  key  stress  strain  relations  for  a  unidirectional  composite 
are  given  in  Tables  2  and  3: 


TABLE  2 

ON-AXIS  STRESS  STRAIN  RELATIONS  -  COMPLIANCE 


TABLE  3 

ON-AXIS  STRESS  STRAIN  RELATIONS  -  MODULUS 


e  ‘e  e 

x  y  s 

ax 

Q  Q 

wxx  vxy 

,ay 

0  Q 

vyx  xyy 

CT 

Q 

'  s 

wss 

In  the  foregoing  matrix  multiplication  table,  each  value  in  the  first 


column  is  equal  to 

the  sum 

of  products  of  corresponding  row  elements 

with  their  column 

headings 

.  This  rule  should  be  self-evident. 

{a  ,  a  .  a.} 

S5 

{Longitudinal,  transverse,  shear) 

x  y  s 

stress  components  in  the  xy  -  plane 

{ev»  e  ,  e  } 

S 

{Longitudinal,  transverse,  shear) 

x  y  s 

stress  components  in  the  xy  -  plane 

= 

Modulus  components 

Sij 

= 

Compliance  components 

S  „ 

sz 

1/E 

XX 

X 

S  ... 

s 

1/EX/ 

yy 

y 

s 

s: 

-v  /E„ 

xy 

y  y 

s 

= 

-v /ew 

yx 

X  X 

SCc 

ss 
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1/Es 

Q 

s 

n£ 

vxx 

X 

Qyy 

s: 

mEw 

y 

^xy 

= 

mE_vw 

x  y 

Qyx 

mE  v 

y  x 

^ss 

“ 

E. 

s 

m 

= 

i/(i-Vx) 

where 

Ex 

= 

Longitudinal  Young's  modulus 

E 

- 

Transverse  Young's  modulus 

y 

e 

vx 

= 

Longitudinal  Poisson's  ratio  =  -  -f- 

fcx 

V 

_ 

ex 

Transverse  Poisson's  ratio  =  -  — 

y 

y 

Es 

= 

Longitudinal  shear  modulus  =  — 

All  the  material  constants  of  the  stress  strain  relation  shown  previously 
are  called  engineering  constants.  They  are  the  familiar  constants  used 
for  conventional  materials  with  subscript  added  to  denote  directionality 
of  properties.  Thus,  the  use  of  engineering  constants  will  often 
facilitate  the  use  of  composites  for  structural  applications.  However, 
it  has  been  found  more  convenient  to  use  compliance  and  modulus  components 
of  multidirectional  composites.  In  Equation  1  S. .  are  compliance  compo- 

■  J 

nents  and  Q.  .  are  modulus  components. 

Subroutine  MODULS  calculates  the  modulus  matrix  Q... 

*  J 

2.  TRANSFORMATION  OF  STRESS  AND  STRAIN 


The  change  of  stiffness  of  unidirectional  composites  as  a  function 
of  ply  orientation  is  a  unique  feature  of  composites.  These  orientational 
variations  of  stress  and  strain  are  the  fundamental  underlying  issues  which 
must  be  understood.  The  relations  governing  these  variations  are  called 
transformation  equations  and  are  given  in  Tables  4  and  5: 


TABLE  4 


STRESS 

TRANSFORMATION 

RELATIONS 

CT 

a1 

a 

1 

2 

6 

°x 

m2 

rt2 

2m  n 

py 

h2 

nl2 

L2mn 

as 

*-mn 

mn 

■  m2  -n2 

,  m  =  cose ,  n  =  sine 
TABLE  5 

STRAIN  TRANSFORMATION  RELATIONS 


e, 

1 

e 

2 

ed 

ex 

m2 

n2 

m^ 

*y 

n2 

ra2 

-mn 

e 

-2mn 

2mn 

m2  -n2 

s 

Where  ,  c^,  Cg  are  off-axis  stress  components,  c?x,  a  ,  ag,  are 

transformed  on-axis  stress  components  and  0  is  the  angle  of  counter¬ 
clockwise  rotation  of  the  on-axis  laminate  (Figure  1). 


Figure  1(a).  Material  Symmetry  Axis  of  a  Unidirectional  Composite 


Figure  1(b).  Off-Axis  Configuration  of  a  Unidirectional  Composite, 
Counter  Clockwise  Rotation  Is  Positive 

Similarly  e, ,  e?,  e,  are  off-axis  and  e  ,  e  ,  e  are  on-axis  strain 
components. 

Subroutine  TRS  calculates  the  transformed  stress  components. 
Subroutine  TRE  calculates  the  transformed  strain  components. 
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3.  OFF-AXIS  MODULUS 


Because  the  composite  laminates  are  made  of  off-axis  and  on-axis 
plys,  the  stiffness  of  off-axis  ply  orientation  must  be  understood.  The 
off-axis  modulus  of  a  ply  can  be  determined  in  three  steps:  (1)  the  off- 
axis  to  on-axis  strain  transformation,  (2)  the  on-axis  stress  strain 
relation,  and  (3)  the  on-axis  to  off-axis  stress  transformation.  This 
process  is  initiated  by  a  given  strain  in  Figure  2(a)  and  eventually 
leads  us  to  the  induced  stress  in  Figure  2(d).  The  three  transformations 
mentioned  above,  when  combined  together,  will  yield  the  required  off-axis 
modulus  and  off-axis  stress  strain  relations  for  arbitrary  angle  of 
rotation. 


tU 

-m* 


<v 


n 


(a) 


TS 


Q 


(b) 


!L 


(c) 


Off-Axis 

Strain 


t,Jsx,y,s 
On- Axis  On- Axis 


\ 

Id) 


Strain 


Stress 


v 


Off-Axis 

Stress 

A 


°n 


Figure  2. 


ij"  1,2,6 

Determination  of  Off-Axis  Modulus 


From  (a)  -  (b) 
From  (b)  -  (c) 
From  (c)  -  (d) 


Use  positive  angle  0  strain  transformations 
Use  the  on-axis  stress  strain  relations  in  modulus 
Use  inverse  stress  transformation. 


Following  the  above  tranformations  and  making  some  mathematical 
simplifications,  off-axis  stress  strain  relations  can  be  obtained, 
i .e. ,  Table  6. 
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TABLE  6 


OFF-AXIS  STRESS  STRAIN  RELATIONS 


£ 

T 
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2 

E 

'  6 

a 

Q 

Q 

Q 

l 

11 

12 

16 

a 

Q 

Q 

Q 

2 

21 

22 

26 

a 

Q 

Q 

Q  ' 

6 

61 

62 

66 

The  off-axis  modulus  components  are  given  in  the  matrix  form  shown  in 
Table  7. 


TABLE  7 

MODULUS  COMPONENTS  Q.  . 


1 

u 

2 

u 

3 

Qn 

,u 

1 

c,os2e 

cos4e 

Q22 

Ui 

-cos2e 

cos4@ 

Ql2 

u4 

_  -cos40 

Q66 

U5 

-  -cos4o 

Ql6 

!  If  sin2e 

sin40 

^26 

sin2e 

-sin40 

Where 

U1  =  Wxx  +  3V  *  2Qxy  +  4QsS) 

^2  "  «xx 

U3  "  «*x  +  ^yy  "  "  4!JSS)  /8 

i-  04  ‘  (Qxx  +  v  +  %  -4<W  /8 

L  U5  *  (Qxx  +  V  '  2V+ 4W  ^ 
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4.  MODULUS  AND  COMPLIANCE  OF  GENERAL  LAMINATES 

According  to  the  classical  laminated  plate  theory,  the  strain 

components  are  considered  to  vary  linearly  along  the  thickness  of  the 

laminate,  i  .e. ,  y 

e  (z)  =  +  zk 

ill 

e  (z)  =  e°  +  Zk  m 

22  2  13) 

......  .  ...  e  (z)  =  e°  +  zk 

6  6  6 

0 

where  e  (i  =  1,2,6)  are  mid  plane  strain  components  and  k.  (1=1 ,2,6)  are 
i  1 

the  corresponding  curvatures. 


A  general  laminate  may  contain  plys  of  the  same  material  or  different 
materials  (hybrid).  The  effective  modulus  of  the  composite  laminate  is  con¬ 
sidered  to  be  an  arithmetic  average  of  the  modulus  of  the  constituent  plys. 
Simple  formulas  can  be  derived  to  compute  the  effective  stress  strain 
relationship  in  terms  of  the  modulus  of  constituent  plys.  The  stress  strain 
relation  for  a  general  laminate  is  given  by: 


where  Nj ,  Ng,  Ng  are  stress  resultants  over  the  thickness  (h)  of  the 
laminate  and  M, ,  M„  and  Mc  are  corresponding  moment  resultants.  These 

12  D 


parameters  are  defined  as  follows: 

!  h/2 

N.j  =  o.  dz 
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or 


N, 


h/2 


-h/2 


h/2 


Qij £j  dz  +  Qijkjzdz  (i  &  4  =  1,2,6) 


-h/2 


or 


and 


A.  .b0.  +  B.  .k. 
J  ij  J 


M. 

1 


h/2 

* 

h/2 


zdz 


(5) 


or 


h/2 


h/2 


M.  = 


-h/2 


+  Vj*1 


•dz 


h/2 


M,  = 


B.  .er:  +  D.  -k. 
ij  J 


h/2 

-h/2 


Q.  .e .  dz 


(6) 


(7) 


h/2 

;  Qij  ej  zdz 

-h/2 


ei  =  ATa.  +  cB-  (i=x,y) 


Where  h  Is  the  thickness  of  the  laminate  and  e..  (i=x,y,s)  are  nonmechanical 
strain  components,  a.  and  g.  are  coefficients  of  thermal  expansion  and 

swelling  coefficients  respectively,  AT  is  the  temperature  difference  and 


c  is  the  moisture  contents.  In  the  foregoing  relations  the  summation  over 
the  range  of  repeated  subscript  will  be  understood.  Figure  3  shows  a 
reference  frame  for  numbering  the  layers.  The  detailed  explanation  to 
these  equations  Is  given  In  Reference  1, 

Subroutine  NMSN  calculates  the  nonmechanical  strain  components  for  each 
layer. 

Subroutine  MOLS  calculates  the  effective  inplane  modulus  matrix  A,  effective 
flexural  modulus  D,  coupling  matrix  B,  nonmechanical  stress,  and  moment 
resultants , 


The  constitutive  equations  (Equation  4)  can  be  rewritten  in  the  following 
form : 


(8) 


e° 

2  B 

N  +  NN 

w-i 

X 

k 

w 

-T  * 
_g  $ 

M  +  MN 

, 

~  J 

where 


-1 

Inplane  compliance  =  a  +  aB  (D-BaB)  Ba 


VTiA  V/V 


Coupling  compliance  =  -aB  (D-BaB) 

-1 

Flexural  compliance  =  (D-BaB) 

-1 


-1 


It  has  been  shown  (Reference  2),  that  the  use  of  normalized  quantities  in 
the  stress  strain  relations  is  very  helpful  in  understanding  the  variation 
in  effective  material  properties  of  composite  laminates  due  to  the 
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f\>{=r 


In  these  equations,  the  multiplying  factors  3  and  1/3  are  consequent  to 
the  normalization  factors  defined  in  Equation  9,  The  well  known  Kirchhoff's 
assumption  in  classical  plate  theory  of  the  linear  strain  distribution  along 
the  ply  thickness,  can  be  stated  as: 


£.=£.+  z*k? 
1  1  1 


(i  =  1,2,6) 


where  z*  =  2z/h  and  -1<z*<1 


Thus,  all  the  equations  required  for  computing  the  effective  modulus  matrix 
are  given. 

5.  INPLANE  AND  FLEXURAL  ENGINEERING  CONSTANTS 
a.  Inplane 

The  inplane  stress  strain  relation  can  be  written  in  the 
fol 1 owi ng  form : 

N*  =  A*  e* 

or  -1 

e  =  a*  N*  ,  a*  =  (A*) 

•V  A,  -v  ^  ^ 

The  typical  effective  engineering  constants  are  given  by: 


Inplane  longitudinal  modulus  =  EQ  =  1/a* 

1  n 


Inplane  transverse  modulus  =  E  =  1/a* 

2  22 


Inplane  shear  modulus  =  E 


1/a* 

66 


Inplane  Poison's  ratio  =  v 


15 


b.  Flexural 


The  moment-curvature  relation  can  be  written  as: 
M*  =  D*  k* 

or  -1 

k*  =  d*  M*  ,  d*  =  (D*) 


The  typical  effective  flexural  engineering 
Longitudinal  Young's  modulus  = 

Transverse  Young's  modulus  =  eJ  = 

Shear  modulus  =  =  1/d* 

6  66 

Poisson's  ratio  =  vf  =  -d*  /d* 

21  12  11 


constants  are  given  by; 

=  1/d* 

1 1 

1/d* 

22 


In  unsymmetric  laminates,  the  effective  stress  strain  relations  and  moment 
curvature  relations  will  not  be  uncoupled.  Therefore,  in  those  situations 
the  effective  inplane  and  flexural  engineering  constants  may  not  give 
much  meaningful  information. 

6.  MATERIAL  PROPERTIES  FOR  FINITE  ELEMENT  ANALYSIS 

The  theory  used  in  this  report  is  based  upon  Kirchhoff's  classical 
plate  theory.  For  the  stress  analysis  of  composite  laminates  by  finite 
element  general  purpose  computer  code,  NASTRAN,  the  effective  modulus 
matrix  A.  ./h  is  used.  In  the  NASTRAN  material  property  input  data  card 

*  J 

MT2,  the  numerical  values  are  given  such  that  the  notation  6^.  in  the 
NASTRAN  manual  represents  the  effective  modulus  A^./h  in  this  report. 
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SECTION  III 
FAILURE  THEORIES 


Figure  4  shows  the  failure  criteria  survey  response  obtained  by  the  AIAA 
composite  materials  subcommittee  [3].  On  the  basis  of  this  response  we  have 
chosen  maximum  stress,  maximum  strain  and  quadratic  polynomial  (Tsai-Wu  [4], 
Chamis  [5],  Hoffman  [6],  and  Hill  [7])  failure  criteria  for  a  comparative 
treatment  in  this  work.  Reference  8  gives  an  extensive  survey  of  failure 
theories.  The  theories  considered  in  this  investigation  are  given  in  {9,10,11]. 
For  completeness,  these  theories  are  restated  here 

Maximum  Stress  Maximum  Strain 


Y 


a 


s 


<  S 


Failure  occurs  when  one  of  the  equalities  is  met. 

Quadratic  Polynomial 

Four  quadratic  polynomial  failure  criteria  given  in  Table  8,  have  been 


Figure  4.  AIAA  Failure  Criteria  Survey. 
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TABLE  8 

QUADRATIC  POLYNOMIAL  FAILURE  CRITERION 


i  t  e  r  i  a 

Par  ame  ter'^>>>*^ 

Tsai-Wu 

Chamis 

Hoffman 

Hill 

Equation 
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xy 
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xy 

Material 
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1  /yy  ' 

"  2  vXX' 

1  Y 

2  X 

1 

1 

1 

1 

SS 

„2 

2 

.2 

2 

S 

S 

S 

S 

1  1 

1  1 

X 

X  X’ 

X  X* 

.  . 

.  ..  ...  _ 

1  1 

1  1 

y 

Y  Y’ 

Y  Y* 

+ 

The  values  within  the  parenthesis  are  used  when  the  corresponding  stress 
component  is  compressive. 

+  + 

In  the  Tsai-Wu  criterion  F*y  is  taken  to  be  -1/2. 

+*f+  ... 

In  Chamis  criterion  F*  varies  from  material  to  material.  For  T3OO/5208 


X  =  Longitudinal  tensile  strength 
X'  =  Longitudinal  compressive  strength 
Y  =  Transverse  tensile  strength 
Y'  =  Transverse  compressive  strength 
S  =  Longitudinal  shear  strength 


’FXX 

Fxy 
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rF  _ 
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F  = 

Fxy 

Fyy 
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,  F  = 

F 
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,  O  = 
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_0 
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Fss- 

0 

-  cs  _ 
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A  superscript  T  denotes  transpose  of  the  matrix  and  the  subscripts  x,  y 
and  s  denote,  respectively,  longitudinal,  transverse  and  shear  directions, 

E  denotes  the  Young's  modulus  of  the  ply.  The  quadratic  polynomial  failure 
criterion,  through  the  use  of  the  stress  strain  relations,  can  be  expressed 
in  strain  space  as  follows: 

eTGe  +  6Te  =  1  (13a) 

-T 

where  the  elements  of  matrix  G  and  vector  G  are  dependent  upon  the  ply 
modulus  matrix  Q. and  F. .  and  F.  of  Table  8  and 

'  J  ■  J  1 


For  the  computation  of  the  failure  envelopes  in  the  qr  space  the  following 
definitions  are  used: 


=  (  G-|  +  ££)/2 

Po  =  (al  +  aZ^2 

=  ( e-j  -  e2)/2 

<l0  =  <°1  -  °2)/2 

=  V2 

ro  =  °6 

where  e  and  a .  are  stress  and  strain  components  in  the  laminate  axes. 

Chamis,  Hoffman  and  Hill  criteria  can  be  deduced  from  Tsai-Wu  failure 
criterion  by  using  different  values  of  F*  and  appropriate  X,  X',  Y  and  Y' 
as  given  in  Table  8. 

According  to  the  quadratic  polynomial  failure  criterion  the  failure 
surface  in  strain  space,  equation  13a,  is  represented  by  the  following 
quadratic  equation: 

G.jC.Ej  +  Giei  =  1  ( i ,3  =  x,y,s)  (13) 

where  nonzero  coefficients  G. .  in  the  foregoing  equation  are  given  by 

'  J 
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G  =  F  Q2  +  2F  Q  Q  +F  Q2 

xx  xxyxx  xyyxxyxy  yy  xy 

G  =  F  Q2  +2F  Q  Q  +F  O2 
yy  xxyxy  xyyxyyyy  yyyyy 

G  =  F  Q  Q  +F  [Q  Q  +Q2  ]  +  F  Q  Q 
xy  xxyxxyxy  xyLyxxyyy  yxyJ  yy  xyvyy 

G  =  F  Q2  =  (Q  /S)2 

ss  sswss  vyss  7 


(14) 


G  =  F  Q  +F  Q 
x  xyxx  yyxy 

G  =  F  Q  +F  Q 
y  xyxy  yyyy 


Subroutine  FAILCO  computes  the  coefficients  G^  and  G^  in  Equation  13, 
Equation  13  can  be  written  in  the  following  expanded  form: 


2  2  2 

jG  £  +2G  ££  +  G  £  +  G  £ 
xx  x  xy  x  y  yy  y  ss  s 

+G  £  +G  £  =  1 

y  ^  y  y 


(15) 


With  the  knowledge  of  applied  stress  or  strain,  the  above  relation  will 
enable  the  designer  to  predict  if  the  structure  is  going  to  survive  or 
fail.  This  relation  becomes  easier  to  use  if  we  introduce  a  strength  ratio 
parameter  R  defined  by: 


M 


M 


-  -  =  Dp 

allowed  imposed 


(16) 


Here,  the  superscript  M  denotes  the  mechanical  strain.  The  strains  that 
must  satisfy  the  failure  criterion  are: 


Gi  Jallowed 
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- 

"  e.J 


£  +  £ 

X  X 

M  ,  N 

£  +  E 

y  y 


i Jallowed 
N 


+  £i 


e. 

i 


-  e 


x 

e 


(17) 


M  N 
es  +  es 


The  superscript  N  denotes  nonmechanical  strain,  ex  and  e^  are  longitudinal 

and  transverse  strain  components  in  the  ply  axis.  Using  Equations  16  and 
.17,  Equation  15  is  given  in  the  following  form: 
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(18) 


G.j  (R  e?  +  -  e. )  (Re j  +  e J  -  ) 

+  G.  (R  £|+  e!  *  e.)  =  1 


With  algebraic  rearrangements  this  equation  reduces  to: 


aR2  +  bR  +  c  =  0 


(19) 


Where  a,  b  and  c  are  coefficients  of  R2,  R  and  1  in  Equation  18  and  are 
dependent  upon  the  material  properties  of  the  laminate  and  applied  loadings. 
The  roots  of  this  equation  give  the  strength  ratios  R  and  R~,  where  R" 
denotes  the  reversed  Toad  strength  ratio. 

Subroutine  FSFTY  computes  the  strength  ratios  R  and  R”,  However, 
in  the  print  output  only  R  has  been  given.  A  parameter  R/h  that  corresponds 
to  the  strength  of  the  ply  (in  case  of  a  unit  applied  stress  resultant 
N. ,  i=l,2,6)  has  also  been  given  in  the  print  output.  Since  on  variation 
of  ply  number  or  orientation  in  a  laminate  the  total  laminate  thickness 
hay  vary,  thus  the  introduction  of  R/h  will  be  helpful  in  understanding 

i 

ithe  strength  of  different  laminates, 

j.  .  . 

For  Chamis'  criterion,  when  nonmechanical  loads  are  included, 
the  signs  of  mechanical  stress  components  govern  the  choice  of  F 

XX 

and  Fyy  of  table  8.  This  can  be  changed  to  incorporate  the  resultant 
;ply  strains  in  subroutines  FSFTY  and  STRNG. 
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FLOW  DIAGRAMS  FOR  DIFFERENT  SUBROUTINES 
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Subroutine  MODULS^  Calculates  ply  on  axis  modulus 


Subroutine  TRS: 
Stress  transformation 


Subroutine  TRE: 
Strain  transformation 
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Subroutine  MODCM:  Ply  off  axis  modulus 


Eng'g 

Const. 

e 

MOGULS 

> 

r 

Ui 

Table  7 

i 

°u 

(0) 

MW 


Subroutine  NMSN;  Ply  nonmechanical  strain 


Flaw  Chart  for  fnplane  nonmechanical  strain 
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Subroutine  MOLS:  Effective  modulus  and  nonmechanical  loads 


NNL  s  Number  of  layers  in  the  laminate. 


Program  Main 


APPENDIX  A 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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INTEGER  TITL(4)*CASE,PARAM*PARAM1,PARAM2,PARAM3,MATRIAL,UNITS 
DIMENSION  PL  NM ( 4  0  >  *  T (3) 

DIMENSION  ZERO (3  *  3 ) *DUM (3  *3 ) *DUM 1 ( 3  * 3 ) 

DIMENSION  U(5)*DU7(3*3)*HS(41)*AN(3)*AM(3) 

DIMENSION  DU1(3*3)*DU2(3,3)*DU3(3*3>*DU4(3,3)*DU5(3*3) 

DIMENSION  ES11C40)*ES22(40) *VS12(40) *GS12(40) , SB1 ( 3 * 3 ) * SB2 ( 3* 3 > 
DIMENSION  SXC(40) *SXL(40) *SYL(40) *SYC(40) *SXLT (40) 

DIMENSION  SB3(3*3)*HLC(40)*HT(40)*TL(40)*TT(40)*P(3)*SH(8) 
DIMENSION  EX(8)*EY(8)*VX(8)*ES(8)  * ALFX ( 8 ) » ALFY ( 8 ) , BT A X ( 8 ) 

DIMENSION  BTAY(8)*X(8) *  XD(8)*Y(8)*  YD ( 8 ) * S ( 8) * LMP I ( 4 0 ) 

DIMENSION  Q(3*3),A(3*3),HL(40)*TH(40)*SFL(40)*QQ(3) 

DIMENSION  N AMES ( 17) 

DIMENSION  NE UN ( 3 ) 

COMMON  /SF/SFL*COUNT 

COMMON  /BLH/  H*HI 

COMMON  /BLM/  P AR AM2 *P A R AM 3 

COMMON /BLK/HLC,HT*TL*TT 

C0MM0N/$BS/SB1*SB2*SB3*P»QQ*DT*C 

COMMON  /UNIT/  IUNIT 

COMMON  /BLK 1 /  SXL , SX C * S YL * S YC * S XL T , HL 
COMMON  /TIT/  TITL 

COMMON  /BKK  9/  PXX * PX X  I , PY Y * P Y Y I  * YUL *  I T 
COMMON  /BLK2/  E SI 1  * E S22 *  VS 12 * GS1 2 * NNL 

COMMON  /DTA/  EX ,E Y * V X *ES * ALFX * ALF Y, B T AX  * BTA Y *X * XD « Y *  YD  * S * SH * NNH 

COMMON  /MPIL/  LMP I 

COMMON  /3LDU/  DU2 *DU4 , DU5 * DU7 

NAMELIST  /LAMDATA/NNM,EX,EY*VX,ES*ALFX* ALFY * BT AX , BT A Y ,X , XD * Y , YD , S , 
*SH 


NAMELIST/LAYERS/NNL*LMPI*TH*PLNM*IUNIT*DT*C*NLDCN 

NAMELIST/LAYER/NML*TH,PLNM,DT*C,NLDCN 

NAMELIST  /STRESS/  AN  *  AM 

DATA  NAMES/7HINITIAL*6HTHEEND*9HTRANSF0RM*6HSTRESS*6HSTRAIN* 
*7HMODULUS*10HCOMPLIANCE*6HMODCOM*8HLAMINATE,7HINPLANE *4HPURE* 
*6HHYBR I D,7HGENERAL*8H STRENGTH, 4H  *10HSTRNGTHPLT*7HENGCPLT/ 
DATA  NEUN/8HNEWMTRLS*2HSI*7HENGLISH/ 
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COMPOSITE  LAMINATE  ANALYSIS  PROGRAM 


*  DEVELOPED  BY 

* 


*  SOM  R.  SONI 

*  UNIVERSITY  OF  DAYTON 

*  RESEARCH  INSTITUTE 

*  300  COLLEGE  PK. 

*  DAYTON  OHIO  45469 

*  TEL.  #  (513)  255-  6809 

* 


FORMAT (5X** 

*- - +) 

URITE  13 
WRITE  1001 
WRITE  13 
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WRITE  1003 
WRITE  13 

1001  FORMAT (  5  X  »  *  A  COMPOSITTE  LAMINATE  ANALYSIS  PROGRAM*) 
CALL  LMDT 

1002  FORMAT <5X**DEVEL0PED  BY**/* 

15X  *  *SOM  R.  SONI**/* 

25X**UNIVERSITY  OF  DAYTON  RESEARCH  INSTITUTE**/, 

35X  »  *  D AY  TON  OHIO  45469*,/, 

45X  » *FOR  QUESTIONS  CALL  513-255-6809*,/) 

1003  F0RMATC5X**MATERI AL  PROPERTY  DATA*) 

WRITE  13 

DO  303  1=1,3 
DO  303  J  =  l,3 
ZERO(I,J)=0. 

303  CONTINUE 

706  FORMAT ( 7F1 0 • 3 ) 

799  CONTINUE 
I  UN I T  =  2 
WRITE  711 
711  FORMAT(lHl) 

READ  700*<TITL(I) *1=1,4) 


IFCTITL(l) 

.EQ. 

NAMES(2) ) 

STOP 

IF(TITLCl) 

.EQ. 

NEWN(l) ) 

READ 

LAMDATA 

IF(TITL(1) 

•  EQ. 

NEWN(l) ) 

GO  TO 

799 

PRINT  LAMDATA 

FORMAT (6A10 ) 

IFCTITLC1) 

•  EQ. 

N AMES  <  2 ) ) 

GO  TO  902 

IF(TITLCl) 

•  EQ. 

NAMESC9) 

•  AND. 

TITL (2) 

.EQ. 

NAMESdO)  ) 

GO 

TO 

722 

IF(TITL(1) 

•  E  Q  « 

NAMES ( 9 ) 

.AND. 

T I TL ( 2  ) 

.  E  Q  . 

NAMESC13) ) 

GO 

TO 

'  723 

IF(TITLCl) 

.EQ. 

NAMES ( 3  > 

.AND. 

TITL (2) 

.EQ. 

NAMES (4 ) > 

GO 

TO 

724 

IF(TITL<1) 

.EQ. 

NAMES (3) 

•  AND. 

T I TL ( 2  ) 

.  E  Q  . 

NAMESC5) ) 

GO 

TO 

725 

IF(TITL<1) 

.EQ. 

NAMES ( 3 ) 

.AND. 

T I TL ( 2  ) 

•  EQ. 

NAMESC6) ) 

GO 

TO 

726 

IF (TITL(l) 

.EQ. 

NAME  S  <  3 ) 

.AND. 

TITL (2) 

•  EQ. 

NAMES(7) ) 

GO 

TO 

726 

IF  (TITLd) 

.EQ. 

NAME  S ( 3 ) 

•  AND. 

TITLC2) 

•  EQ. 

NAMESC8)) 

GO 

TO 

726 

724  CONTINUE 

READ  701,  (AN (I), 1=1, 3),  THETA 

701  FORMAT ( 4F10. 3 ) 

WRITE  702, THETA 

702  FORMAT (5X,*STRESS  TRANSFORMATION  THRU  * ,F7 .2 , 2 X , *DEGREE S* , / , 1  OX , 

.*  SIGMA1  SIGMA2  SIGMA6*) 

WRITE  13 

WRITE  703,  <AN<I)*I=1,3) 

703  FORMAT  <5X,*GIVEN  *,3F10.3) 

CALL  TRS< AN, T, THETA) 

WRITE  704, (T( I) ,1=1,3) 

704  F0RMAT(5X,*TRANSF0RMED  **3F10.3) 

WRITE  13 

GO  TO  799 

725  CONTINUE 

READ  701,  <AN(I),I=1,3),  THETA 
WRITE  705, THETA 

705  F0RMAT(5X**STRAIN  TRANSFORMATION  THRU  * , F7 . 2 , 2 X , * DEGREES* , / , 1 0 X , 

.*  EPSLN1  EPSLN2  EPSLN6*) 

WRITE  13 

WRITE  703,  (ANCI) ,1=1,3) 

CALL  TRE( AN, T, THETA) 

WRITE  704, (TCI), 1=1, 3) 

WRITE  13 
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GO  TO  799 

726  CONTINUE 

READ  700, MATRIAL, UNITS 
IFCUNITS  .NE.  NEWN ( 3  )  )  IUNIT=1 
CALL  MINP(MATRIALtLI) 

IF (LI  .EQ.  8)  READ  7 06 , E X ( L I ) * E Y ( LI ) * VX ( L I ) ,ES < L I )  , TH ( 1  ) 

I F  < L I  .LT.  8)  READ  706, TH(1> 

NNL  =  1 

LHP  I ( 1) =L I 
PLNM  <!>  =  !• 

CALL  MTDM (PLNM,PM  *GM  *PMI , TP  I ) 

CALL  MODULS(ESlld),ES22(l)*VS12(l),GS12(l)*Q,U) 

CALL  M0DCM(TH( 1) ,U,SB1) 

CALL  MTAD (ZERO,  SB1,SB2,PMI) 

CALL  MTAD(ZER0»Q*SB3,PMI) 

IF(TITL(2)  .EQ.  NAMES(7)>  GO  TO  727 
WRITE  707  *  TH ( 1 ) 

707  F0RMAT(5X,*M0DULUS  TRANSFORMATION  THRU  *,F7.2,*  DEGREE  ANGLE*) 

IF ( IUNI T  .NE.  1)  WRITE  713 

IFdUNIT  .EQ.  1)  WRITE  712 

712  FORMAT  (30X,*(GPA)*) 

713  FORMAT (30X,*(l.F+06  PSI)*) 

WRITE  13 

WRITE  709 

709  FORMAT (5Xi*M0DULUS  OF  THE  MATER  I AL* *  1  OX  * *TR ANS FORMED  MODULUS*) 
WRITE  13 

CALL  WITE1(SB3,SB2) 

WRITE  13 

IF(TITL(2)  .EQ.  NAMESC6))  GO  TO  799 

727  CONTINUE 

CALL  INVRS(SB3»DUM) 

CALL  INVRS(SB2,DUM1) 

CALL  MTAD(ZERO,DUM,DU1,1000.  ) 

CALL  MTAD(ZERO,DUM1*DU2»100Q.) 

WRITE  708,  TH ( 1 ) 

708  FORMAT (5X,*C0MPLIANCE  TRANSFORMATION  THRU  *,F7.2,*  DEGREE  ANGLE* 

.) 

IF(IUNIT  .EQ.  1)  WRITE  714 
IFdUNIT  .NE.  1)  WRITE  715 

714  FORMAT(30X,* (1/TPA)*) 

715  FORMAT(30X,*(1.0E-09/PSI)*> 

WRITE  13 

WRITE  710 

710  FORMAT (5X,*C0MPLIANCE  OF  MATER  I AL*, 13X, TRANSFORMED  COMPLIANCE*) 
WRITE  13 

CALL  WITE1(DU1,DU2) 

WRITE  13 
GO  TO  799 

722  READ  700, CASE, PARAM,PARAM1,PARAM2,P ARAM3 
TITLC4)=CASE 

I F (CASE  .EQ.  NAMES(12))  GO  TO  731 

IF (CASE  .NE.  NAMES ( 1 1) )  STOP  "ERROR  IN  MIND" 

CALL  LMNT(NNL,NLDCN,TH,LMPI ,DT,C,PLNM) 

CALL  MTDM(PLNM,PM,GM,PMI,TPI> 

CALL  MOLS (TH , HL ,HS ) 

CALL  INVRS(SB1,A) 

CALL  MTAD(ZERO,A,DUM,H) 

CALL  MTAD(ZER0,SB1,DUM1,HI) 

CALL  N0RM(PMI,ZER0,DU1,DU2,DU3,DU4,PARAM) 

CALL  N0RM1(TPI,A,PARAM1) 
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IF (PARAM3  .EQ.  NAMES<17))  CALL  PLTEN < TH *HL *HS *  PM  I  * ZER 0  * PLNM > 
GO  TO  799 

731  CONTINUE 

CALL  LMNT  (NNL*NLDCN»TH*LMPI *DT  »C*PLNM) 

CALL  MTDMCPLNH*PHf6M*PMI*TPI) 

CALL  MOLS(TH,HL,HS) 

CALL  INVRS(SB1*A) 

CALL  MTAD(ZERO*A*DUM»H) 

CALL  MTAD(ZER0,SB1,DUM1,HI> 

CALL  NORM (PM  I f ZER 0 , DU  1  * DU2 * DU3 ♦ DU4 1 P ARAM) 

CALL  N0RM1(TPI*A,PARAM1) 

IF (PARAM2  .EQ.  NAMES<14)>  CALL  ST RNG ( TH « HS »  PM  * NLDCN ) 

IF ( PAR AM2  .EQ.  NAME  S ( 16 ) )  CALL  S TPLT < TH , HS *  PM , PLNM > 

IF(PAR AM3  .EQ.  NAMES<17>)  CALL  PLTEN ( TH , HL * HS ,PM I  , Z ERO ,PLNM ) 
GO  TO  799 

723  READ  70 0« CASE *P AR AM f P AR AM  1 ,P AR AM  2 
T  I TL ( 4 ) =  C  ASE 

I F  ( C ASE  .EQ.  NAMES(12)>  GO  TO  732 
I F  <  C ASE  .NE.  NAMES(ll))  STOP  "ERROR  IN  MIND" 

CALL  LMNT ( NNL * NLDCN ♦ TH *LMP I ,DT*C*PLNM> 

CALL  MTDM(PLNM*PM*GM»PMI»TPI> 

CALL  MOLS(TH,HL*HS) 

CALL  INVRS(SB1*A) 

CALL  MTAD(ZERO,A,DUM.H> 

CALL  MTAD(ZER0»SB1*DUM1*HI) 

CALL  NORM(PMItZERO.DUl*DU2*DU3*DU4*PARAM) 

CALL  N0RM1(TPI»A*PARAM1  ) 

I F ( P  AR  AM  2  .EQ.  NAMESC14))  CALL  STRN G ( TH , HS t PM * NLDCN ) 

I F ( P  AR  AM  2  .EQ.  NAMES<16>)  CALL  STPL T ( TH * HS *  PM , PLNM ) 

IF ( PARAM3  .EQ.  NAMES(17)>  CALL  P L TE N ( TH , HL * HS , PM  I  * ZE RO, PLNM > 
GO  TO  799 

732  CONTINUE 

CALL  LMNT(NNL»NLDCN»TH»LMPI*DT*C*PLNM) 

CALL  MTDM(PLNM*PM*GM*PHI*TPI) 

CALL  MOLS(TH»HL*HS) 

CALL  INVRS(SB1,A) 

CALL  MTAD(ZERO, A,DUM,H) 

CALL  MTAD(ZER0»SB1»DUM1»HI  ) 

CALL  N0RM<PMI*ZER0,DU1*DU2* DU3 , DU4 * P AR AM ) 

CALL  N0RM1(TPI,A,PARAM1) 

IF (PARAM2  .EQ.  NAMES(14)>  CALL  STRNG ( TH ,HS *PM « NLDCN ) 

IF (PARAM2  .EQ.  NAMES(16>>  CALL  STPL T C TH * HS *  PM * PLNM ) 

I F ( P  AR AM  3  .EQ.  NAMES<17>)  CALL  PLTE N ( TH , HL* HS t PM  I , Z ER 0 « PLNM ) 
GO  TO  799 
902  CONTINUE 
WRITE  13 
STOP 
END 
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SUBROUTINE  AD JUST  < I CND , SFLN , SPLC , AN ) 

INTEGER  CRITRIA*SPACE»PLANE 

DIMENSION  STX1(200)»STY1(200)*SCX1(200)*SCY1(200) 
DIMENSION  STX2C 200) »STY2 (200 ),SCX2< 200 ),SCY2 (200) 
DIMENSION  STX3(200) fSTY3(200  > ,SCX3(2  00) ♦SCY3C20  0) 
DIMENSION  SFLN(40)fSPLC(A0) , AN ( 3 > ,N AME < 2 > 

COMMON  /S T 1 /  STX1»STY1*SCX1»SCY1»NSF1 
COMMON  /ST2/  STX2 *STY2 * SCX2 * SCY2 *NSF2 
COMMON  /ST3/  STX3 ,ST Y 3 ♦ SCX 3  *  SC Y 3  * NS F 3 
COMMON  /CRIT/  CR I TR I  A * SPAC E , PL ANE » FS 1 2 
DATA  NAME/2HQR,6H STRESS/ 

IF ( PLANE  .EQ.  NAME(1>)  GO  TO  1 
ST XI <ICND)  =  SFLN(NSF1 >  * AN  < 1 > 
STY1(ICND)=SFLN(NSF1)*AN(2> 

SC XI <ICND)=SPLC<NSF1)*AN< 1 ) 
SCY1(ICND)=SPLC<NSF1)*AN<2> 

STX2< ICND)=SFLN(NSF2)*AN(1) 
STY2(ICND)=SFLN(NSF2)*AN<2) 
SCX2(ICND)=SPLC(NSF2)*AN(1) 
SCY2<ICND)=SPLC(NSF2)*AN<2) 
STX3<ICND)=SFLN(NSF3)*AN(1) 
STY3(ICND)=SFLN(NSF3)*AN<2) 
SCX3(ICND)=SPLC(NSF3)*AN(1> 
SCY3(ICND)=SPLC(NSF3)*AN<2) 

RETURN 

CONTINUE 

SQ  2  =  SQR  T ( 2  •  0  ) 

SQ3=SQ2 

IF (SPACE  .EQ.  NAME(2) )  SQ3=1./SQ2 

STX1(ICND)=SQ2*SFLN(NSF1)*ANC1) 

STX2(ICND)=SQ2*SFLN(NSF2)*ANC1) 

STX3<ICND)=SQ2*SFLN<NSF3)*AN<1> 

SCX1(ICND)=SQ2*SPLC<NSF1)*AN<1) 

SCX2(ICND)  =  SQ2*SPLC(NSF2) *AN(1 ) 
SCX3(ICND)=SQ2*SPLC(NSF3)*AN(1) 

STY1 <ICND)=SFLN(NSF1)*AN(3)/SQ3 

STY2(ICND)=SFLN(NSF2)*AN(3)/SQ3 

STY3<ICND)=SFLN(NSF3)*AN(3)/SQ3 

SCY1(ICND)=SPLC(NSF1)*AN(3)/SQ3 

SCY2(ICND)=SPLC(NSF2)*AN(3)/SQ3 

SCY3(ICND)=SPLCCNSF3)*ANC3)/SQ3 

RETURN 

END 


SUBROUTINE  AM AX < XX , NN, AMX ) 
DIMENSION  XXtlO) 

AMX=XX( 1) 

DO  1  1=1, NN 

IF(XXtl)  .GE.  AMX)  AMX=XXCI) 

CONTINUE 

RETURN 

END 


SUBROUTINE  AM  I N < XX , NN ♦ AMN , I K > 
DIMENSION  XX ( 2 0 ) 

AMN=XXC1) 

DO  1  1=1, NN 

IF(XXCI)  .LE.  AMN)  GO  TO  2 
GO  TO  1 

2  AMN=XX(I) 

IK  =  I 

1  CONTINUE 

RETURN 
END 


SUBROUTINE  COEF(G,GB,GS,GSB,THET) 

DIMENSION  G<3*3),GB<3)»GS(3,3)*GSB<3)*U<5) 
CALL  US ( G , U ) 

CALL  MODCM(THET*U»GS) 

Hl=(GB(l)+GB(2))/2. 

H2=(GB(l)-GB(2))/2. 

TH2=2.*THET 
C  =  COSM  <  TH2 ) 

S=SINM(TH2) 

GSB ( 1 ) =H1+H2*C 

GSB  <2>=H1-H2*C 

GSBC3)=H2*S 

RETURN 

END 
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SUBROUTINE  F A ILCO C E * Q , X ,XD , Y * YD  * XLT * G ,GB ) 

INTEGER  CRITRIA, SPACE, PLANE 
DIMENSION  E(3),E1(3) 

DIMENSION  Q(3,3)*G(3»3),F(3,3),GB(3),FB(3)*C(3*3) 

DIMENSION  NAME(fe) 

COMMON  /CRIT/  CRI TR I  A  *  SPACE , PL ANE ,F S 1 2 

DATA  NAME  /7HTSAI  WU , 7HH0FFM AN  * 4HHI LL , 6H CH AM  I S , 6HSTR A IN , 6H STRE SS/ 

XL=X 

XC  =  XD 

YL  =  Y 

YC  =  YD 

IFICRITRIA  .EQ.  N  AME  < 1 )  .OR.  CRITRIA  .EQ.  NAME(2)>  GO  TO  1 
IFCCRITRIA  .EQ.  NAME<4)>  GO  TO  2 
XC  =  X 
YC  =  Y 
GO  TO  1 

2  CALL  MVM(Q»E*E1) 

IF(E1(1)  . L  T • 0 •  )  XL  =  XD 
IF(E1<1)  .GT  .  0  .  )  XC  =  X 
IF ( E 1 < 2 )  »L T  •  0  »  )  YL=YD 
IF(E1C2)  • G  T • 0 • )  YC=Y 
1  CONTINUE 

FC1,1)=1.0/(XL*XC) 

F(2,2)=1.0/(YL*YC) 

C  FC1,2)=0. 

C0E=SQRT(F(1,1)*F(2,2)) 

F< 1,2)=FS12*C0E 

IF(CRITRIA.EQ.  NAME ( 2 )  .OR. CRITRIA  .EQ.  NAME(3)>  F< 1  *  2  >  =-. 5*F (1,1) 

IF (CRITRIA  .EQ.  NAME (4 ) )  F < 1 ,2 )=FS1 2*COE 

FC3,3)=1.0/(XLT*XLT) 

FB<1)=(XC-XL)*FC1,1) 

FB<2)=(YC-YL)*F<2,2> 

FB<3)=0.0 

F(2,1)=F(1,2) 

F  <  3 , 1 ) =0. 0 
F(3,2)=0.0 

F (1,3) =0.0 
F(2,3)=0.0 
CALL  MATM(Q,F,C) 

CALL  MATM(C,Q,G) 

CALL  MVM( Q,FB,GB) 

RETURN 

END 


SUBROUTINE  F SFTY ( ET ,EN ,G* GB , SM , SMN ) 

DIMENSION  ET<3),EN(3),G(3,3),GB(3),TS(3),T(3) 
CALL  MVM (G,ET»T) 

CALL  MVM ( G,EN,TS) 

AC  =VVM ( T , ET  ) 

BC=  VVM(T,EN)+VVM(TS,ET)+VVM(ET,GB) 

CC  =  VVM(EN,TS)  +  VVM(GB,EN)-1. 0 
CALL  ROOTS< A C ,BC ,CC , SM , SMN ) 

RETURN 

END 
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SUBROUTINE  INVRSCQt A) 

DIMENSION  Q(3,3) ♦ A<3,3) 

DET=Q(1*1>*QC2*2>*QC3,3)+2.0*Q{1,2>*Q<1,3>*Q(2,3)-G<2,2>*Q<1,3)*« 
*-Q(l,l)*Q<2,3>**2-Q<3,3)*QCl,2>**2 
A<1,1>=CQ<2,2)*Q<3,3>-Q(2,3)**2>/DET 
A(3*3)=CQCl*l)*Q(2t2>-Q<l»2)**2)/DEI 
A(2,2)=(Q<1»1>*G<3,3>-Q<1,3)**2>/DET 
A(1,3)=(Q(1»2}*Q(2»3)-Q(2,2)*Q(1,3) )/DET 
A<1,2>  =  <Q(1,3>*G(2,3>-Q<1,2>*Q<3,3>  >/DET 
A<2,3)=<Q<1»2)*Q<1,3)-Q<1,1>*Q<2,3>>/DET 
DO  1  1  =  1,3 
DO  1  J=  1 »  3 
A( J,I)  =  A< I* J) 

1  CONTINUE 

RETURN 
END 


C 

C 

c 

c 


c 

c 

c 

c 

c 


SUBROUTINE  LMDT 

DIMENSION  EX(8),EY<8) ,VX<8) ,ES<8>  , ALF X ( 8 > , ALF Y C 8 ) , BT AX < 8 > 

DIMENSION  BTAY<8) ,X(8)  ,  XD(8),Y(8),  YD ( 8  )  ,S (8) ,SH (8 ) 

COMMON  /DTA/  EX * E Y • V X , ES , A LFX , ALF Y, BT A X , BTA Y « X , XD * Y , YD , S « SH « NNM 
DATA  EX/145. *2 04. *138. ♦ 38 .6 »76 . ♦ 0 . ♦ 69 . , 0 . / 

DATA  EX/ 1 81., 2 04. *  13 8., 38. 6 *76. *0.* 69. *0./ 

DATA  EY /ll. 0,18. 5, 8. 96, 8. 27, 5. 5*0., 69. *0./ 

DATA  EY/10.3$18. 5*8.96*8. 27*5.5»0.»69.*0./ 

DATA  VX/. 25*. 23, .3, .26, .34,0. *.3,0./ 

DATA  VX/.28, . 23  * . 3, . 26  * . 34 , 0 . * . 3 *  0 . / 

DATA  ES/5. 50, 5. 59 *7. 1,4. 14, 2. 3,0., 26. 5,0./ 

DATA  ES/ 7. 17, 5. 59, 7. 1,4. 14, 2. 3,0., 26. 5,0./ 

DATA  ALFX/.02,6.1,-.3,8.6,-4.0f.O,0.,0./ 

DATA  ALFY/22. 5, 30. 3, 28. 1,22. 1,79.0,. 0,12. 5,0./ 

DATA  BTAX/0.,0.,0.,0.,0.,0.,0.,0./ 

DATA  BTAY/0.6,0.6,.44,.6,.6,.0,0.,0./ 

DATA  X /I 448., 1260. ,144 7., 1062., 1400.,. 0, 400. ,0./ 

DATA  X/l 500., 1260. ,144 7., 1  062., 1400.,. 0, 400., 0./ 

DATA  XD/1448. ,2500.,1 447. ,610., 235. ,.0,400.,  0./ 

DATA  XD/1500.,2500., 1447., 610., 235. *.0, 400., 0./ 

DATA  Y/ 52. *61. ,51. 7, 31., 12., 0., 400., 0./ 

DATA  Y/ 40. ,61. ,51. 7, 31., 12., 0., 400., 0./ 

DATA  YD/207.,202.,206.,118.,53.,0.,400.,0./ 

DATA  YD/246.,202.,206.,118.,53.,0.,400.,0./ 

DATA  S/93.0,67.,93.,72.,34.,0.,230.,0./ 

DATA  S/68. 0,67. *93. ,72., 34., 0., 230., 0  ./ 

DATA  SH/.125E-03,.125E-03».125E-03,.125E-03,.125E-03».001,1.,0./ 
DATA  NNM/7/ 

RETURN 

END 
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SUBROUTINE  LMNT CNNL , NLDCN»TH»LHPI»DT»C *PLNM ) 

INTEGER  TITL(4),MATRIAL, UNITS 

DIMENSION  T  H  ( 4  0  ) ,NAMES(4) ,PLNM(40),LMPI (40)»DESCRP(8) 

COMMON  /UNIT/  IUNIT 
COMMON  /TIT/  TITL 

NAMELlST/LAYERS/NNL,LMPI,TH,PLNM,IUNIT*DT,C,NLDCN 

NAMELIST/LAYER/NNL,TH,PLNM,DT»C,NLDCN 

DATA  NAMES  /4HPURE , 6HH Y BR I D , 2HS I , 7HENGL I SH/ 

444  FORMAT(SAIO) 

WRITE  13 

445  FORMAT! 1H0,8A10) 

WRITE  13 

READ  444  »  DESCRP 
PRINT445»DESCRP 
WRITE  13 

IF  (TITLC4)  .EQ.  NAMES ( 1 )  )  GO  TO  609 
READ  LAYERS 
GO  TO  610 

609  READ  700»MATRIAL»UNITS 

IF (UNITS  .NE.  NAMES ( 4 ) )  IUNIT=1 
READ  LAYER 
700  FORMAT ( 2 A10 ) 

CALL  MINP(MATRIAL,LI) 

DO  711  IM=1,NNL 
LMPI <IM)=LI 
711  CONTINUE 

610  CONTINUE 
WRITE  6  0  4  *NNL 

604  F0RMAT(6X»* NUMBER  OF  PLYS  =*  *  1 3> 

WRITE  605 

605  FORMAT (6X**ANGLES  OF  PLY  ORIENTATION  IN  DEGREES*  BOTTOM  L  AYER  =  1*  > 
WRITE  6051* (TH( I) *I=1*NNL) 

6051  FORMAT(6X,12F5.1) 

WRITE  665 

665  FORMAT ( 6 X,* NUMBER  OF  LAYERS  FOR  CORRESPONDING  PLY  ORIENTATION*) 
WRITE  6051, (PLNM(I), 1=1, NNL) 

WRITE  607 

607  F0RMAT(6X,*MATERIAL  PROPERTY  IDENTIFICATION  NUMBER  FOR  CORR.  PLY*) 
WRITE  6071, (LMPI(I)*I=1, NNL) 

6071  F0RMAT(6X,12I5) 

WRITE  608  ,DT,C 

608  FORMAT (6X,*TEMPERATURE  DT  =  *,F10.4,*  MOI STURE  =  * *F 1 0 . 4  ) 

13  FORMAT  (5X,* - - 

* - ★  ) 

RETURN 

END 
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SUBROUTINE  MATM ( A*B,C) 
DIMENSION  A(3,3),B(3,3),C(3*3) 
DO  2  1=1,3 
DO  2  J= 1  * 3 
SUM  =  0  •  0 
DO  1  K  =  1 , 3 

SUM=SUM+A<J,K)*B<K,I> 

1  CONTINUE 
C(J,I)=SUM 

2  CONTINUE 
RETURN 
END 


SUBROUTINE  MINP(MAT,LI) 

DIMENSION  NAMES  (8  ) 

DATA  NAMES/9HT300/5208*7HB4/5505,7HAS/3501,9HSCOTCHPLY,8HKEVLAR49, 
*4HCORE»8HALUMINUM*3HNEW/ 

DO  1  1=1,8 

1  IF (MAT  .EQ.  NAME  S  C I ) )  GO  TO  2 
STOP  "ERROR  IN  MIND" 

2  L  I  =1 
RETURN 
END 
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SUBROUTINE  MODCM ( TH • U ,Q ) 
DIMENSION  U ( 5 ) 

DIMENSION  Q  <  3 ♦ 3 ) 

TH2=2.*TH 
C2=COSM ( TH2 ) 

S2=SINMCTH2) 

TH4=4.*TH 

C4=C0SM(TH4> 

S4=SINM ( TH4 ) 

G(  ltl)=U< 1>+U(2)*C2  +  U(3)*C4 

Q(2»2>=U< 1)-U(2)*C2+U(3)*C4 

QC1*2>=U(4>-U(3)*C4 

Q(3»3>=U<5)-U(3>*C4 

Q<l»3)=U(2>*S2/2.+U<3>*S4 

Q<2,3)=U<2)*S2/2.-U(3)*S4 

CALL  SYM(Q) 

RETURN 

END 


SUBROUTINE  MODULS (EL *E T * VLT *GLT * Q*U > 

DIMENSION  U ( 5) 

DIMENSION  Q ( 3 « 3 ) 

DO  1  I=l»3 
DO  1  J  =  1  *  3 
Q  ( I  «  J )  =  0  • 

1  CONTINUE 

IF  (EL. EG).  0.)  GO  TO  2 

AM=1./(1.0-VLT*VLT*ET/EL> 

Q(1*1)=EL*AM 

Q(1»2)=AM*ET*VLT 

Q(2*1)=Q(1*2> 

Q<2*2)=AM*ET 

Q(3,3)=GLT 

2  CONTINUE 
CALL  US ( Q  *U  ) 

RETURN 

END 
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SUBROUTINE  MOLS ( THET A *HL « HN ) 

DIMENSION  EC3>*El<3>,P<3>,QGt3>,F(3>,ESll<40>tES22<40>*GS12<40) 
DIMENSION  VS12(40),SB1  (3*3)  »SB2(3,3>  «SB3<3,3)*HL<40>*HN<41)*QC3*3} 
DIMENSION  THETA(40)*HLC(40>*HT(40>*TL<40)*TT(40> 

DIMENSION  U<5>*QX<3*3> 

C0MM0N/SBS/SB1,SB2,SB3,P,QQ*DT*C 

COMMON/BLK/HLC,HT,TL,TT 

COMMON  /BLK2/  ES 1 1 *E S22 *  VS 1 2 * 6S 1 2 ,NNL 

COMMON  /BLH/  H*HI 

TH 1 =0 • 

DO  1  1=1* NNL 
TH1=TH1+HL( I  ) 

1  CONTINUE 
HTT=THl/2. 

HH=TH1*TH1 
NN=NNL+ 1 
HN (1 ) = - H  T  T 
DO  2  J=2,NN 
HN(J)=HN< J-1)+HL(J-1> 

2  CONTINUE 

DO  66  1=1*3 
P ( I  )  =0  . 

QQ ( I ) =0 • 

DO  66  J=1 *3 
SB2(I*J)=0. 

SB3(I,J)=0. 

66  CONTINUE 

DO  6  NL= 1 *NNL 
H1=HN(NL+1)-HN(NL) 

H2=H1*(HN(NL+1)+HN(NL) )/2. 

G=HN(NL+1)*HN(NL) 

H3=Hl*(Hl*Hl+3.*G)/3. 

THT  =  THETA(NL  ) 

CALL  MOGULS (ES11 <NL ) * E S22 < NL ) * VS 12 < NL > * GS 1 2 ( NL ) *QX,U> 

CALL  MODCM ( THT  *U*Q ) 

CALL  NMSN(C*DT*E*HLC*HT  *TL*TT*NL) 

CALL  MVM(QX*E*E1) 

CALL  TRS(E1,F,-THT) 

DO  4  1=1,3 
DO  3  J  =  I  *  3 

SB1(I,J)=SB1(I,J)+Q(I,J)*H1 
SB2( I* J)=SB2< I, J)+Q( I ,J) *H2 
SB3(I*J)=SB3<I,J)+Q(I*J)*H3 

3  CONTINUE 
P(I)=P(I)+FCI)*H1 
QQ(I)=QQ<I)+F<I)*H2 

4  CONTINUE 
CALL  SYM(SBl) 

CALL  SYM  <  SB2 ) 

CALL  SYM(SB3) 

6  CONTINUE 

H  =  TH1 
HI=1 ,/H 
RETURN 
END 
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SUBROUTINE  MTAD ( A *B « C, CN) 

DIMENSION  A(3»3)»B(3,3),C(3,3) 

DO  1  1  =  1*3 
DO  1  J  =  l,3 

C(I»J)=A(I,J)+CN*B(I*J) 

1  CONTINUE 

RETURN 
END 

SUBROUTINE  MTDM ( PLNM , PM ,GM , PHI , TP  I ) 

DIMENSION  PLNM(40)*HL(40) 

DIMENSION  ES11(40)*ES22(40)»VS12(40),GS12(40) 

DIMENSION  SXC(40) , SXL ( 4 0 ) , S YL ( 4 0 ) ,S Y C ( 4 0 ) , SXLT ( 4 0 ) 

DIMENSION  HLC (40) ,HT ( 4 0 ) , TL ( 4 0 ) , TT ( 4 0 ) , SH ( 8 ) 

DIMENSION  EX(8),EY(8),VX(8) ,ES(8)  * ALFX ( 8 ) * ALF Y ( 8 ) ,BTA X ( 8 ) 

DIMENSION  BTAY(8)fX(8) ,  XD(8),Y(8),  YD (8 ) ,S ( 8) ,LMPI  (4 0  ) 

COMMON  /UNIT/  IUNIT 

COMMON  /BLK 1 /  S XL , S XC , S YL , S Y C * SXLT, HL 

COMMON/ BLK/HLC,HT,TL»TT 

COMMON  / B L K 2 /  ESll,ES22»VS12*GS12»NNL 

COMMON  /DTA/  EX , E Y , V X , ES , AL FX , ALF Y, BT AX ,BTA Y , X , XD, Y , YD, S , SH , NNM 

COMMON  /MPIL/  LMPI 

CON1=1.0 

CON2=l .0 

PM=1.0E+06 

IF( IUNIT  .EQ.  1)  GO  TO  4455 
PM=PM/6895. 

C0Nl=5./9. 

C0N2=39.4 
4455  CONTINUE 

GM=1000.*PM 
DO  1  1=1, NNL 
1 1 =LMPI ( I  ) 

ES11(I)=EX(II)*GM 
ES22(I)=EY<II)*GM 
V  S 1 2  < I ) =VX (II) 

GS12(I)=ES(II)*GM 

TL(I)=ALFX(II)*CON1*1.0E-06 

TT(I)=ALFY(II)*CON1*1.0E-06 

HLC(I)=BTAX(II) 

HT( I )=BTAY( II ) 

SXL ( I ) =X ( I I ) *PM 
SXC(I)=XD(II)*PM 
S YL ( I ) =Y ( 1 1 ) *PM 
SYC(I)=YD(II )  *PM 
SXLT(I)=S(II)*PM 
HL(I)=SH(II)*C0N2*PLNM(I) 

1  CONTINUE 

IF (IUNIT •  NE.  1)  PM=  1 • 0E+03 

GM=1 000. *PM 

PMI=1./GM 

TPI=1000.*GM 

IF(IUNIT.EO.l)  WRITE  24 

IF(IUNIT.NE.l)  WRITE  25 

24  FORMAT(20X,*SI  UNITS*) 

25  FORMAT C20X**ENGLISH  UNITS*) 

RETURN 

END 
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SUBROUTINE  MVM(X*Y,Z> 
DIMENSION  X<3*3)*Y(3)*Z(3> 
DO  1  1=1*3 
SUM  =  0  •  0 
DO  2  J= 1  *  3 
SUM=SUM+X(It J)*YCd) 

2  CONTINUE 

Z(I)=SUM 
1  CONTINUE 

RETURN 
END 


SUBROUTINE  M X STRN < K * EN » A AN , Q * PM ) 

DIMENSION  EPST(3) ,EPSC(3) *SI6R(3) *SIGA( 3) *EN(3) *AANC3) , AK (3) 
DIMENSION  SXL<40)*SXC<40>*SYL<40>,SYC<40>*SXLT<40>*A(3*3)»G(3*3> 
DIMENSION  HL ( 4  0  ) 

COMMON  /BLK1 /  SXL *SX C * S YL * S YC * SXLT * HL 

COMMON  / BLH/  H*HI 

COMMON  /UNIT/  IUNIT 

COMMON  /SF/ SFL*  COUNT 

CALL  INVRS(Q*A) 

EPST (1)=SXL(K)*A(1*1) 

EP$C(1>=-SXC(K)*AC1,1> 

EPST (2)=SYL(K)*A(2*2) 

EPSC(2)=~SYC<K)*A(2*2> 

EPST (3)=SXLT  <K)*A(3*3) 

EPSC<3)=-SXLT(K)*A<3,3> 

DO  1  1=1*3 
DEN=EPST(I) 

IF  «  EN ( I >  .LT.  0.)  DEN  =  EPSC(I> 

AK ( I )=EN ( I ) /DEN 

1  CONTINUE 

CALL  AMAX(AK,3,AMX) 

DO  2  11=1*3 

SIGR(II)=AAN(II)/AMX 

SIGA(II)=SIGR(II)/H/PM 

2  CONTINUE 
SFL=1./(AMX*H*PM) 

IF  C COUNT • EQ • 1  •  )  RETURN 
WRITE  6»(SIGR(I)*I=1*3) 

6  FORMAT (5X**STRESS  RESUL T ANT S* * 9X * *N 1  =  *  * E8 . 2 *  4 X  * *N2  =  * < E8 . 2 « 4 X  * 
U*N6=*,E8.2) 

IF< IUNIT  .EQ.  1)  WRITE  7* < S I G A ( I ) *  I  =  1  * 3 ) 

IF ( IUNIT  .NE.  1)  WRITE  8 * C S I G A C I ) *  I  =  1 , 3 ) 

7  FORMAT (5X»*AVERAGE  STRESSES  MPA  *  *  * S IGM A 1  =  * *F6 . 1 , 2X * *S I GMA 2  =  *  * 
»F6.1*2X**SIGMA6=*,F6.1> 

8  F0RMAT(5X,*AVERAGE  STRESSES  KSI  * *  * S I GM A 1  =  * * F6 . 1  *  2 X  * *S I GM A2  =  * * 
*F6.1*2X**SIGMA6=*,F6.1) 

RETURN 

END 
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SUBROUTINE  MXSTRS (K, AN*AAN*PM> 

DIMENSION  SIGT(3) *S  I GC  (  3  >  *  S  I GR (3) *SIGAC3)*AN<3> *AAN(3>* AKC3) 
DIMENSION  SXL(40)*SXC(40)*SYL(40)*SYC(40)*SXLT(40) 

DIMENSION  HL  <  4  0 ) 

COMMON  /BLK 1 /  SXL *SXC « S YL * S YC * SX LT* HL 

COMMON  /BLH/  H,HI 

COMMON  /UNIT/  IUNIT 

COMMON  /SF/SFL  *  COUNT 

SlGT<l)=SXL<K> 

SIGT<2)=SYL(K) 

SIGT(3)=SXLT(K) 

SIGC(l) =-SXC(K) 

SIGC<2)=-SYC(K> 

SIGC(3)=-SXLT  (K) 

DO  1  1  =  1*3 
DEN  =  S I GT ( I ) 

IFCAN(I)  .LT.  0.)  DEN=S I GC ( I ) 

AK ( I )=AN( I) /DEN 

1  CONTINUE 

CALL  AMAX  <  AK*  3*  AMX) 

DO  2  11=1*3 

SIGR(II)=AAN(II)/AMX 

SIGA(II)=SIGR(II)/H/PM 

2  CONTINUE 
SFL=1./(AMX*H*PM) 

I F  <  COUNT • EQ • 1  •  )  RETURN 
URITE  6*(SIGR(I),I=1,3) 

6  FORMAT (5X**STRESS  RESULTANTS* *9X * *N 1  =  *  * E8. 2  *  4 X * *N2  =  * * E8 . 2 *4 X, 
U*N6=**E8.2> 

IF ( IUNIT  .EQ.  1)  WRITE  7, ( S I G A C I ) *  I  =  1  *  3 ) 

I F ( I UN I T  .NE.  1)  WRITE  8 , < S I G A C I ) *  I  =  1  * 3 ) 

7  F0RMAT(5X**AVERAGE  STRESSES  MPA  * *  * S IGMA 1  =  * ,F 6 . 1  * 2X * *S IGMA2  =  * * 
,F6.1*2X**SIGMA6=**F6.1  ) 

8  F0RMAT(5X»*AVERAGE  STRESSES  KSI  * *  * S IGM A 1  =  * * F 6 • 1  * 2X * *S I GM A2= *  * 
*F6.1*2X,*SIGMA6=*,F6.1) 

RETURN 

END 


SUBROUTINE  NMSN (C ,DT ,EN1 , HL , HT *TL*T T, K ) 
DIMENSION  EN1 (3) *HL(40) ,HT(4  0) *TLC40>  *TTC40) 
EN1 (1)=C*HL(K)+DT*TLCK) 
EN1<2)=C*HT(K)+DT*TT(K> 

EN1 (3)=0. 

RETURN 

END 
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SUBROUTINE  NORM < PMI * ZER 0* DU 1 *DU2 * DU3 * DU4 ,P AR > 

INTEGER  PAR 

DIMENSION  NAMES(6)*A(3*3) 

DIMENSION  DU1<3*3)*DU2<3,3)*DU3(3*3)*DU4<3*3),ZER0<3*3) 

DIMENSION  SB1 <3*3)*SB2<3*3>  *SB3(3*3)*P(3)  *  QQ ( 3 ) *HLC(40)  *HT ( 4  0 ) 

DIMENSION  TL(40) ,TT<40) *DUM(3*3) ♦ DUM 1 ( 3 * 3 ) * E CON C 4 ) * E CONF C 4 ) 

COMMON  /UNIT/  IUNIT 

COMMON  /BLH/  H*HI 

C0MM0N/SBS/SB1*$B2*SB3*P*QQ*DT*C 

COMMON/BLK/HLC*HT*TL*TT 

DATA  NAMES  / 1 OHDI MENS  I ONL * 1 OHNOR M AL I Z ED  * 4HB0TH * 8HENGC0NST * 3H ALL * 
*4H  / 

13  F0RMATC5X** - 

—  -  —  —  —  — 

HSQ=H*H 

HIS=1./HSQ 

PMN 1=PM I / H 

PMN2=2.*PMI/HSQ 

PMN3=12.*PMN1/HSQ 

CALL  INVRS(SBlfA) 

CALL  MTAD(ZERO,A*DUM,H) 

EC0N(1)=PMI/DUM<1,1) 

EC0N<2)=PMI/DUM(2*2) 

EC0N(3)=-DUM(1*2) /DUM (1,1) 

EC0N(4)=PMI/DUM<3,3) 

CALL  INVRS(SB3*DUM1 ) 

HQTW=HSQ*H/12. 

CALL  MTAD (ZERO, DUM1, DUM. HQTU) 

EC0NF(1)=PMI/DUM<  1*1) 

ECONF (2)=PMI/DUM(2*2) 

EC0NF<3)=-DUM(1,2)/DUM(  1*1) 

EC0NF(4 )=PMI/DUMC3*3) 

CALL  MTAD(ZER0,SB1*DU1*PMN1  ) 

CALL  MTAD(ZER0,SB2*DU2*PMN2) 

CALL  MTAD(ZER0,DU2*DU3»3. ) 

CALL  MTAD(ZER0,SB3*DU4,PMN3) 

WRITE  13 


I F  (  P  AR 

.EQ. 

N  AMES  <  6 ) ) 

GO 

TO 

725 

I F ( P AR 

.EQ. 

NAMES<4> ) 

GO 

TO 

724 

I F ( P AR 

.EQ. 

NAMES (2) ) 

GO 

TO 

723 

WRITE  14 
WRITE  15 

14  FORMAT ( 35  X*  *  A  B**//> 

15  FORMAT  <  35  X  *  *  B  D*) 

WRITE  13 

CALL  WITE (SB1,SB2) 

CALL  WITE (SB2,SB3) 

WRITE  13 

IF (PAR  .EQ.  N AMES ( 1 )  )  GO  TO  725 
723  CONTINUE 
WRITE  42 

42  FORMAT  C  35X, *  A#  B#*,//) 

IF ( IUNIT  .EQ,  1)  WRITE  43 
IF( IUNIT  .NE.  1)  WRITE  443 


4G 


GPA*  ) 

l.E+06  PSI*) 


TO  725 


NG  CONSTANTS  *  ES  IN  GPA*) 

IFCIUNIT.NE.l)  WRITE  4461 

4461  FORMAT (6X**SOME  ENGINEERING  CONSTANTS*  ES  IN  l.E+06  PSI*) 

WRITE  447*(EC0N<IE) *IE=1,4) 

447  F0RMAT(6X**INPLANE  E1=*,F7.3**  E2=*,F7.3,*  V21=**F7.3**  E6=* 

*,F7.3) 

WRITE  448,<EC0NFCIE),IE=1,4) 

448  F0RMAT(6X**FLEXURAL J *EF 1  =  * * F 7 . 3 , *  EF2=*,F7.3,*  VF2 1  =  *  * F 7. 3 *  *  EF6  =  * 
*  *F7. 3 ) 

725  CONTINUE 
RETURN 
END 


43  FORMAT (35X**3B#  D# 

443  FORMAT  C35X**3BU  DU 

WRITE  13 

CALL  WITE1(DU1*DU2) 

CALL  WITE1(DU3*DU4) 

WRITE  13 

IF(PAR  .NE.  NAMESC5))  60 
724  CONTINUE 

IFCIUNIT.EQ.l)  WRITE  446 
446  FORMAT(6X**SOME  ENGINEERI 
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SUBROUTINE  NORM  1 C TP I  * A *P AR AM ) 

INTEGER  PARAM 

DIMENSION  DU1(3*3)*DU2(3*3)*DU3(3*3) 
DIMENSION  SB1 (3*3) *SB2(3*3) *SB3(3*3) 
DIMENSION  TL(40) *TT (40) *DU5(3*3) *DU6 
DIMENSION  DU9(3*3)*DU10(3*3) *PNVC3) * 
DIMENSION  NAMES<4)*D7( 3,3) *DM<3*3) 
COMMON  /UNIT/  IUNIT 
COMMON  /BLH/  H,HI 
COMMON/SB S/SB1,SB2,SB3*P,QQ,DT,C 
COMMON/BLK/HLC,HT,TL*TT 
COMMON  /BLDU/  DU2 ,DU4 , D U5 * D U7 
DATA  NAMES  / 1 OHD I  MENS  I ONL *  1 OHNORM AL I 
FORMAT  (5X,* - - 


,DU4<3,3) ,ZER0C3,3) 

♦  PC  3) ,QQ(3),HLCC40),HT(40) 

C3,3)*DU7C3,3),DU8C3,3) 

0NV(3),DUM(3*3) 


ZED,4HB0TH,4H  / 


HSG=H*H 


HI=1./H 

HIS=1./HSQ 

TPN1=TPI*H 

TPN2=HSQ*TPI/2. 

TPN3=TPN1*HSQ/12. 

TPNN=TPN2/3. 

DO  75  1=1,3 
DO  75  J  =  1 ,3 
ZEROCI*J)=0. 


47 


CALL  MATM(A,SB2*DU1> 

CALL  MATM(SB2,DU1,DU2) 

CALL  MTAD(SB3,DU2,DU3,-1.) 

CALL  INVRS(DU3,DU4> 

CALL  MATM<DU1,DU4,DU5) 

CALL  MATM(SB2»A,DU1) 

CALL  MATM(DU5*DU1,DU6) 

CALL  MATM(DU4,DU1,DU2> 

CALL  MTAD(A,DU6,DU7,1. > 

CALL  MTAD < ZERO ,DU5 ,DUM , -1 . > 

CALL  MTAD(ZER0,DU7,D7, 1.) 

CALL  MTAD(ZERO,OUM,DM,l.) 

CALL  MTAD(ZER0,DU7,DU1*TPN1> 

CALL  MTAD(ZER0,DUM*DU8,TPNN) 

CALL  MTAD(ZER0»DU2,DUM»~1«) 

CALL  MTAD <ZER0,DUM,DU9,TPN2 ) 

CALL  MTAD(ZERO,DU4,DU10,TPN3) 

I P ( PAR  AM  «EQ • NAMES ( A ) )  GO  TO  74 
IF (PARAM  .EQ.NAMES(l))  GO  TO  71 
IF (PARAM  .EQ.NAMES<2) >  GO  TO  72 

71  CONTINUE 
WRITE  13 
WRITE  16 
WRITE  17 

16  FORMAT(35X,*ALPHA  BETA  *,//> 

17  F0RMAT(35X»*TRSBTA  DELTA*) 

WRITE  13 

CALL  WITE(D7,DM) 

CALL  WITE <DUM,DU4> 

WRITE  555 

555  F0RMAT(6X,*N0NMECH.  STRESS  AND  MOMENT  RESULTANTS  N*M*> 
CALL  WRT (P,5) 

CALL  WRT(QQ,5) 

WRITE  13 

IF(PARAM  ,NE.NAMES(3) )  GO  TO  74 

72  CONTINUE 
WRITE  44 

44  FORMAT (35X,*  ALPHA#  BETA#/3*,//) 

I F  <  I  UNI T  .EQ.  1)  WRITE  45 

I F (  I UN I T  .NE.  1)  WRITE  4451 

45  FORMATC35X,*  TRSBTA#  DELTA#  (l./TPA)  *) 

4451  FORMAT(35X,*  TRSBTA#  DELTA#  1 . OE - 09/PS I  * ) 

WRITE  13 

CALL  WITE1(DU1,DU8) 

CALL  WITE1(DU9,DU10) 

DO  557  I ST= 1 »  3 
PNV (IST)=P( I  ST ) *H I 
QNV< IST)=6.*QQ< IST)*HIS 
557  CONTINUE 

WRITE  556 

556  FORMATC 6X,*N0NMECH.  EFFECTIVE  STRESS  AND  MOMENT  N#,M#*> 
CALL  WRT (PNV, 5) 

CALL  WRT  <QNV*5> 

WRITE  13 
74  CONTINUE 

RETURN 
END 


48 


SUBROUTINE  PLT(XyXlyXXyX2yXXXyX3yXCyXClyXXCyXC2yXXXCyXC3yNyTHyIC) 
INTEGER  SMBL  y IMMy PL ANE 

DIMENSION  X(200)yXl(200)yXX(200)yX2(200)yXXX(200)yX3(200)yXC(200)y 
*XC1  (200) yXXC(200) yXC2(200) y XXXC(  2  00) yXC3(200) 

DIMENSION  TH(40)yNAME(2) 

COMMON  / BKK 9/  PX X »PX X I • P Y Y y P Y Y I y YUL « I T 
COMMON  /SMB/  SMBLyIMM 

COMMON  /CRIT/  CR I TR I  A , SP ACE y PL ANE y F S 1 2 
DATA  NAME  /9HSUPERPOSE y 2HQR / 

N2=N+1 
N3=N+2 
X(N2)=PXX 
X(N3)=PXXI 
XI ( N2 ) =  PYY 
X1(N3)=PYYI 
P YS  =PYY/P  Y Y I 
PXS=PXX/PXXI 
DD  Y  =  2  »  55 

IF (PLANE  .E  Q •  NAME(2)  )  DDY=3.55 

DDYl=DDY-.3 

DDY2  =  DDY-  .6 

IF (SMBL  .EQ.  NAME ( 1 )  .AND.  IMM  .GT.  1)  GO  TO  3 

CALL  SPAXIS(0.yPYSylH  y 1 y YUL y 9 0 . y XI ( N2 > y XI ( N3 ) y 0 . 2 y 3 . y . 2 y 0 . y 2 . y I T y 
*.15) 

CALL  SPAXIS(PXSyO*ylH  y-ly7.0y0.yX(N2)yX(N3)y2.25y.ly.25y0.y2.yITy 
*.15) 

3  CONTINUE 

IF ( IC  .EQ.  1)  GO  TO  1 
I F ( I C  .EQ.  2)  GO  TO  2 
XC(N2)=0. 

XC(N3)=PXXI 
XC1 (N2)=0. 

XC1 (N3)=PYYI 
XXC(N2)=0. 

XXC(N3)=PXXI 

XC2(N2)=0. 

XC2(N3)=PYYI 

XXXC(N2)=0. 

XXXC(N3)=PXXI 

XC3(N2)=0. 

XC3 ( N3 ) =PYY I 

IF  (SMBL  .EQ.  NAME(D)  GO  TO  4 
CALL  NUMBER (1.0 y DDY y . 15 y TH ( 1 )  y 0 . y 1 ) 

CALL  FLINE(XCyXCly-Nylyl0y5) 

CALL  NUMBER ( 1.0 yDDYl y . 15y TH( 2) yO.yl) 

CALL  FLINE(XXCyXC2y~Nylyl0y6) 

CALL  NUMB ER ( 1.0 y DDY2 y . 1 5 y TH ( 3  )  y 0  .  y 1  ) 

CALL  FLINE(XXXCyXC3y-Nylyl0y7) 

GO  TO  1 

4  CONTINUE 

CALL  FLINE(XCyXCly-Nyly0y5) 

CALL  FLINE(XXCyXC2y-Nyly0y6) 

CALL  FLINE(XXXCyXC3y-Nyly 0y7) 

1  CONTINUE 
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X3<N2)=0. 

X3<N3)=PYYI 

XXX<N2)=0. 

XXX (N3) =  PXXI 
X2(N2)=0. 

X2<N3)=PYYI 
XX (N2)=0. 

XX(N3)=PXXI 

I F ( SMBL  .EQ.  NAME ( 1 ) )  60  TO  5 
CALL  SYMBOL  < .5,  DDY2 , . 15 , 7 « 0 . , - 1  ) 
CALL  FLINECXXX,X3,-N,1,10,7) 

CALL  SYMBOL  < ,5*  DDY 1 , . 15 , 6 , 0  .  , - 1 ) 
CALL  FLINE(XX,X2,-N,1,10,6) 

GO  TO  2 

5  CONTINUE 

CALL  FLINE<XXX,X3,-N, 1,0,7) 

CALL  FLINE(XX,X2,-N, 1,0,6) 

2  CONTINUE 

X1(N2)=0. 

X<N2)=0. 

I F ( SMBL  .EQ.  NAME(l))  GO  TO  6 
CALL  SY MBOL ( • 5 ,  DDY  , . 15 , 5 , 0 . , -1 ) 
CALL  FLINECX,X1,~N,1,10,5) 

CALL  PLOT (13.,0.*3) 

CALL  PLOT (13., 0., -2) 

RETURN 

6  CONTINUE 

CALL  FLINE (X,X1,-N,1, 0,5) 

RETURN 

END 
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SUBROUTINE  PLTEN ( TH , HL * HS * PM I  * ZE RO, PLNM > 

INTEGER  P AR AM2*PARAM  3 
DIMENSION  NAMEC3) ,PLNM(20) 

DIMENSION  SB1<3,3),SB2<3,3) , SB3 < 3 ,3 ) , P ( 3 ) ,QQ(3) »  DUM  (  3 , 3  )  ,  X  (  2  0  0  ) 
DIMENSION  ZERO (3, 3) * ECON1 < 2 0 0 ) , E CON2 ( 2 0 0 ) , ECON3 ( 20 0 ) * ECON A C 2 0 0 ) 
DIMENSION  TH(4O),HL(40),HLC(40)*HT(4O),TL(4O)»TT<4O)»HS<41) 
DIMENSION  ES11(40),ES22(40),GS12(40)*VS12(40),A<3,3) 

COMMON  /BLK 2/  ES 1 1 , ES22 , VS 1 2 ♦ GS 1 2 ,NNL 

COMMON  /BLH/  H,HI 

COMMON  /UNIT/  IUNIT 

COMMON  /SBS/SB1,SB2,SB3,P,QQ,DT,C 

COMMON  /BLK/HLC,HT,TL,TT 

COMMON  /BKK9/  P XX ,PXX I , P Y Y * P Y Y I , YUL « I T 

COMMON  /BLM/  P AR AM2 ,P AR AM 3 

DATA  NAME/6HPLTEND,8HPLTSTART»6HPLT0NE/ 

IF (PARAM2  .EQ.  NAME ( 3 ) )  GO  TO  7 
IF(PARAM2  . NE •  NAME  C  2) )  GO  TO  6 
7  CONTINUE 

CALL  PL0TS<DUM,DUM,99) 

CALL  PLOT (5«,2.*3) 

CALL  PL0T<5.,2.,-2) 

6  CONTINUE 

DO  1  11=1*91 
X< II)=TH(3) 

CALL  MOLS (TH»HL»HS) 

CALL  INVRS(SB1,A) 

CALL  MTAD(ZERO*A*DUM*H) 

ECONKII  )=PMI/DUM(1,1) 

EC0N2(II)=PMI/DUM<2,2) 

EC0N4<II)=PMI/DUM(3,3) 

EC0N3(II)=-DUM(1,2)/DUM<1,1) 

TH(3)=TH<3)+1. 

TH(4)=-TH<3) 

TH(5)=TH<4) 

TH(6)=TH(3) 

1  CONTINUE 

YUL=6.0 
PXX=0. 

PXXI=15.0 

PYY=0. 

CONT  =  1 • 

I F ( I  UN  I T  .EQ,  1)  CONT  =  10.0 
PYYI=5.*CONT 

IF (IUNIT  .EQ.  1)  CALL  SYMB OL ( . 1 , 4 .5 , . 15 , 3HGP A , 0 . , 3 ) 

IFdUNIT  .NE.  1)  CALL  SYMBOL C . 1, 4 .5  * . 15 ,3HMSI *  0. *3) 

CALL  SYMBL<3. ,5. * TH * PLNM , NNL * 0) 

CALL  SYMBOL(0.2,5.0,.2,1HE,0.,1) 

CALL  SYMBOLS. 2, 0.2, .2*140, 0.*-l> 

CALL  SYMBOL (1  .0 0 , 2 .55 , .15 , 1  HE  ,  0 . , 1 ) 

CALL  SYMBOL (1.15, 2.4 5,.15,2H 11,0. ,2) 

CALL  SYMBOL(1.00,2.25,.15,1HE,0.,1> 

CALL  SYMB0L(1.15*2.15,.15,2H22,0.,2) 

CALL  SYMBOL < 1.0  0,1.95, .15,1 HE,0.  ,1) 

CALL  SYMB0L(1.15,1.85,.15,2H6  ,0.*2> 

CALL  PLT(X,ECON1,X,ECON2,X,ECON4,X,ECON1,X,ECON2,X,ECON4,91,TH,1) 
PYY 1=0.2 

CALL  SYHBL(3.,5.*TH,PLNM,NNL,0) 

CALL  SYMBOL (6. 2, 0.2, .2,140,0.,-!) 

CALL  SYMBOL (0.2, 5.0,. 25,1 33, 0.,-l) 
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CALL  SYMBOL  <0.35*4.9,.10.2H21*0.,2> 

CALL  S YMBOLC 1.00* 2. 55*. 15,133,0. *-l  ) 

CALL  SYMBOL<1.15*2.45*.15,2H21,0.,2) 

I  T  =  1 

CALL  PLT<  X*EC0N3*X,EC0N2*X,EC0N4,X*EC0N1,X,EC0N2*X,EC0N4,91,TH,2> 
I  T=  0 

IF(PARAM2.EQ.NAME(3) )  GO  TO  8 
IF ( PAR AM2  .NE.  NAME  C 1 ) )  GO  TO  5 
8  CONTINUE 

CALL  PLOTE(N) 

5  CONTINUE 

WRITE  3 

3  FORMAT(10X**PLOTS  FOR  ENGINEERING  CONSTANTS  DRAWN*) 

RETURN 

END 


SUBROUTINE  R OOT S ( A , B * C , SM , SMN ) 
A2=2. 0*A 
BS-B  *B 

DC=BS-2.0*A2*C 
DS^SQRT ( DC ) 

SM=<DS-B)/A2 

SMN=-(R+DS)/A2 

RETURN 

END 
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SUBROUTINE  SE TAN < AN t DCN ) 

INTEGER  PLANE»SPACE,CRITRIA 
DIMENSION  AN(3) ,NAME(2> 

COMMON  /CRIT/  CRI TR I  A  *  SPACE *PL ANE *FS1 2 
DATA  NAME/2HQR*6HSTRAIN/ 

IF ( PLANE  .EQ.  NAME ( 1 ) )  GO  TO  1 
AN(l)=COSM(DCN) 

AN(2)=SINMCDCN) 

AN ( 3 )  =  0  » 

RETURN 

1  CONTINUE 

SQ2=SQRTC2. ) 

AN(1)=SQ2*C0SM<DCN) 

AN(2)=-ANC1) 

AN(3)=SQ2*SINM(DCN> 

IF (SPACE  .EQ.  N AME ( 2  )  )  AN < 3 ) =AN ( 3 > /2 . 

RETURN 

END 


SUBROUTINE  STPLT ( TH * HS • PM ♦PLNM ) 

INTEGER  PARAM2,PARAM3*CRITRIA,SPACE*PLANE*SMBL*MULTI 
DIMENSION  VN ( 3 ) *  PLNM (40) 

DIMENSION  GS(3,3)*GSB(3) 

DIMENSION  NAME ( 14  ) 

DIMENSION  Q(3*3)*HL(40)*TH(4Q)*SFLN(40> 

DIMENSION  SPLC(40),STX1  (200) *  STY 1 C 2 0 0 ) , SCX 1 C 20 0 ) *SCY1(200) 
DIMENSION  STX2(200) *STY2(200) »SCX2(200) ,SCY2(200) 

DIMENSION  STX3(200> ♦STY3<200  )*SCX3<200>  *SCY3(200> 

DIMENSION  GB (3) »G<3*3> ,SFLL<40> *SFU<40> *EN1C3) ♦ENDC3)* QQ(3) 
DIMENSION  SB3(3,3) » HT < 4 0 > » TL < 4 0 > *TT<40> ,PC3) *HLC(40> 

DIMENSION  SXC(40) ♦SXL<40> ,SYL<40) *S YC ( 4 0 ) , SXLT C 4 0 ) 

DIMENSION  ESll(40)tES22(40),VS12(40),GS12(40),SBl(3»3)»SB2(3*3) 
DIMENSION  DU2C3*3)*0U4(3*3> » DU5 < 3 t 3 > * E 1 < 3 > » E2 < 3 ) *E3C3)*E4(3) 
DIMENSION  ETR<3> ,ENT(3) *U(5) *DU7<3*3) tHS<21)  ,  AN < 3), AM  (3) 

COMMON  /BLK 2/  ESI 1 ♦ ES2 2 f VS1 2  * GS 1 2 »NNL 

COMMON /BLK /HLC»HT,TL,TT 

COMrtON/SBS/SBl* SB2tSB3»P*QQ*DT*C 

COMMON  /BLK 1 /  SX L » SXC * S YL ♦ S YC * SXLT * HL 

COMMON  /UNIT/  IUNIT 

COMMON  /BLH/  H,HI 

COMMON  /BLDU/  DU2 »DU4 »DU5 *DU7 

COMMON  /BKK9/  PXX * PX X  I  * P Y Y » P Y Y I , YUL *  I T 

COMMON  /SMB/  SMBL  *  I  MM 

COMMON  /BLM /  P AR AM2 « P AR AM  3 

COMMON  /SF/SFL*  COUNT 

COMMON  /CRIT/  CRITRI A , SP ACE ♦ PL  AN E » F S 1 2 
COMMON  /ST1 /  STX1*STY1*SCX1*SCY1*NSF1 
COMMON  /ST2/  STX2 ,ST Y 2 *SC X2 ♦ SC Y2 *NS F2 
COMMON  /ST3/  ST X3 *ST Y3 * SCX3 *  SC Y3 » NS F3 
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DATA  NAME/6HPLTEND,8HPLTSTART,7HTSAI  UU ,9HM A XSTR ESS » 
*9HMAXSTRAIN*6HSTRAIN» 6HSTRESS»6HCHAMIS*7HH0FFMAN*4HHILL*6HPLT0NE, 
*2HQR,9H SUPER  POSE ♦9HMULTICURV/ 

COUNT=l. 

FCTRS=1 • 

FC  T  RF  =  1 • 

READ  705,CR I JRIA* SPACE, PLANE *  SMB L *F S 1 2 ♦ MULT  I »FC TRS *FCTR F 
7  05  FORMAT  (4A10»F10.3 «A10  *2F10.3  ) 

IF(PARAM3  .EQ.  NAME(ll))  GO  TO  716 
IF (PARAM3  .NE.  N AME ( 2  )  >  GO  TO  701 
716  CONTINUE 

CALL  PL0TS(DUM»DUM,99> 

CALL  PL0TC5. ,5. ,3) 

CALL  PLOT (5 • «5 • *-2  ) 

I F ( FCTRF  .EQ.  0.)  FCTRF=1. 

CALL  FACTOR (FCTRF) 

701  CONTINUE 
YUL=8 • 

I  T=  0 

IF (FCTRS  .EQ.  0.)  FCTRS=1. 

CONT=FCTRS*l  . 

IFCIUNIT  .EQ.  1)  CONT=FCTRS*10.0 

PXX=-400.0*CONT 

PXXI=100. 0*CONT 

PYY=-500.0*CONT 

PYYI=PXXI 

I F ( PL  AN  E  .NE.  NAME(12) )  GO  TO  726 
PXX=-75.0*CONT 
C  PXX=-150. 0*CONT 

C  PXXI=50.0+CONT 

PXXI=25.0*CONT 
C  PYY=-200. 0*CONT 

PYY=-100.0*CONT 
PYYI=PXXI 
726  CONTINUE 

NPT=91 
MM=NNL/2 

I F ( SMB  L  .NE.  NAME(13) )  MM  =  1 
DO  722  I MM= 1  *  MM  »  3 
725  CONTINUE 

DO  601  I CND=1 »NPT 
DCN=18  0* (ICND-1 )/(NPT-l ) 

CALL  SETAN<  AN*OCN) 

AM ( 1 ) =0 • 

AM ( 2 ) =0 • 

AM ( 3 ) =0 • 

NSF1 =1 
NSF  2  =  2 
NSF3=3 

DO  600  I K  =N SF 1 » NSF3 

K= IMM+ I K-  1 

I F ( K  .GT.  NNL)  K=1 

IF(ES11(K)  .EQ.  0.0)  GO  TO  600 

DO  109  I  =  1 »  3 

SUM  1  =  0 • 

SUM2=0. 

SUM3=0. 

SUM4=0. 

DO  209  J  =  1 »  3 

SUM1=SUM1+DU7< I, J)*AN< J)-DU5(It J)*AM( J > +HS ( K ) * ( -DU2 ( I t J ) * AN ( J > + 
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*DU4CIt J)*AMCJ>) 

SUM2=SUM2+DU7CItJ)*ANCJ)-DU5CItJ)*AMCJ)+HSCK+l)*C-DU2CItJ)*ANCJ)+ 
*DU4CItJ)*AMC J) ) 

SUM3=SUM3+DU7C It J>*  P C J ) -DU5 ( I t J ) *QQ C J) +  HS ( K ) * ( -DU2 ( I t J ) *  PCJ)  + 
*DU4CltJ)*QQCJ)) 

SUM4=SUM4+DU7C It J)*  PCJ)-DU5CItU)*QQCJ)+HSCK+l)*C-DU2CItJ)*  PC  J)  + 
*DU4 Cl ,J)*QQCJ)  ) 

209  CONTINUE 

El C I ) =SUM1 
E2CI)=SUM2 
E3CI)=SUM3 
E4C I)=SUM4 
109  CONTINUE 

C  E 1 ( J )  =MECHANI CAL  STRAIN  COMPONENTS  LOWER  SURFACE 

C  E 2 ( J )  =MECHANICAL  STRAIN  COMPONENTS  UPPER  SURFACE 

C  E4  C  J )  =NONMECHANICAL  STRAIN  COMPONENTS  UPPER  SURFACE 

C  E3( J)  =NONMECHANICAL  STRAIN  COMPONENTS  LOWER  SURFACE 

THT=TH(K) 

CALL  TRE(EltETR,THT> 

CALL  TRE  <  E3  tENTt  THT ) 

CALL  NMSN(CtDT,ENltHLCtHTtTLtTTtK) 

CALL  VDI (ENTtENltEND) 

CALL  M0DULS(ES11 (K) tES22(K ) tVS12(K>  tGS12CK) tQtU) 

IF CCRITRIA  .EQ.  NAME  C  4 )  .AND.  SPACE  .EQ.  NAME  C  7) >  60  TO  702 

IFCCRITRIA  .EQ.  NAMEC5)  .AND.  SPACE  .EQ.  N AMEC  7) )  GO  TO  703 

I F ( SPACE  .EQ.  NAMEC6)  )  GO  TO  708 
IF ( SPACE  .EQ.  NAMEC7) >  GO  TO  709 

708  CALL  TRE (ANtETRtTHT) 

IFCFCTRS  EQ.  0.0)  FCTRS=1. 

PXX=-.04*FCTRS 

PXXI=.01*FCTRS 

PYY=~.05*FCTRS 

PYYI=PXXI 

I  T  =  2 

IFCCRITRIA  .EQ.  N AME ( 4 ) )  GO  TO  702 
IFCCRITRIA  .EQ.  NAMEC5) )  GO  TO  703 

CALL  FAILCOCETRtQtSXLCK) tSXCCK) tSYL(K) tSYCCK) tSXLTCK) tGtGB) 

CALL  COEFCGtGBtGStGSBtTHT) 

CALL  FSFTYC ANtENDtGStGSBtSMtSMN) 

SFLN  C IK ) =SM 
CALL  VNFCAN,VN) 

CALL  TRE  CVNtETRtTHT) 

CALL  FAILCOCETRtQtSXLCK ) tSXCCK) tSYL(K) tSYC(K) tSXLTCK ) tGtGB) 

CALL  COEFCGtGBtGStGSBtTHT) 

CALL  FSFTYC ANtENDtGStGSEtSMtSMN) 

SPLC  C IK )  =  SMN 
GO  TO  600 

709  CONTINUE 

CALL  FAILCOCETRtQtSXLCK) tSXCCK) tSYLCK) tSYCCK) tSXLTCK) tGtGB) 

CALL  FSFTYCETRtENDtGtGBtSMt SMN) 

SFLNCIK)=SM/CH*PM) 

IFCCRITRIA  .NE.  NAME  C  8 )  )  GO  TO  710 
CALL  VNFCETRtVN) 

CALL  FAILCO  CVNtQtSXL  CK) tSXC  CK) tSYLCK ) t S YC C K ) t SX LT C K ) t G t GB  ) 

CALL  FSFTYCETRtENDtGtGBtSMtSMN) 

710  CONTINUE 
SFLLCIK)=SM 
SPLCCIK)=SMN/CH*PM) 

GO  TO  600 

702  CONTINUE 
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CALL  MVM ( Q , ETR  *  El ) 

CALL  MXSTRS(K,E1,AN,PM> 

IF ( SPACE  .EG.  NAME ( 6  )  )  SFL=SFL*H*PM 
SFLN( IK ) =SFL 
CALL  VNF  <  El  *  VN ) 

CALL  MXSTRS(K»VN,AN»PM) 

IF  <  SPACE  .EQ.  N AME( 6  > )  SFL=SFL*H*PM 
SPLCC IK)=-SFL 
GO  TO  600 
703  CONTINUE 

CALL  MXSTRN(K,ETR,AN*Q,PM> 

IF ( SPACE  .EQ.  NAME(6>)  SFL=SFL*H*PM 
SFLN ( IK ) =SFL 
CALL  VNF  (ETR  »VN ) 

CALL  MXSTRN ( K , VN , AN , Q , P M ) 

I F ( SP AC  E  .EQ.  N AME  <  6 ) )  SFL=SFL*H*PM 
SPLC(IK) =  -SFL 

600  CONTINUE 

CALL  ADJUST (ICND»SFLN,SPLC,AN) 

601  CONTINUE 

I F ( I  MM  .  GT .  1)  GO  TO  721 

IF ( PLANE  .EQ.  NAME<12)>  GO  TO  718 

IF ( SMBL  .NE.  NAME (13 ) )  CALL  S YMBL ( - 3 . 0 , 2 . 75 , TH ,PLNM , NNL , 1 > 

C  CALL  RECT(-4.,-5.,8.0,6.0,0.,3> 

IF (SPACE  .EQ.  N AM E ( 7 ) )  GO  TO  711 
CALL  SYMBOL(-.3*2.4,.25,129,0.,-1) 

CALL  SYMBOL (1.6 *.35,. 25* 129 ♦0.*-l> 

GO  TO  712 

711  CONTINUE 

CALL  SYMBOL (-.3, 2. 4*. 25, 108,0. *-l) 

CALL  SYMBOL(1.6,.35,.25,108,0.*-1> 

IF  < IUNIT  .EQ.  1)  CALL  S YMB OL (- . 5 » 2 . , . 15 , 3HMP A , 0 . , 3 > 

IFdUNIT  .NE.  1)  CALL  S  YMBOL  <  -  .5 , 2.  ,  .  15 , 3HK  S  I  ,  0  .  ,  3 ) 

712  CONTINUE 

CALL  SYMBOL (-.1,2.30, .10,1 H 2,0.* 1) 

CALL  SYMBOL(-.1,2.55».10*1HO*0.,1> 

CALL  SYMBOL (1.8*0.25, .10, 1H1  ,0.,1  ) 

CALL  SYMBOL(1.8,0.50,.10,1HO,0.,1> 

GO  TO  719 

718  CONTINUE 

I F ( SMBL  .NE.  NAME (13 ) )  CALL  S YMB L ( - 2 . 3 , 3 . 75 , T H , PLNM , NNL , 1 ) 

C  CALL  RECT (-3., -4., 8.0, 6.0,0., 3) 

CALL  SYMBOL (-.5,3.4 ,.2,2H2R,0.  ,2) 

CALL  SYMBOL(-.7,3.4*.2,25,0.,-1) 

CALL  SYMBOL(2.4,.35,.2,2H2G,0.,2> 

CALL  SYMBOL(2.2,.35,.2,25,0.,-1) 

IF  ( SPACE  .EQ.  NAME(7>)  GO  TO  720 
CALL  SYMBOL (2. 8*0. 25,. 15, 12 9, 0.,-l) 

CALL  SYMBOL (-.1,3. 3,. 15, 12 9,0. ,-l> 

GO  TO  719 
720  CONTINUE 

CALL  SYMBOL(2. 8, 0.25, .15,108,0.,-!)  _ ... 

CALL  SYMBOL(-.1,3.3,.15,108,0.,-1> 

IFdUNIT  .EQ.  1)  CALL  S  YMB  OL  ( -.  5 , 3.  ,  .  15 , 3HMPA  *  0  .  ,  3  ) 

IFdUNIT  .NE.  1)  CALL  S  YMBOL  (  -  .  5 , 3  .  ,  .  1 5 , 3HKS  I  ,  0  .  ,  3  > 

719  CONTINUE 
0DX=-2.85 
DDY=-3.7 

IF (CRITRIA  .EQ.  NAME ( 3 )  .OR.  CRITRIA  .EQ.  NAME ( 8 )  )  GO  TO  900 
GO  TO  901 
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900  CALL  SYMBOL (0»8*-3»5*0»15*4HF  =,0.,4) 

CALL  NUMBER (1.5 ,-3. 5*0. 15  ,F SI 2*0. *1) 

CALL  SYMBOL< 0.95* -3 .6* 0 • 12* 2HXY, 0. *  2 > 

CALL  SYMBOL <1.05* -3 .3*0. 12*11*0. *-l) 

901  CONTINUE 

IFtCRITRIA  .EQ.  NAME(3) )  CALL  S YMBOL ( DDX * DDY *  0 . 1 5 *  1 6HTS A  I  WU  CRITE 
*RIA,0.,16> 

IFCCRITRIA  .EQ.  NAME(4))  CALL  SYMBOL ( DDX , DDY, . 15, 1 OHMAX  STRESS, 0.* 
*10) 

IF(CRITRIA  .EQ.  N AME ( 5 ) >  CALL  S YMBOL C DDX , DDY * . 1 5 , 1  OHM AX  STRAIN, 0., 
*10) 

IF(CRITRIA  .EQ.  NAME<8)>  CALL  S YMBOL ( DDX , DD Y *  0 . 1 5, 15HCH AM  I S  CRITER 
*IA,0.,15> 

I F ( CR I TR I  A  .EQ.  NAME ( 9 ) )  CALL  S YMBOL ( DDX , DD Y , 0 . 1 5  *  1 6HH0FFMAN  CRITE 
*RIA,0.,16) 

IF (CRITRIA  .EQ.  NAME(IO))  CALL  S YMB OL ( DDX , DD Y , 0 . 15 » 1 6HH  ILL  CR I  TER  I 
* A  , 0 • , 16) 

721  CONTINUE 

CALL  PLT(STX1,STY1,STX2,STY2,STX3,STY3,SCX1,SCY1,SCX2,SCY2,SCX3, 
*SCY3,NPT ,TH*9) 

722  CONTINUE 

IF ( MULT  I  .EQ.  NAME ( 1 4  )  )  60  TO  729 
IFCSMBL  .EQ.  NAME ( 13  )  )  GO  TO  728 
GO  TO  729 

728  CALL  PLOT (13., 0., 3) 

CALL  PLOT (13., 0., -2)  , 

729  CONTINUE 

IF(PARAM3  .EQ.  NAME (11))  GO  TO  717 
IF(PARAM3.NE.NAME(1) )  GO  TO  6 
717  CONTINUE 

CALL  PLOTE(N) 

6  CONTINUE 

WRITE  11 

11  FORMAT (5X,*FAILURE  SURFACE  FOR  THIS  LAMINATE  HAS  BEEN  PLOTTED*) 
RETURN 
END 


SUBROUTINE  STRNG(TH,HS,PM,NLDCN) 

INTEGER  CRITRIA, SPACE, PLANE 

DIMENSION  Q(3*3)»HL(40),TH(40),SFLN(40),SFUN(40) 

DIMENSION  GB ( 3 ) *6(3,3)* SFLL (40) ,SFU(40)*EN1(3)*END(3),GQ(3> 
DIMENSION  SB3(3,3)*HT(40)*TL(40),TT(40),P(3) ,HLC (40) 

DIMENSION  SXC(40) »SXL(40),SYL(40),SYC(40),SXLT (40) 

DIMENSION  ES11(40),ES22(40),VS12(40)»GS12(40)»SB1(3*3)»SB2(3»3) 
DIMENSION  DU2(3,3),DU4(3,3),DU5{3*3),E1(3),E2(3),E3(3)*E4(3) 
DIMENSION  ETR(3),ENT(3),U(5) *DU7 ( 3, 3 ) , HS ( 4 1 ) , AN (3), AM (3) 
DIMENSION  NAMES(3),ETR2(3) 

COMMON  /BLK2/  ES 1 1 ,E S2 2 ,VS 1 2 , GS 1 2 ,NNL 
COMMON /BLK/HLC,HT,TL,TT 
C0MM0N/SBS/SB1,SB2»SB3,P»QQ,DT,C 
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COMMON  /BLKl/  SXL,SXC,SYL,SYC,SXLT,HL 

COMMON  /UNIT/  IUNIT 

COMMON  /BLH/  H,HI 

COMMON  /BLDU/  DU2 ,DU4 * DU5 , DU7 

COMMON  /SF/ SFL, COUNT 

COMMON  /CRIT/  C R I TR I A , S PACE  *  PL ANE  ,F S 1 2 
NAMELIST  /STRESS/  AN, AM 

DATA  NAMES/7HTSAI  WU , 9HM AX S TR A  IN , 9H M A XSTRES S/ 

C0UNT=2. 

READ  700, CRITRIA, SPACE, PLANE, SMBL,FS12 
IF ( NLDCN  .EQ.  1)  READ  STRESS 
700  FORMAT (4A10*F10.3) 

DO  601  ICND=1, NLDCN 
IF(NLDCN.EQ.l)  GO  TO  602 
READ  STRESS 

602  CONTINUE 
WRITE  13 
WRITE  603 

603  FORMAT (6X»*APPLIED  MECHANICAL  STRESS  AND  MOMENT  RESULTANTS  N,M*> 
CALL  WRT1 (AN) 

CALL  WRTKAM) 

WRITE  13 
DO  600  K= 1 ,NNL 

IF(FS11(K)  • E Q •  0.0)  GO  TO  600 
DO  109  1=1,3 
SUM  1 =0 • 

SUM2=0 • 

SUM3=0. 

SUM4  =  0 • 

DO  209  J=l,3 

SUM1=SUM1+DU7(I,J)*AN< J)-DU5<I* J)*AM( J)+HSCK)*<-DU2( I,J)*ANC J>+ 
*DU4  < I »  J ) *  AM ( J ) > 

SUM2=SUM2+DU7(I,J)*AN(J)-DU5(I,J)*AM(J>+HS(K+1)*(-DU2(I,J)*AN(J)+ 

*DU4<I,J)*AM(J>> 

SUM3=SUM3+DU7< I, J>*  P ( J ) -DU 5 C I ,U ) *Q Q ( J ) +HS C K ) * ( -DU2 C I , J ) *  PfJ>  + 
*DU4(I,J)*QQ(J)) 

SUM4=SUM4+DU7( I ,J>*  P C J ) -DU5 ( I , J ) *G G C J ) +  HS C K+ 1 ) * ( -DU2 C I , J ) *  P( J)+ 
*DU4(I,J>*QG(J>> 

209  CONTINUE 

El < I) =SUM1 
E2(I)=SUM2 
E3 ( I ) =SUM3 
E  4 < I)  =  SUM4 
109  CONTINUE 

C  E 1 <  J )  =MECHANICAL  STRAIN  COMPONENTS  LOWER  SURFACE 

C  E  2  <  J )  =MECHANICAL  STRAIN  COMPONENTS  UPPER  SURFACE 

C  E4 ( J)  =NONMECHANICAL  STRAIN  COMPONENTS  UPPER  SURFACE 

C  E3(J>  =NONMECHANICAL  STRAIN  COMPONENTS  LOWER  SURFACE 

WRITE  4494 

4494  FORMAT (5X**EFFECTIVE  NONMECHANICAL  STRAIN  COMP.  AT  LOWER  SURF.*) 

CALL  WRT(E3,5) 

WRITE  13 
WRITE  18,TH(K> 

13  FORMAT (5X  ,* - 

★ - *  ) 

18  FORMAT  <5X,*0N  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN.  =*,F5.1> 
WRITE  19 

19  FORMAT (5X,*MECH/NONMECH  LOWER  SURF.  UPPER  SURF.  *) 

WRITE  20 

20  F0RMATC5X,*  MECH.  1  2  *) 
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21 


WRITE  21 

FORMAT (5X»*  NONMECH.  3  4  *> 

WRITE  13 
THT=TH(K) 

CALL  TRE<E1,ETR*THT) 

CALL  TRE(E3,ENT*THT) 

CALL  NMSN<C*DT*EN1*HLC*HT,TL*TT*K> 

CALL  VOI (ENT,EN1,END) 

CALL  MODULS<ESll(K) *ES22<K) ,VS12<K) , G SI  2 < K > * G* U ) 

CALL  FAILCO(EltQtSXL(K) *SXC(K) ,SYLCK) *SYC(K) tSXLT(K) *6tGB) 
CALL  FSFTY<ETR,END,G*GB,SM,SMN) 

WRITE  701 
WRITE  13 
CALL  WRT(ETR*1) 

CALL  MVM  <  Q*ETR  *E1 ) 

SFLL (K) =  SM 
SFLN(K)=SM/<H*PM) 

CALL  TRECE2,ETR2,THT) 

CALL  WRT <  ETR  2  *  2 ) 

CALL  MVM ( Q*  ETR2*E2) 

CALL  WRT (  EN  D  *  3 ) 

CALL  MVM ( Q*END*E3 ) 

CALL  TRE(E4,ENT*THT> 

CALL  VDI <ENT*EN1*END> 

CALL  WRT ( END*  4) 

CALL  MVM (Q*END*E4) 

701  FORMAT  (10X**STRAIN  COMPONENTS*) 

WRITE  13 

703  FORMAT  <10X,*STRESS  COMPONENTS*) 

WRITE  703 
WRITE  13 
CALL  WRT (El  *  1 ) 

CALL  WRT(E2*2) 

CALL  WRT (E3*3) 

CALL  WR  T  <  E4  *  4 ) 

CALL  FSFTY(ETR2*END*G*GB»SM*SMN) 

SFUCK)=SM 

SFUN(K)=SM/(H*PM) 

WRITE  13 

IFCCRITRIA  .EQ.  N  AM  E  S ( 2  ) )  GO  TO  502 
IF (CRITRIA  .EQ.  NAMES( 3) )  GO  TO  503 
501  WRITE  510, CRITRIA 


510  FORMAT(5X*A10**FAILURE  CRITERIA*) 
WRITE  606,SFLLCK) ,SFUCK) 


606 

FORMAT(5X**STRENGTH  RATIO  R: 
*E 1 0 • 3 ) 

LOWER 

SURF.=*,E10.3,* 

UPPER 

SURF •  =  *  * 

IF ( IUNIT  .EQ.l)  GO  TO  4492 
WRITE  4491*SFLN<K)*SFUNCK) 

GO  TO  4493 

4492 

WRITE  449*SFLN(K) *SFUN(K) 

4493 

CONTINUE 

449 

FORMAT (5X,*STRENGTH  R#  MPA  I 
*F1 0.3) 

LOWER 

SURF .=**F10.3** 

UPPER 

SURF •  =  *  * 

4491 

F0RMAT(5X,*STRENGTH  R#  KSI  : 

LOWER 

SURF.=**F10.3,* 

UPPER 

SURF •=** 

*F  1 0 • 3 ) 

WRITE  13 
GO  TO  504 
502  WRITE  505 

505  FORMAT< 5X,*MAXIMUM  STRAIN  FAILURE  CRITERIA*) 
WRITE  506 
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506  FORMAT  ( 5X  **M  AX  APPLIED  STRESS  RESULTANTS/  AVERAGE  STRESSES*) 
WRITE  13 

WRITE  507 

507  FORMAT (5X**L0WER  SURFACE*) 

CALL  MXSTRN(KtETRtAN*Q,PM) 

WRITE  508 

508  FORMAT (5X**UPPER  SURFACE*) 

CALL  MXSTRN(K*ETR2*AN*Q*PM) 

GO  TO  5 0 A 

503  WRITE  509 

509  F0RMAT(5X**MAXIMUM  STRESS  FAILURE  CRITERIA*) 

WRITE  506 

WRITE  507 

CALL  MXSTRS<KfEl*ANtPM) 

WRITE  508 

CALL  MXSTRS(K*E2*AN*PM) 

5  0 A  CONTINUE 
WRITE  702 

702  FORMAT (1H0//) 

600  CONTINUE 

601  CONTINUE 
RETURN 
END 


SUBROUTINE  SYM(X) 
DIMENSION  X  (  3  *  3  ) 
X(3,1)=X<1*3) 

X (3*2)=X<2«3) 

XC2*1)=X(1*2) 

RETURN 

END 


SUBROUTINE  S YMBL C X « Y « TH * PLNM, N « I ) 

DIMENSION  TH ( 2  0 ) *  PLNM (20) 

CALL  SYMBOL (X,Y». 15f AHNNL=»0.,A) 

X 1 =X+ • 6 
HNL  =  N 

CALL  NUMBERCXl*Y,.15tHNLt0.*-l) 

X2=X-.25 

Yl=Y-.25 

CALL  SYMBOL(X2,Y1*.15,1AA,0.*-1) 

CALL  SYMBOL(Xl*Yl**15t5HPLIES»0**5) 

XSl=X-.25 

XS2=X+.8 

YSl=Y-.55 

NHL=N/2 

DO  1  INL=1»NHL 

I F  < I  • EQ .  0  .AND.  INL.  GT.  2)  GO  TO  2 

IF  <  TH  < INL )  .GE.  0.)  CALL  NUMBER ( XS1 ♦ Y SI , . 1 5  * TH C I NL ) ♦ 0 . *  1 ) 
XS3=XS1-.15 

IFCTHCINL)  .LT.  0.)  CALL  NUMBER ( XS3 * YS 1  * .15 *TH C I NL ) * 0 . 1 1 ) 
2  CONTINUE 

60 


IF ( I  .FQ.  0  • AND*  INL.  6T.  2>  CALL  SYMBOL C XS1 * YS1 *. 15* 1 40  *  0 . *-l > 
I F { I  .EG.  0  .AND.  INL  .EG.  4)  CALL  S YMB OL ( XS3 * YS 1 * . 15 *  1 H- *  0 . *  1 ) 
CALL  NUMBER  CXS2* Y SI *.15, PLNM< INL)* 0.*-l) 

YSl=YSl-.3 

CONTINUE 

RETURN 

END 


SUBROUTINE  TRE  <Z*T*THE  > 
DIMENSION  Z(3)*TC3) 
TH2=2.*THE 
C=C0SM(TH2) 

S=SINM(TH2> 

P=0.5*(Z(1)+Z(2)> 

Q=0.5*(Z(1)-Z<2>> 

R=0.5*Z<3> 

T(1 )=P+Q*C+R*S 

T  <2  )=P-Q*C-R*S 

T<3)=-2.*<Q*S-R*C> 

RETURN 

END 


SUBROUTINE  TRS ( Z  *  T*  THE) 
DIMENSION  Z(3)*TC3) 
TH2=2.*THE 
C=COSM(TH2) 

S=SINM<TH2> 
P=0.5*(Z(1)+Z<2)> 
Q=0.5*<Z(1)-Z<2) > 

R=Z<3> 

T(1)=P+Q*C+R*S 

T(2)=P-Q*C-R*S 

T(3  )=-Q*S+R*C 

RETURN 

END 


SUBROUTINE  US(Q*U) 

DIMENSION  Q(3,3)*U<5) 

SQ12=Q(1*1)*Q<2*2) 

SQ24=2.*Q(1*2)+4.*GC3*3> 

U( l)=(3.*SQ12+SG24)/8. 

U(2)=(G(1*1 >-Q(2,2>)/2. 
U<3)=(SQ12-SQ24)/8. 

U(4)=<SQ12+6.*Q(l*2)-4.*Q<3,3))/8. 

U(5>=(SQ12-2.*Q(l*2)+4.*QC3,3>>/8. 

RETURN 

END 
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SUBROUTINE  VDI(A*B*C) 
DIMENSION  A (3 )*B (3) *C (3) 
DO  1  1=1*3 
CCI)=A<I>— BCI) 

1  CONTINUE 

RETURN 
END 


SUBROUTINE  VNF(V.VN) 
DIMENSION  V(3),VNC3> 
DO  1  1=1*3 
VN(I)=-V< I) 

1  CONTINUE 

RETURN 
END 


SUBROUTINE  WITECQ*A) 

DIMENSION  Q<3*3)*A<3*3> 

DO  1  1=1*3 

WRITE  2,(Q(I»J)*J=1*3)*(A(I*J)*J=1*3) 
2  F0RMAT(2X*3E11»3*1X*3E11«3> 

1  CONTINUE 

RETURN 
END 


SUBROUTINE  UITE1(X,Y) 

DIMENSION  X(3*3)*YC3*3) 

DO  1  1=1*3 

WRITE2*(X<I*J)*J=lf3)*(Y(I*K)*K=l*3) 

1  CONTINUE 

2  FORMAT<2X,6F10.1> 

RETURN 

END 


SUBROUTINE  URITECQ) 
DIMENSION  Q ( 3 * 3 ) 

DO  1  1=1*3 

WRITE2*  CQ(I* J> *  J  =  1  *  3 ) 

2  FORMAT ( 10X*3C10X*E16»3> ) 

1  CONTINUE 

RETURN 
END 
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SUBROUTINE  WRT<X*INDX) 
DIMENSION  X  <  3 ) 

IF< INDX.GT.4)  60  TO  1 
WRITE2*INDX*(X<I)*I=1*3) 

2  FORMAT C10X*I4»3E16»3) 
RETURN 

1  URITE3+ (X(I) *1=1*3) 

3  FORMAT(10X*3E16.3) 

RETURN 

END 


SUBROUTINE  WRT1CX) 
DIMENSION  X ( 3 ) 
WRITE2*(X(I)*I=1*3) 
2  FORMATC 10X» 3F16.3) 

RETURN 
END 


FUNCTION  COSM ( ALPHA  > 

DATA  CONV/. 017453292519943/ 

IF ( ALPHA  .EG.  90.0  .OR.  ALPHA  .EQ.  270.)  60  TO  10 
IF ( ALPHA  .EQ.-90.0  *0R.  ALPHA  .EQ.-270.)  GO  TO  10 
COSM=COS(CONV*ALPHA) 

RETURN 

10  C  0  S  M  =  0  • 

RETURN 

END 


10 


FUNCTION  SINM(ALPHA) 

DATA  CONV/. 017453292519943/ 

IF  C  ALPHA  .EQ.  180.  .OR.  ALPHA 
IFtALPHA  .EQ«-*180»  .OR.  ALPHA 
S I NM=SIN( CONV* ALPHA) 

RETURN 

SINM=0. 

RETURN 

END 


.EQ.  360.)  GO  TO  10 
.EQ.-360.)  GO  TO  10 


FUNCTION  VVM ( X* Y ) 
DIMENSION  X(3)*Y<3) 
SUM  =  0  •  0 
DO  1  1=1*3 
SUM=SUM+X(I)*Y(I) 

1  CONTINUE 

V VM=SUM 
RETURN 
END 
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USER'S  INSTRUCTIONS 
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TRANSFORMATIONS 


Objective 


Stress  Trans. 


Strain  Trans. 


Modulus  Trans. 


Input 

Card  # 

Data 

Format 

1 

TRANSFORMb  STRESS 

(2A10) 

2 

°1'  a2'  °6'  6° 

(4F10.3) 

1 

TRANSFORMb  STRAIN 

(2A10) 

2 

^1 1  ^  g  t 

(4F10.3) 

1 

TRANSFORMb  MODULUS 

(2A10) 

2 

Material+  Units* 

(2A10) 

3 

0° 

(5F10.3) 

Compliance  Trans .  1 

2 

3 

Modulus  &  Compliance 

Trans .  1 

2 

3 


TRANSFORMb  COMPLIANCE  (2A10) 


Material+  Units*  (2A10) 

0°  (5F10.3) 

TRANSFORMb  MODCOM  (2A10) 

Material+  Units*  (2A10) 

0°  (5F10.3) 


*Units  -  either  SI  or  English 

If  the  name  of  the  material  being  used  is  not  given  in  Table  1, 

give  NEW  and  .replace  the  card  #3  by  E  ,  E  ,  v  ,  E  ,  9° .  E  ,  E 

and  E  are  in  GPa.  x  y  x  s  _x  y 

s 

b  denotes  a  blank  column  between  two  commands. 

Input  commands  are  punched  strictly  according  to  the  format.  The 
first  letter^  starting  from  the  first  column  of  the  assigned  columns. 
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Material  Properties 


The  material  properties  of  seven  commonly  used  materials, 
given  in  Table  1,  are  included  in  the  program.  These  properties 
can  be  used  by  providing  the  corresponding  material  name  for  pure 
laminates  and  the  material  property  identification  number  for 
hybrid  laminates.  If  the  problem  under  investigation  requires 
material  properties  different  from  those  given  in  Table  1,  the 
existing  values  can  be  replaced  through  appropriate  input  data. 
The  following  data  cards  are  used: 

(1)  NEWMTRLS  Format  (AlO) 

(2)  b$LAMDATAb  NNM=  ,  EX=  ,  EY=  ...  etc.  such  that 


NNM 

EX 

EY 

VX 

ES 

ALFX 

ALFY 

BTAX 

BTAY 

X 

XD 

Y 

YD 

S 

SH 


=  Number  of  materials  being  read  in 


E  for  each  material 

X 


=  E 


x 


=  E, 


=  a . 


=  a 


=  3 

=  3, 


x 


X 

X' 

Y 

Y« 

S 

h 


-do- 

-do- 

-do- 

-do- 

-do- 

-do- 

-do- 

-do- 

-do- 

-do- 

-do- 

-do- 


-do- 

o 

For  the  use  of  these  properties  in  the  mechanics  problems,  the 
corresponding  names  or  material  property  identification  numbers 
given  in  Table  1  are  used.  The  units  of  these  input  quantities 
are  the  same  as  those  of  Table  1,  i.e.  SI  units.  The  output  results 
can  be  obtained  either  in  SI  units  or  in  English  units,  by  providing 
suitable  instructions . 
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Laminate  Analysis 


On  the  basis  of  lamination  theory,  the  effective  material 
properties  and  strength  characteristics  of  multidirectional  lam¬ 
inates  are  studied  by  using  the  following  input  data: 

(1)  LAMINATE  bb  TYPE  Format  (4A10) 

(2)  CASE  PARAM  PARAMl  PARAM2  P ARAM 3 

Format  (5A10) 

The  various  choices  for  the  parameters  -  TYPE,  CASE, 
PARAM,  PARAMl-3  are  given  below: 


TYPE 


INPLANE 

GENERAL 


DIMFNSIONL 

NORMALIZED 

PARAM 

BOTH 

(Effective 

Modulus) 

ENGCONST 

ALL 

Blank  for  none 

DIMENS IONL 

PARAMl 

NORMALIZED 

(Effective 

Compliance) 

BOTH 

Blank  for  none 
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STRENGTH 

STRNGTHPLT 

PARAM2 

PLT START 

(Strength) 

PLTEND 

PLTONE 

Blank  for  none 

ENGCPLT 

PLTSTART 

PARAM3 

(Engg.  Constants) 

PLTEND 

PLTONE 

Blank  for  none 

(3)  Title  of  the  problem 


Format  (8A10) 


(4)  If  the  CASE  is  pure  then: 

(a)  MATERIAL*  UNITS+  (2A10) 

*name  of  the  material  from  Table  1 
+SI  or  English 

(b)  b$LAYERbNNL=  ,  TH=  ,  PLNM=  , 

DT=  ,  C=  ,  NLDCN=  $ 


NNL 

TH 

PLNM 

DT 

C 

NLDLN 


Total  number  of  layers 

8^,  02 ,  ...  0NNLf  bottom  ply  is  the  first  ply 

Number  of  plies  corresponding  to  each  ply 
orientation 

Temperature  difference,  AT 
Moisture  content 

Number  of  loading  conditions  considered 


In  all  cases  when  PARAM2  is  not  STRENGTH,  NLDCN  is  not  active  and 
is  given  equal  to  1.  When  P ARAM 2  is  STRENGTH,  NLDCN  is  active 
and  will  require  like  number  of  STRESS  input  cards  following  card 
#  ( 5)  . 
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If  the  case  is  HYBRID  then: 


(c)  b*$LAYERSbNNL= 

,  LMPI= 

,  TH= 

PLNM= 

,  IUNIT = 

DT= 

,  c= 

,  NLDCN= 

bb$ 

NNL  =  Total  number  of  layers 

LMPI  =  Layer  material  property  identification 

number  (Table  1) 

TH  =  ®i'  ®2'  ®3  ***  ®NNL  starting  from  bottom 

I  UNIT  =  1  for  SI  and  2  for  ENGLISH  units 

PLNM  =  Number  of  plies  corresponding  to  each 
orientation 

DT  =  Temperature  difference 

C  =  Moisture  content 

NLDCN  =  Number  of  loading  conditions  considered 


(3)  CRITERIA  SPACE  PLANE  SMBL  FS12  MULTICURV  FCTRS  FCTRF 
(IF  APPLICABLE)  (4A10,  F10.3,  A10,  2F10.3) 


TSAI  WU 

CHAMIS 

HOFFMAN 

CRITERIA 

HILL 

MAXSTRAIN 

MAXSTRESS 

SPACE 


STRESS 

STRAIN 


PLANE 


QR 

PRINCIPAL 


*b  denotes  a  blank  column. 


FCTRF  remains  the  same  for  all  cases  in  a 


computer  run. 
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SMBL 


SUPERPOSE 

YES 


FS12: 


Quadratic  polynomial  failure  criteria  interaction 
term. 


MULTICURV:  If  the  failure  surfaces  following  the  failure  surfaces 

being  generated  are  going  to  be  drawn  in  the  same 
figure,  use  this  command.  Otherwise,  leave  blank. 

FCTRS :  Factor  by  which  the  scale  for  strength  plot  to  be 

altered . 


FCTRF :  Factor  by  which  the  figure  size  in  the  strength  plot 

to  be  altered. 

(6)  b$STRESSb  AN  =  N1 ,  N2 ,  Ng,  AM  =  ,  M 2,  Mgbb$ 

where  AN(I),  I  =  1,  2,  3  =  Mechanical  stress  resultants 
AM (I) ,  1=1,  2,  3  =  Mechanical  moment  resultants 

(7)  THEEND  Format (A10) 

This  card  denotes  the  end  of  input  instructions . 

Description  of  Each  Command 

TYPE:  What  type  of  laminate  is  being  considered,  INPLANE  or 

GENERAL?  If  only  inplane  analysis  is  conducted,  one 
can  still  use  general  option,  but  some  of  the  outputs 
may  not  be  of  interest.  Also  because  in  the  inplane 
case  stacking  sequence  does  not  matter,  the  input 
instructions  become  a  bit  simpler. 


CASE:  Which  case  of  the  laminate  material  system  is  under 

consideration,  PURE  or  HYBRID?  For  a  pure  laminate 
the  data  card  stating  the  name  of  the  material  is 
sufficient,  whereas  in  hybrid  laminates  the  layer 
material  property  identification  number  (LMPI)  for 
each  ply  orientation  is  required.  This  parameter 
(LMPI)  is  given  in  $LAYERS  card.  The  description  of 
LAYERS  card  is  given  in  #(4)c,  page  69. 
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PARAM:  What  quantities  in  the  effective  modulus  of  the  lam¬ 

inate  are  required? 

DIMENSIONL:  Dimensional  effective  modulus  matrices  A,  B, 

and  D. 

Normalized  effective  modulus  matrices  A*1,  B* , 
and  D* . 

Both  dimensional  and  normalized  modulus 
matrices  A,  B,  D  and  A*,  B*,  D* . 

Effective  inplane  engineering  constants  of 

the  laminate  E? ,  E2 ,  vS, ,  E?  and  effective 

I  z  Zi  o  f  f  f  f 

flexural  engineering  constants  E-^ ,  E2>  v2i'  E6 

All  the  aforementioned  quantities. 

If  you  don't  want  any  of  these  quantities 
leave  this  column  blank. 

PARAMl :  What  quantities  in  the  effective  compliance  of  the 

composite  laminate  are  required? 

DIMENSIONL:  Dimensional  effective  compliance  matrices  a, 

6,  5. 

NORMALIZED:  Normalized  effective  compliance  matrices  a*, 

6*,  6*. 

BOTH:  Both  dimensional  and  normalized  compliance 

matrices  a,  3,  6  and  a*,  3*/  6*. 

PARAM2:  This  command  controls  the  strength  predictions  and 

plotting  the  failure  surface. 

STRENGTH:  Computes  the  strength  ratio  R  for  each  ply 

orientation  and  corresponding  parameter  R/h 


NORMALIZED: 

BOTH: 

ENGCONST : 

ALL: 

Blank: 
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STRNGTHPLT : 


PLTSTART : 


PLTEND : 


PLTONE : 


Blank: 


for  a  given  loading  condition.  These  parameter 
are  calculated  on  the  basis  of  six  different 
failure  criteria  viz.  Tsai  Wu,  Chamis,  Hoffman, 
Hill,  maximum  stress  (MAXSTRESS)  and  maximum 
strain  (MAXSTRAIN)  failure  criteria.  The 
criteria  to  be  used  has  to  be  assigned  at  an 
appropriate  place  in  the  input  data. 

This  command  generates  failure  surfaces  for 
each  ply  orientation  of  the  laminate.  The 
program  is  set  such  that  it  will  give  failure 
surfaces  for  the  first  three  ply  orientations, 
or  all  the  ply  orientations.  For  obtaining 
failure  envelopes  for  all  the  plies  in  tne 
laminate  if  the  number  of  plies  is  more  than  6, 
the  SMBL  parameter  in  card  number  5  is  SUPER¬ 
IMPOSE,  otherwise  blank.  There  is  a  choice  of 
failure  criteria  space,  plane  and  the  failure 
criteria  interaction  term  F*^,  as  given  in 
card  5 . 

If  the  plot  assigned  by  this  card  is  the  first  / 
strength  (engineering  constant)  plot  in  the 
Job  P ARAM 3  (P  ARAM  2)  should  be  PLTSTART. 

If  the  plot  assigned  by  this  card  is  the  last 
strength  (engineering  constant)  plot  in  the 
Job  P ARAM 3  (P ARAM 2)  should  be  PLTEND. 

If  there  is  only  one  strength  (engineering 
constant)  plot  in  the  Job  PARAM3  (PARAM2) 
should  be  PLTONE . 

If  no  strength/  effective  engineering  constants 
plot  is  required  leave  these  columns  blank. 
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PARAM3:  This  command  controls  the  plots  of  effective  engi¬ 

neering  constants  versus  the  angle  ply  <j>  in  the  lam¬ 
inate  (0  /8  /+d>  )  . 

n/  p' -Yq' s 

ENGCPLT:  This  command  will  give  the  plot  of  the  effec¬ 

tive  engineering  constants  E|,  E E£  and 

PLTSTART 

PLTEND  :  As  mentioned  in  PARAM2 . 

PLTONE 

I 

Blank:  If  no  effective  engineering  constants/strength 

plot  is  required. 

CRITERIA:  There  is  a  choice  of  failure  theories  that 

can  be  used  for  finding  the  successive  ply 
strengths  of  the  laminate.  Six  failure  theo¬ 
ries  [9]  have  been  included  in  the  program. 

In  the  quadratic  polynomial  failure  criteria, 
the  hygrothermal  effects  are  included  whereas  in 
the  Max.  stress  and  the  Max.  strain  failure  theories 
these  effects  are  not  included. 

SPACE:  The  failure  surfaces  can  be  obtained  either 

in  stress  space  or  in  strain  space. 

PLANE:  Stress  space  and  strain  spaces  can  be  repre¬ 

sented  either  in  principal  stress  or  strain 
loadings  or  in  qr  plane  loading  in  which: 

Pe  =  Ce  1  +  e2)/2  =  0 

q£  =  (ex  -  e2)/2  '  re  =  e6//2 

Pa  =  Ccr1  +  a2)  /2  =  0 

ga  =  (al  “  a2)/2 
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By  assigning  the  PLANE  parameter  as  QR  a 
strength  plot  in  qr  space  can  be  obtained. 
Otherwise  the  strength  plot  will  be  in  princi¬ 
pal  strain  or  stress  space  as  the  case  may  be. 

SMBL :  This  input  parameter  controls  the  plotting  of 

failure  surfaces  w.r.t.  the  number  of  plies 
for  which  the  failure  surfaces  are  plotted. 
SUPERPOSE  will  plot  the  failure  envelopes  of 
all  the  plies  if  the  number  of  plies  is  up  to 
50.  If  the  corresponding  columns  are  left 
blank,  failure  envelopes  for  the  first  three 
layers  will  be  plotted.  The  laminate  layer  data 

should  be  given  such  that  the  first  three  orientations 
are  of  prime  Interest? 

FS12:  The  failure  criteria  interaction  term  (F*  ) 

xy 

in  the  case  of  Tsai  Wu  and  Chamis  failure 
criteria. 

MULTICURV:  In  case  the  failure  envelopes  generated  by 

the  current  data  input  is  going  to  be  plotted 
on  the  same  figure  as  those  generated  by  the 
following  data  set  provide  the  command 
MULTICURV  in  the  corresponding  columns,  other¬ 
wise  leave  these  columns  blank. 


FCTRS :  For  changing  the  scale  of  the  strength  plot, 

i.e.  for  reducing  the  scale  to  1/2,  FCTRS  is 
given  as  0.5. 

FCTRF :  For  changing  the  figure  size  in  the  strength 

plot,  i.e.  for  increasing  the  figure  size  by 
25%,  FCTRF  is  given  as  1.25.  FCTRF  <  1.4. 

^tFor  laminates  with  a  number  of  ply  orientations  less  than  3,  laminate  data 
should  be  arranged  such  that  NNL<s3.  This  precaution  is  taken  to  avoid  the 
use  of  unspecified  values  of  0.,  (i=2,3)  in  case  NNL=1 ,  stored  in  the 
available  spaces. 


74 


ILLUSTRATIONS 
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I.  Transformation  of  in-plane  stress  (a1f  a_,  a r)  through  an 

JL  Z  D 

angle  0 . 

o1  =  10,  MPa,  o 2  =  15.  MPa,  og  =  12.  MPa,  0  =  27.5° 
Input  data: 

TRANSFORM  STRESS 

10.  15.  12.  27.5 

THEEND 


Output : 


STRESS  TRANSFORMATION  THRU 

SIGMA1 


27.5  DEGREES 
SIGMA2  SIGMAS 


GIVEN  10.000  15.000  12.000 

TRANSFORMED  20.856  4.104  6.531 


II.  Transformation  of  in-plane  strain  (e^,  e0,  e^)  through  an 
angle  0 . 

e  =  2.0,  £0  =  .05,  £,  =  1.0,  0  =  15.0 
}  z  o 

Input  data: 


TRANSFORM  STRAIN 

2.0  .05  1.0  15.0 

THEEND 


Output : 


STRAIN  TRANSFORMATION  THRU  15.  DEGREES 

EPSLN1  EPSLN2  EPSLN6 


GIVEN 

TRANSFORMED 


2.000  .050  1.000 

2.119  -.069  -.109 
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III.  Modulus/Compliance  Transformation 

1.  Transformation  of  modulus  matrix  of  T300/5208  material  stored 
in  the  program  through  30°  angle. 

2.  Transformation  of  compliance  matrix  of  AS/3501  material  stored 
in  the  program  through  30°  angle. 

3.  Transformation  of  modulus  matrix  and  compliance  matrix  of 
a  new  material  through  30°  angle. 


E  =  190.  GPa,  E  =12.0  GPa,  v  =  0.3,  E  =8.1  GPa. 
x  y  '  x  s 


Input : 


TRANSFORM 
T3  00/52  08 
30.0 

TRANSFORM 

AS/3501 

30.0 

TRANSFORM 

NEW 

190. 

THEEND 
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Output 


ENGLISH 

UNITS 

V" 

MODULUS 

TRANSFORMATION 

THRU 

30.  DEGREE 

ANGLE 

tl. 

E  +  06  PSD 

MODULUS 

OF  THE  MATERIAL 

TRANSFORMED  MODULUS 

26. A 

.4 

0.0 

15.9 

4.7 

7.9 

.4 

1.5 

0.0 

4.7 

3.4 

2.9 

0.0 

0.0 

1.0 

7.9 

2.9 

5.3 

Case. 


SI  UNITS 

COMPLIANCE  TRANSFORMATION  THRU  30.  DEGREE  ANGLE 

Cl/TPA) 


COMPLIANCE 

OF  MATERIAL 

TRANSFORMED 

COMPLIANCE 

7.2 

-2.2 

0.0 

36.6 

-5.5 

-49.0 

-2.2 

111.6 

0.0 

-5.5 

88.8 

-41.4 

0.0 

0.0 

140.8 

-49.0 

-41.4 

127.6 

ENGLISH  UNITS 

MODULUS  TRANSFORMATION  THRU  30.  DEGREE  ANGLE 

tl.E+06  PSI) 


MODULUS  OF 

THE  MATERIAL 

TRANSFORMED  MODULUS 

27.7 

.5 

0  .0 

16.8 

5.0 

8.2 

.5 

1.8 

0.0 

5.0 

3.8 

3.1 

0.0 

0.0 

1.2 

8.2 

3.1 

5.6 

COMPLIANCE 

TRANSFORMATION  THRU 
(1 .0E- 

30. 

09/PSI) 

DEGREE  ANGLE 

COMPLIANCE 

OF  MATERIAL 

TRANSFORMED  COMPLIANCE 

36.3 

-10.9 

0.0 

211.8 

-51.9 

280.4 

-10.9 

574.6 

0.0 

-51.9 

481.0 

185.8 

0.0 

0.0  851  .2 

-280.4 

-185.8 

687.3 
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IV.  Effective  material  properties  of  laminates  and  strength  of 
sandwich  laminates: 

1.  Modulus  matrix  A  for  a  (0/90)  laminate  -  dimensional. 

~  b 

2.  Modulus  matrix  A*  for  a  (0/90)  laminate  -  normalized. 

—  s 

3.  Compliance  matrix  a  for  a  (0/90)  laminate  -  dimensional. 

~  s 

4.  Compliance  matrix  a*  for  a  (0/90)  laminate  -  normalized. 

<v  S 

5.  Engineering  constants  for  a  (0/90)  laminate. 

s 

6.  Strength  of  a  sandwich  laminate.  One  layer  thickness  of 
sandwich  core  is  equal  to  the  thickness  of  eight  plies. 


Input :  ’ 


LAMINATE  INPLANE 
PURE  DIMENSIONL 

MODULUS  MATRIX  FOR  < 0 /90 > -S YMM .  LAMINATE 
T300/520B  SI 

SLAYER  NNL=4,TH=G.*2*9Q.,  Q.,PLNM=4*1.  *DT=0.* 

C=0.,NLDCN=1  S 
LAMINATE  INPLANE 
PURE  NORMALIZED 

MODULUS  MATRIX  FOR  C0/90)-SYMM.  LAMINATE 
T300/5208  SI 

SLAYER  NNL=4*TH=0.*2*90.t0.»PLNM=4*l. *DT=0.* 

C=0.vNLDCN=l  $ 

LAMINATE  INPLANE 

PURE  DIMENSIONL 

COMPLIANCE  MATRIX  FOR  (0/90)-SYM.  LAMINATE 
T30G/5208  SI 

SLAYER  NNL  =  4*TH  =  0. *2*90. *0. *PLNH  =  4*1.  tDT  =  0 .* 

C=0.*NLDCN=1  S 
LAMINATE  INPLANE 
PURE  NORMALIZED 

COMPLIANCE  MATRIX  FOR  (0/90J-SYM.  LAMINATE 
T300/5208  SI 

SLAYER  NNL=4»TH=0.f2*90.*0.»PLNM=4*l.*DT=0.» 

C=0.»NLDCN=1  $ 

LAMINATE  INPLANE 
PURE  ENGCQNST 

ENGINEERING  CONSTANTS  FOR  C0/90)-SYM.  LAMINATE 
T3  0  0/52  08  SI 

SLAYER  NNL=4*TH=0.*2*90. t 0 . » PLNM=4* 1 . *DT=0.» 

C=0.,NLDCN-1  $ 

LAMINATE  INPLANE 

HYBRID  STRENGTH 

STRENGTH  OF  C 0 / 0 /9 0 /90/4CO RE ) -S YMM.  L  AMI •  TSAI-WU  FXYS=-.5 
SLAYERS  NNL=5,LMPI=2*1,6*2*1*TH=0.,90. *0.*90.*0. * 
PLNM=2*2.*1.*2*2.*IUNIT=1*DT=0.»C=0.*NLDCN=1  $ 

TSAI  WU  STRESS  -0 .5 

SSTRESS  AN=1.*0.*0.*AM=3*0.  $ 

THEEND 
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Output : 


MODULUS  MATRIX  FOR  (0/90>-SYMM.  LAMINATE 


NUMBER  OF  PLYS  =  A 

ANGLES  OF  PLY  ORIENTATION  IN  DEGREES*  BOTTOM  LA YER=1 
0.0  90.0  90.0  0.0 

NUMBER  OF  LAYERS  FOR  CORRESPONDING  PLY  ORIENTATION 

1.0  1.0  1.0  1.0 

MATERIAL  DROPERTY  IDENTIFICATION  NUMBER  FOR  CORR.  PLY 
1  1  1  1. 

TEMPERATURE  DT=  0.0000  MOISTURE=  0.0000 
SI  UNITS 


A  B 


B  D 


480E  +  08 

.145E+G7 

0. 

145E+07 

•  4  8  0  E+  0  8 

0. 

0. 

•358E+07 

291E-10 

0. 

0. 

•637E-11 

0. 

0. 

.909E-12 

*  29  1 E  -1  0 

0. 

0. 

0. 

. 637E-11 

0. 

0.  - 

0. 

. 909E-12 

. 167E+01 

• 302E-01 

0. 

•302E-01 

•  33 1 E  +  0  0 

0. 

0. 

0. 

• 747E-0 1 

MODULUS  MATRIX  FOR  <0/90)-SYMM.  LAMINATE 


NUMBER  OF  PLYS  =  A 

ANGLES  OF  PLY  ORIENTATION  IN  DEGREES*  BOTTOM  LA YER=1 
0.0  90.0  90.0  0.0 

NUMBER  OF  LAYERS  FOR  CORRESPONDING  PLY  ORIENTATION 

1.0  1.0  1.0  1.0 

MATERIAL  PROPERTY  IDENTIFICATION  NUMBER  FOR  CORR.  PLY 
1111 

TEMPERATURE  DT=  0.0000  MOISTURE=  0.0000 
SI  UNITS 


A#  B# 


3B#  D# 

GP  A 

98.1 

2.9 

0.0 

.0 

0.0 

0.0 

2.9 

96.1 

0.0 

0.0 

.0 

0.0 

0.0 

0.0 

7.2 

0.0 

0.0 

.0 

.  0 

0.0 

0.0 

160.4 

2.9 

0.0 

0.0 

.0 

0.0 

2.9 

31.8 

0.0 

0.0 

0.0 

.0 

0.0 

0.0 

7.2 
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COMPLIANCE  MATRIX  FOR  <0/90)-SYM.  LAMINATE 


NUMBFR  OF  PLYS  =  4 

ANGLES  OF  PLY  ORIENTATION  IN  DEGREES,  BOTTOM  LA YER=1 
0.0  90.0  90.0  0.0 

NUMBER  OF  LAYERS  FOR  CORRESPONDING  PLY  ORIENTATION 

1.0  1.0  1.0  1.0 

MATERIAL  PROPERTY  IDENTIFICATION  NUMBER  FOR  CORR.  PLY 
11  1  1 

TEMPERATURE  DT  =  0.0  0  0  0  MOISTURE=  0.0  0  0  0 

SI  UNITS 


ALPHA  BETA 


TRSBTA  DELTA 


•  2  0  8E-0  7 

-.628E-09 

o.- 

-.364E-18 

• 452E-19 

0. 

-.628E-09 

•208E-07 

0. 

•182E-19 

-•402E-18 

0. 

0. 

0. 

•279E-C6 

0. 

0. 

-.340E-17 

-.364E-18 

•182E-19 

0. 

•  60  0E  +  0  0 

-.547E-01 

0. 

•452E-19 

-.4Q2E-18 

0  . 

, -.547E-01 

•3D3E+01 

0. 

0. 

0. 

-.340E-17 

0. 

0. 

•  1 3  4  E+  0  2 

NONMECH. 

STRESS  AND 

MOMENT  RESULTANTS  N * M 

0. 

0. 

0. 

0. 

0. 

0. 

COMPLIANCE 

MATRIX  FOR 

( 0/90 )-SYM. 

LAMINATE 

© 

NUMBER  OF  PLYS  =  4 

ANGLES  OF  PLY  ORIENTATION  IN  DEGREES,  BOTTOM  LA YER=1 
0.0  90.0  90.0  0.0 

NUMBER  OF  LAYERS  FOR  CORRESPONDING  PLY  ORIENTATION 

1.0  1.0  1.0  1.0 

MATERIAL  PROPERTY  IDENTIFICATION  NUMBER  FOR  CORR.  PLY 
1111 

TEMPERATURE  DT=  0  .  0  0  0  0  MOISTUR  E=  0.000  0 
SI  UNITS 


ALPHA# 

TRSBTA# 

BE  TA  # /3 

DELTA# 

Cl./TPA) 

10.4 

-.3 

0.0 

-.0 

.  0 

CD 

• 

o 

-.3 

10.4 

0.0 

.0 

-.0 

0.0 

n  .o 

0.0 

139.5 

0.0 

0.0 

-.0 

-.0 

.  0 

0.0 

6.2 

-.6 

0.0 

.  0 

-.0 

0.0 

“  •  6 

31.5 

0.0 

0.0 

0.0 

-.0 

0.0 

0.0 

139.5 

NONMECH. 

EFFECTIVE 

STRESS 

AND 

MOMENT  N# 

,  M# 

0  . 

0. 

0. 

0. 

0. 

0. 
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ENGINEERING  CONSTANTS  FOR  CP/90)-SYM.  LAMINATE 


NUMBER  OF  PLYS  =  A 

ANGLES  OF  PLY  ORIENTATION  IN  DEGREES*  BOTTOM  LA  YE  R  =  1 
0.0  90. C  90.0  0.0 

NUMBER  OF  LAYERS  FOR  CORRESPONDING  PLY  ORIENTATION 

1.0  1.0  1.0  1.0 

MATERIAL  PROPERTY  IDENTIFICATION  NUMBER  FOR  COP.R.  PLY 
1111. 

TEMFERATURE  DT=  0.0000  MOISTURE=  0.0000 
SI  UNITS 


SOME  ENGINEERING  CONSTANTS,  ES  IN  GP A 

INPLANE  El =  95.991  E2  =  95.991  V21  =  .030  Efe=  7.170 

FLEXURAL:, EF1=160.114  EF2  =  31.727  VF21=  .091  EF6=  7.170 
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STRENGTH  OF  < 0  /  0 /9G/ 9 0/ 4C ORE > - S YMM  .  LAMI.  TSAI-WU  FXYS  =  -.5  (?) 


NUMBER  OF  PL YS  =  .5 

ANGLES  OF  PLY  ORIENTATION  IN  DEGREES*  BOTTOM  L AYER=1 
0.0  Q0.0  0.0  90.0  0.0 

NUMBER  OF  LAYERS  FOR  CORRESPONDING  PLY  ORIENTATION 

2.0  2.0  1.0  2.0  2.0 

MATERIAL  PROPERTY  IDENTIFICATION  NUMBER  FOR  CORR.  PLY 
1  1  6  11 

TEMPERATURE  DT=  0.0000  MOISTURE=  0.0000 
SI  UNITS 


APPLIED  MECHANICAL  STRESS  AND  MOMENT  RESULTANTS  N * M 

1.000  0.000  0.000 

0.000  0.000  0.000 


EFFECTIVE  NONM ECH AN  I C A L  STRAIN  COMP. 

0.  0. 


0. 


ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN.  =  0.0 

MECH/NONMECH  LOWER  SURF.  UPPER  SURF. 


MECH. 

1 

2 

NONMECH. 

3 

4 

STRAIN 

COMPONENTS 

1 

.104E-07 

-.314E-09 

0. 

2 

•104E-07 

-.'314E-09 

0. 

3 

0.  . 

0. 

0. 

A 

0. 

0. 

0. 

STRFSS  COMPONENTS 


1 

♦  1 89E  +  0  A 

.  269E+02 

0. 

2 

•  1 80  E  +  0  4 

*  269E+  02 

0. 

■ 

-Z 

0. 

0. 

0. 

4 

c. 

0. 

0. 

TSAI-WU 

FAILURE  CRITERIA 

STRENGTH 

RATIO  r:  LOWER 

SURF  •  = 

.682E+06 

UPPER 

SURF.r 

•682E+06 

STRENGTH 

R #  MPA  :  LOWER 

SURF.= 

340. 941 

UPPER 

SURF • = 

340.941 

EFFECTIVE 

NONMECHANICAL 

STRAIN 

COMP. 

0. 

0. 

0. 
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ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN, 
MECH/NONMECH  LOWER  SURF.  UPPER  SURF. 


=  90.0 


MECH. 

1 

2 

NONMECH. 

3 

4 

STRAIN 

COMPONENTS 

1 

-.314E-09 

•  1  04E-07 

0. 

2 

-.314E-09 

.  104E-07 

0. 

7 

0. 

0. 

0. 

4 

0. 

0. 

0. 

STRESS 

COMPONENTS 

1 

-.269E+02 

.  1  07E+03 

0. 

2 

-•269E+02 

•  1  07E  +  03 

0. 

3 

0. 

0. 

0. 

4 

0. 

0. 

0. 

TSAI  WU  FAILURE  CRITERIA 

STRENGTH  RATIO  RJ  LOWER  SURF.=  .373E+06  UPPER  SURF.=  .373E+06 
STRENGTH  R«  MPA  :  LOWER  SURF.=  186.698  UPPER  SURF.=  186.698 


EFFECTIVE  NONMECHANICAL  STRAIN  COMP. 

0  •  0  .  0  . 


ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN. 

=  90.0 

MECH/NONMECH 

LOWER  SURF. 

UPPER  SURF. 

MECH. 

1 

2 

NON MECH. 

4 

STRAIN 

COMPONENTS 

1 

-  .  3  1 4  E  -  0  9 

.104E-07 

0. 

2 

-•314E-09 

.104E-07 

0. 

3 

0. 

0. 

0. 

4 

C. 

0. 

0. 

STRESS 

COMPONENTS 

1 

-•269E+02 

•  1  07E+03 

0. 

2 

-.269E+Q2 

•  1 07  E+  0  3 

0. 

3 

0. 

0. 

0. 

4 

0. 

0. 

0. 

TSAI  WU  FAILURE  CRITERIA 

Strength  ratio  r:  lower  surf.=  .373E+06  upper  surf.=  .373E+06 
STRENGTH  R  8  MPA  :  LOWER  SURF.=  186.698  UPPER  SURF.=  186.698 
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EFFECTIVE  NONMECHANICAL  STRAIN  COMP. 


• 

o 

• 

o 

0. 

ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN.  = 

0.0 

MECH/NONMECH 

LOWER  SURF. 

UPPER  SURF. 

MECH. 

1 

2 

NOMMECH. 

3 

4 

STRAIN 

COMPONENTS 

1 

•104E-07 

-*  3 14  E- 09 

0. 

2 

.1 04E-07 

-.314E-09 

0. 

3  • 

0. 

0. 

c. 

4 

0. 

0. 

0  . 

STRESS 

COMPONENTS 

1 

■ 189E+04 

• 269E+02 

0. 

2 

.189E+04 

.  2.69E  +  02 

0  . 

3 

0. 

0. 

0. 

4 

0. 

0. 

0. 

TSAI  WU  FAILURE  CRITERIA 

STRENGTH  RATIO  R :  LOWER  SURF.=  .F82F+06  UPPER  SUPF.=  .682E+06 
STRENGTH  R#  MPA  I  LOWER  SURF.=  340. 941  UPPER  SURF.=  340.941 
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V.  Unsymmetric  Laminates: 

Replace  the  properties  of  first  two  materials  by  new  material 
properties,  calculate  the  effective  modulus  matrices  and  engineering 
constants  and  conduct  the  strength  analysis  of  tvro  unsymmetric 
laminates  in  the  presence  of  curing  stress 

1.  Laminate  (30/-30/30/-30/90/90)  total. 

2.  Laminate  (0/90/45/-45)  total. 

Input  data: 

NEUMTRLS  SI 

SLAMDATA  NNM=2*EX=113  .0 8 *1 37 .9 1 E Y  =  9 .24  *  1 1  .03  * VX=2 * . 3  * 
ES=2*5.86*ALFX=2*-.9f  ALF Y  =  2 *23. 04  * BT AX=2*0 . *BT A Y  =  2*0 . * 

X=2*1448. fXD=1448.*Y=2*51.7*YD=2*207*S=2*93*SH=. 000125  $ 

““LAMINATE  GENERAL 

PURE  ALL BOTH  STRENGTH 

< 3 0/-30/3Q/-3 0/90/90) -TOTAL  LAMINATE 
T300/5208  ENGLISH 

SLAYER  NNL=6,TH=2*90.,-30.»30.*-30.,30.*PLNM=6*1.*DT=-175.* 
C=0.»NLDCN=1  $ 

TSAI  UU  -0.5 

_  IS TRESS  AN=1.*2*0.,AM=3*0.  S 
LAM  IN  ATE  GENERAL 

PURE  ALL  BOTH  STRENGTH 

( 0/90/45/-45)-TOTAL  LAMINATE 
T300/5208  ENGLISH 

SLAYER  NNL=4,TH=-45., 45. *90.» 0. »PLNM=4*1 .,DT=-175.* 

C=0.*NLDC?Y=1  $  - .  . 

TSAI  WU  -0.5 

_  SSTRESS  AN=l.t2*0.*AMr3*0.  S 
THEEND 


Note:  The  stacking  sequence  is  such  that  the  bottom  ply  corresponds 

to  TH(1)  in  the  input  data  card  "$LAYER",  where  as  in 

customary  notation  symbols  (0/90/+45)  ,  0  denotes  top  ply 

—  t  ' 

angle  and  -45  denotes  bottom  ply. 

Output : 
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( 30/-3 0 / 30 /-3 0/9 0/90) -TOTAL  LAMINATE 


© 


NUMBER  OF  PLYS  =  6 

ANGLES  OF  PLY  ORIENTATION  IN  DEGREES,  BOTTOM  LAYER=1 
90.0  90.0-30.0  30.0-30.0  30.0 
NUMBER  OF  LAYERS  FOR  CORRESPONDING  PLY  ORIENTATION 
1.0  1.0  1.0  1.0  1.0  1.0 
MATERIAL  PROPERTY  IDENTIFICATION  NUMBER  FOR  CORR.  PLY 
111111 
TEMPERATURE  DT=  -175.0000  MOISTURE=  0.0000 
ENGLISH  UNITS 


A  B 


B  D 


•214E+06 

.624E+05 

0. 

•855E+03 

• 249E+03 

•231E+03 

•  624  E  +  05 

•  2  1 4  E+  0  6 

0. 

•  249E+  0  3 

- . 135E+04 

•876E+02 

• 

0. 

•756E+05 

•231E+03 

•  876E  +  0  2 

•249E+03 

•855E+03 

•249F+03 

.231E+03 

•127E+02 

•373E+01 

•228E+01 

• 249E+03 

-•135E+04 

•  876E  +  0  2 

•373E+01 

•200E+02 

• 863E+00 

.231F+03 

•876E+02 

•  2  4  9E+  0  3 

•228E+01 

•863E+00 

•  46  8E+  0 1 

A# 

B  # 

3B# 

D# 

l.E+06  PSI 

7.2 

2.1 

0.0 

2.0 

•  6 

.5 

2.1 

7.2 

0.0 

•  6 

-3.1 

.2 

0  .0 

0.0 

2  .6 

.5 

.2 

•  6 

5.9 

1.7 

1.6 

5.9 

1.7 

i.i 

1.7 

-9.3 

.6 

1.7 

9.3 

.4 

1.6 

•  6 

1.7 

1.1 

.4 

C\J 

• 

CM 

SOME  ENGINEERING  CONSTANTS, 

ES  IN  l.E+06 

PSI 

INPLANE  : 

,  E 1 =  6.611 

E  2  = 

6.611  V  2 1  = 

.292 

E6  = 

2.558 

FLEXURAL : 

,  EF 1  =  5.151 

EF2  = 

8.782  VF 2 1  = 

.167 

EF6  = 

1.988 

ALPHA 

BETA 

TRSBTA 

DELTA 

•767E-05 

-•328E-05 

•210E-05 

-•382E-03 

- • 247E- 03 

-•198E-03 

-•328F-05 

•114E-04 

•  4  9  3E  -  0  7 

-.238E-03 

•  86  2  E-0  3 

- • 973E- 0  4 

•21 0E-05 

•493E-07 

•170E-04 

- . 295E-0  3 

-•505E-05 

- • 86 1 E- 0  3 

-•382E-03 

-.238E-03 

-•295E-03 

•  127E  +  0  0 

-  • 33  0  E -0 1 

-. 166E-01 

-•247E-03 

•862E-03 

-•505E-05 

- • 33  OE-O 1 

. 118E+00 

- • 937E-0  2 

-. 198E-03 

-•973E-04 

-•861E-03 

-•166E-01 

r • 93  7E-02 

• 281E+00 

NONMECH. 

STRESS  AND 

MOMENT  RESULTANTS  N, M 

-.362E+02 

362E  +  02 

•119E-11 

•  257E+00 

**  • 

257E+00 

•741E-01 

37 


ALPHA#  BETA#/3 


TRSBTA#  DELTA#  1.0E-09/PSI 


226.7 

-96.9  62. 

2 

-55.6 

-35.9 

-28.8 

-96.9 

337.3  1. 

5 

-34.7 

125.4 

-14.2 

62.2 

1.5  501. 

5 

-42.9 

-.7 

-125.3 

-166.8 

-104.0  -128. 

8 

272.7 

-71.0 

-35.7 

-107.8 

376.3  -2. 

2 

-71.0 

253.8 

-20.1 

-86.3 

-42.5  -375. 

9 

-35.7 

-20.1 

603.3 

NONMECH. 

EFFECTIVE  STRESS 

AND 

MOMENT  NtUM# 

-. 123E+04 

”  • 

123E+04 

.404E-10 

•176E+04 

**  • 

176E+04 

•509E+03 

APPLIED 

MECHANICAL  STRESS 

AND 

MOMENT  RESULTANTS 

N  *  M 

1.000 

0.  000 

0.000 

0.000 

0.000 

0.000 

EFFECTIVE 

NONMECHANICAL  STRAIN 

COMP. 

-. 113E-02 

• 

328E-03 

-•654E-03 

ON  AXIS  STRAIN  /STRESS  COMPONENTS*.  PLY  ORIEN.  = 

90.0 

MECH/NONMECH  LOWER  SURF 

• 

UPPER  SURF 

• 

MECH. 

1 

2 

NONMECH 

.  3 

4 

STRAIN  COMPONENTS 

1 

•367E-06 

. 133E-04 

-•503E-05 

2 

-•848E-06 

•  1 14E-04 

- • 405E- 05 

3 

•240E-03 

.  111E-02 

•  654  E- 03 

4 

-•636E-04 

.  142E-02 

. 5  08  E- 0  3 

STRESS  COMPONENTS 


1  .115E+02 

2  - • 939E+0 1 

3  .442E+04 

A  -.476E+03 


. 181E+02 
•151E+02 
•160E+04 
. 189E+04 


-•427E+01 
-.344E+01 
•  556  E  +  03 
•432E+03 


TSAI  WU  FAILURE  CRITERIA 

STRENGTH  RATIO  Rt  LOWER  SURF.=  .330E+03  UPPER  SURF.=  .366E+03 
STRENGTH  R#  KSI  :  LOWER  SURF.=  11.174  UPPER  SURF.=  12.396 
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EFFECTIVE  NONMECHANICAL  STRAIN  COMP. 

- ♦821E-0?  .239E-D4  -.508E-03 


ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN.  =  90.0 
MECH/NONMECH  LOUER  SURF.  UPPER  SURF. 

MECH.  1  2 

NONMECH.  3  A 


STRAIN 

COMPONENTS 

1 

-•848E-06 

.  11AE-0 A 

-.A05E-05 

2 

-.206E-05 

•955E-05 

-.308E-05 

3 

-•636E-04 

.  1A2E-02 

•508E-03 

A 

-.368E-03 

. 173E-02 

•362E-03 

STRESS  COMPONENTS 


1  -.939E+01  •  151 E  +  0  2 

2  -.302E+02  .121E+02 

3  -• 476E  +  03  .  1 89  E+  0  A 

A  -.537E+04  .218E+0A 


•  3AAE+  01 

•  262  E  +  0 1 
• A32E+03 

•  308E  +  03 


TSAI  WU  FAILURE  CRITERIA 

STRENGTH  PATIO  RJ  LOWER  SURF.=  .366E+03  UPPER  SURF.=  .A15E+03 
STRENGTH  R#  KSI  t  LOWER  SURF.=  12.396  UPPER  SURF.=  1A.0A8 


EFFECTIVE  NONMECH AN  I C AL  STRAIN  COMP. 

-•515F-03  -.280E-03  -.362E-03 


ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN.  =-30.0 
MECH/NONMECH  LOWER  SURF.  UPPER  SURF. 

MECH.  1  2 

NONMECH.  34 


STRAIN  COMPONENTS 


1 

•532E-05 

•217E-05 

•116E-0A 

2 

•A02E-05 

•369E-06 

•  1 05E- 0  A 

3 

-.387E-03 

.  17AE-02 

-•38AE-03 

A 

-.296E-03 

•166E-02 

'  .217E-03 

STRESS 

COMPONENTS 

1 

•  8  87E  +  0  2 

•509E+01 

•986E+01 

2 

•  6  66E  +  0  2 

•213E+01 

•896E+01 

3 

-.569E+04 

•220E+04 

-•327E+03 

A 

-•A22E+0A 

.212E+0A 

.  185E  +  03 

TSAI  WU  FAILURE  CRITERIA 

STRENGTH  RATIO  R*.  LOWER  SURF.= 

. 7  A  7E  +  0  3 

UPPER 

SURF. =  . 

103E  +  04 

STRENGTH  R # 

KSI  :  LOWER  SURF . = 

25.278 

UPPER 

SURF • = 

34.700 

B9 


EFFECTIVE  NONMECHANICAL  STRAIN  COMP. 

-•208E-03  - • 584E-03  -.216E-03 


ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN.  =  30.0 
MECH/NONMECH  LOWER  SURF.  UPPER  SURF. 

MECH.  1  2 

NONMECH.  3  A 


STRAIN  COMPONENTS 


1 

2 

3 

4 


•584E-05 

•371E-05 

-•483E-03 

-.266E-03 


-. 145E-05 
-.241E-05 
. 184E-02 
. 163E-02 


-.843E-05 
- .834E-05 
-.433E-03 
-• 889E- 03 


STRESS  COMPONENTS 


1 

2 

T 

4 


960E+02 

603E+02 

724E+04 

374E+04 


•405E+00 
-• 176E+01 
•225E+04 
.209E+04 


-•717E+01 
- • 7  0  9E+  0 1 
- .368E+03 
-.756F+03 


TSAI  WU  FAILURE  CRITERIA 

STRENGTH  RATIO  Rt  LOWER  SURF.=  .127E+04  UPPER  SURF.=  .164E+04 
STRENGTH  R  8  KSI  :  LOWER  SURF.=  42.843  UPPER  SURF.=  55.621 


EFFECTIVE  NONM ECH AN  I C AL  STRAIN  COMP. 

• 983E- 04  -.888E-03  -.703E-04 


ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN.  =-30.0 


MECH/NONMECH 

LOWER  SURF. 

UPPER  SURF. 

MECH. 

1 

2 

NONMECH  . 

3 

4 

STRAIN  COMPONENTS 


1 

2 

3 

4 


•273E-05 

•144E-05 

-•205E-03 

-.115E-03 


-•144E-05 
-.324E-05 
. 157E-02 
.148E-02 


•947E-05 
•841E-05 
•819E-03 
• 142E-02 


STRESS  COMPONENTS 


1 

2 

3 

4 


.445E+02 
•224E+0  2 
-.276E+04 
-.129E+04 


-.832E+00 
-.379E+01 
•203E+04 
•  1 95  E+  04 


•805E+01 
•715E+01 
•  696  E+  03 
• 121E+04 


TSAI  WU  FAILURE  CRITERIA 

STRENGTH  RATIO  R:  LOWER  SURF.=  .145E+04  UPPER  SURF.=  .198E+04 
STRENGTH  R#  KSI  :  LOWER  SURF.=  49.147  UPPER  SURF.=  67.112 


•755E-04 


EFFECTIVE  NONPECHANICAL  STRAIN  COMP, 

•  4  05F- 0  3  -.119E-02 


ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN.  =  30.0 
MECH/NONMECH  LOWER  SURF.  UPPER  SURF. 


•  MECH. 

1 

2 

NONMECH. 

3 

4 

STRAIN 

COMPONENTS 

1 

•157F-05 

-• 337F-05 

-.825E-05 

2 

-.565E-06 

-•  434E-05 

-.817E-05 

3 

-.492E-04 

. 141E-02 

-•  135E-02 

4 

•168E-03 

•  120  E- 02 

-.180E-02 

STRESS 

COMPONENTS 

1 

•246E+02 

-•392E+01 

-•702E+01 

2 

-.111E+02 

-• 608E+01 

-•694E+01 

3 

-.240E+03 

•  189E+  04 

-•114E+04 

4 

•326E+04 

.169E+04 

-.153E+04 

TSAI  WU  FAILURE  CRITERIA 

STRENGTH  RATIO  R:  LOWER  SURF.=  .202E+04  UPPER  SURF.=  .226E+04 
STRENGTH  R#  KSI  :  LOWER  SURF.=  68.«33  UPPER  SURF.=  76.599 


91 


(  0/90/45/-45  )  -  TOT  AL  LAMINATE! 


f 


NUMBER  OF  PLYS  =  4 

ANGLES  OF  PLY  ORIENTATION  IN  DEGREES,  BOTTOM  L A  YE  R=  1 
-45.0  45.0  90.0  0.0 

NUMBER  OF  LAYERS  FOR  CORRESPONDING  PLY  ORIENTATION 

1.0  1.0  1.0  1.0 

MATERIAL  PROPERTY  IDENTIFICATION  NUMBER  FOR  CORR.  PLY 
1111 

TEMPERATURE  DT  =  -175.0000  MOISTURE^  0.0  0  0  0 

ENGLISH  UNITS 

A  B 

B  D 

•142F+06 

•416E+05  - 

•582E-09  .350E+03 

-  .  166E-»  03 

•  92  0  F+  0  2 

•416F+05 

•  1 4  2  E+  06 

•582E-09  -.166E+03 

-•183E+02 

-92QE+02 

•582E-09 

•582E- 09 

•5G4E+05  •  92  0  E  +  0  2 

•920E+02 

-. 166E+03 

.350E+03 

-•166E+03 

•920E+C2  .642E+01 

•135E+01 

-. 906E+00 

.  166E  +  03 

-.183E+02 

•920E+02  .135E+01 

•279E+01 

-.  906F+00 

.920E+02 

.  9  2  0  E+  0  2  - 

•  1 66  E+  0  3  -.906E+00 

-•906F+00 

. 1  joF+01 

A#  B# 

3B#  D# 

l.E+06 

PSI 

7.2 

2.1 

-.0  1.8 

-.9 

.5 

2.1 

7.2 

.0  -.9 

-.  1 

.5 

-.0 

.0 

2.6  .5 

.5 

-.9 

5.4 

-2.6 

1.4  10.1 

2.1 

-1.4 

-2.6 

-.3 

1.4  2.1 

4.4 

-1.4 

1.4 

1.4 

-2.6  -1.4 

1.4 

2.6 

SOME  ENGINEERING  CONSTANTS,  ES  IN  l.E+06  PSI 

INPLANE  : 

,  E 1 =  6.611 

E  2=  6.611  y  21-  . 

292  E6  = 

2.558 

flexural: 

,  EF  1  =  8.796 

E  F  2  =  3.404  VF  2 1  = 

368  EF6= 

2.036 

ALPHA  BETA 

TRSBTA  DELTA 

• 128E-04  ■ 

-•444E-05  - 

.233E-05  -.111E-02 

. 111E-02 

-•709E-03 

•444E-05 

.9 16E-05  - 

. 132E-05  • .558E-03 

-  .559E-03 

-.401 E- 03 

•233E-05 

-.132E-05 

•308E-04  .155E-03 

-•153E-03 

•  334E-0  2 

•111E-02 

•558E- 03 

•155E-03  .279E+00 

-•166E+00 

•110E+00 

•111E-02  ■ 

-.559E-03  - 

. 153E-03  - .166E+00 

. 562E+  00 

.  173E+00 

•709E-03  ■ 

-  .401E-03 

•334E-02  .110E+00 

. 173E+00 

.  117E+01 

NONMECH.  : 

STRESS  AND  MOMENT  RESULTANTS  N , M 

-•241E+02 

-. 241E+02 

.587E-12 

•428E- 01 

-•428E-01 

•428E-01 
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ALPHA# 


RETA#/3 


TRSBTA# 

DELTA# 

1.0E-09/PSI 

252.3 

-87.4 

-46.0  -72.0 

72.0 

-45.9 

-87.4 

180.4 

-25.9  36.1 

-36.1 

-25.9 

-46.0 

-25.9 

607.5  10.0 

-9.9 

216.2 

-215.9 

108.2 

30.1  177.8 

-106.0 

70.0 

216.1 

-108.4 

-29.7  -106.0 

358.0 

110.2 

-137.6 

-77.8 

648.6  70.0 

110.2 

747.7 

NONMECH. 

EFFECTIVE  STRESS  AND  MOMENT  N# 

♦  M# 

-  •  123E  +  04 

-.  123E+04 

•  298E- 

10 

•661E+03 

-•661E+03 

•661E+03 

APPLIED  MECHANICAL  STRESS  AND  MOMENT  RESULTANTS  N,M 

1.000  0.000  0.000 

0.000  0.000  0.000 


EFFECTIVE  NONMECHANICAL  STRAIN  COMP. 

-•693E-03  .283E-03  -.487E-03 


ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN.  =-45.0 
MECH/NONMECH  LOWER  SURF.  UPPER  SURF. 

MECH.  1  2 

NONMECH.  3  4 


STRAIN  COMPONENTS 


1  .185E-05 

2  .36OE-05 

3  -.492E-04 

4  -.232E-03 


•  651 E- 05 
•476E-05 
. 179E-02 
. 197E-02 


•392E-04 

•282E-04 

-.976E-03 

-•610E-03 


STRESS  COMPONENTS 


1  .333E+02 

•  953  E  +  0 1 

•  333E  +  02 

2  •  6 15E ♦ 0  2 

•  789E  +  0 1 

•240E+02 

3  -.878E+02 

.240E+04 

- • 83  0  E+0  3 

4  -.304E+04 

•  257E+  04 

-.519E+03 

TSAI  WU  FAILURE  CRITERIA 

STRENGTH  RATIO  R:  LOWER  SURF.= 

•279E+03 

UPPER 

SURF  •=  •  357E  +  0  3 

STRENGTH  R#  KSI  :  LOWER  SURF.= 

14.160 

UPPER 

SURF. =  18.138 
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effective:  nonmechanical  strain  comp. 

-•511E-03  .996E-04  -.122E-03 


ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN.  =  45.0 


MECH/NONMECH 

LOWER  SURF. 

UPPER  SURF 

MECH. 

1 

2 

NONMECH. 

3 

4 

STRAIN  COMPONENTS 


1 

2 

3 

4 


•476E-C5 
•  3  0  2  E  -  0  5 
-•354E-03 
-.171E-03 


•  36  0  E- 05 
•535E-05 
•210E-02 
. 191E-02 


-.282E-04 
-.172E-04 
.610E-03 
•244E- 03 


STRESS  COMPONENTS 


1  .801E+02 

2  . 520E+02 

3  -.5C0E+04 

4  -.205E+04 


679E+01 
845E+01 
269E+04 
2  5 1 E+  0  4 


-•240E+02 
-.147E+02 
•519 E +^3 
•208E+03 


TSAI  WU  FAILURE  CRITERIA 

STRENGTH  RATIO  R:  LOWER  SURF.=  .371E+03  UPPER  SURF.=  .450^+03 
STRENGTH  R#  KSI  :  LOWER  SURF.=  18.835  UPPER  SURF.=  2J.P28 


EFFECTIVE  NONMECHANICAL  STRAIN  COMP. 

-.328E-03  - • 833E-C4  .244E-03 


ON  AXIS  STRAIN  /STRESS  COMPONENTS:  PLY  ORIEN.  =  90.0 


MECH/NONMECH 

LOWER  SURF. 

UPPER  SURF 

MECH. 

1 

2 

NONMECH. 

3 

4 

STRAIN  COMPONENTS 


1  -.444E-05 

2  .1 05E-05 

3  -.171E-G3 

4  -.354E-03 


. 128E-04 
. 733  E- 0  5 
•191E-02 
.210E-02 


•233E-05 
• 583E- 05 
-.244E-03 
-•610E-03 


STRESS  COMPONENTS 


1  -.681E+02 

2  ,2  03E  +  02 

3  -.205E+04 

4  -.500E+04 


. 155E+02 
•103E+02 
•251E+04 
.269E+04 


•19RE+01 

.495E+01 

-•208E+03 

-•519E+03 


TSAI  WU  FAILURE  CRITERIA 

STRENGTH  RATIO  R:  LOWER  SURF.=  .297E+03  UPPER  SURF.=  .462E+03 
STRENGTH  R #  KSI  :  LOWER  SURF.=  15.095  UPPER  SURF.=  23.466 


610E-03 


EFFECTIVE  NONMECHANICAL  STRAIN  COMP. 

-•145E-03  - • 266E-03 


ON  AXIS  STRAIN  /STRESS  COMPONENTS*.  PLY  ORIEN 
MECH/NONMECH  LOWER  SURF.  UPPER  SURF. 

MECH.  1  2 

NONMECH.  3  4 


0.0 


STRAIN  COMPONENTS 


1 

.733E-05 

.  105E-05 

-.583E-05 

2 

.185E-05 

•  653E-*  05 

-.932E-05 

3 

-.232E-03 

.  197E-02 

•  61 0E- 03 

4 

-.492E-04 

* 179E-02 

•976E- 03 

STRESS 

COMPONENTS 

1 

•121E+03 

• 438E+0 1 

-.495E+01 

2 

•332E+02 

•957E+01 

-.792E+01 

3 

-.304E+04 

•  257E  +  04 

•519E+03 

4 

-.878E+02 

•240E+04 

•  83  0  E  +  0  3 

TSAI  WU  '  FAILURE  CPITERIA 
STRENGTH  RATIO  R:  LOWER  SURF.= 

•934E+03 

UPPER 

SURF.=  .517E+03 

STRENGTH  R# 

KSI  :  LOWER  SURF . = 

47.422 

UPPER 

SURF.=  26.242 
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VI.  Engineering  constants  plot: 

To  replace  material  properties  of  the  first  two  materials  in  the 
program  and  use  the  first  material  to  compute  effective  modulus 
matrices  for  (0/90/0^) g  laminate,  and  plot  engineering  constants 
for  (0/90/+  <t>2)s  laminates,  $  varies  from  0°  to  90°. 


Input  data: 

NEUMTRLS 

SLAM  DATA  NNM=2,EX=185  .«220.*EY  =  11. 2*10. 5*VX  =  0. 29*0. 31* 
ALFX=12.5tll.,ALFY=-.5*-.3,BTAX=2*O.,BTAY=2*O.»X=140O., 
13OO.»XD=14OO.*13O0.*Y=35.,32*YD=230.*21O.*S=2*75.,SH=2*.125E-3  $ 

LAMINATE  INPLANE 

PURE  DIMENSIONL  PL  TONE  ENGCPLT 

PLOT  FOR  ENGINEERING  CONSTANTS  OF  (0/SQ/F  /-F  )  -  i 

T300/5208  SI  2  2' 

SLAYER  NNL=8*TH=0.*90.*0.*0.*0.*0.*90.*0.*OLNM=2*l.*4fc2., 
2*1.*DT=0.*C=0.*NLDCN=1  $ 

THEEND 


Output : 


PLOT  "OR  ENGINEERING  CONSTANTS  OF  (  Q  /  9  0  /  Fj  /  -  F*  >g  -LAMINATE 
NUMBER  OF  PLYS  =  8 

ANGLES  OF  PLY  ORIENTATION  IN  DEGREES  »  BOTTOM  L  A  Y  E  R  =  1 
0.090.0  G . G  0.0  0 . G  0.090.0  0.0 

NUMBER  OF  LAYERS  FOR  CORRESPONDING  PLY  ORIENTATION 
1.0  1.0  2 . G  2.  0  2.0  2.0  1.0  1  .  C 

MATERIAL  PROPERTY  IDENTIFICATION  NUMBER  FOR  CORR.  PLY 
iiiiiiii 
TEMPERATURE  QT=  0.0000  MOISTURE-  0.0000 
SI  UNITS 


A  B 


b 


•  235  E>09 

•490 E +07 

u  • 

.407r-G9 

•546E-1I 

0  » 

•490 E +07 

.606E+G8 

0  • 

.546r-ll 

•116 E -09 

0  . 

• 

0  . 

.1 08E+08 

0  . 

0  . 

•182E-1Q 

. 4  G  7l -0  9 

.546E-11 

0  . 

.3B4F+C2 

•91SE+00 

0. 

. 116E- 09 

0  . 

•  9 1 8E  +  0 0 

.  1  7  0  E  +  0  2 

0  . 

* 

0. 

• 182E-10 

0  . 

0  . 

.  20  2E+  01 

PLOTS  FOR  ENGINEERING  CONSTANTS  DRAWN 
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Vii.  Failure  envelopes  for  a  pure  or  hybrid  symmetric  laminate  on  the 
basis  of  Tsai  Wu  failure  criteria  in  strain  space,  =  -.5. 

1.  (0/90/+30)  pure  T300/5208  laminate. 

*■“  s 

2.  (0/90/+30)  hybrid  laminate  0°  and  90°  plies  T300/5208  and 

“  s 

+30°  plies  KEVLAR  49. 


Input  data: 

LAMINATE  INPLANE 

pUR£  STRNGTHPLTPLTSTART 

(0/90/30/-30)  SYMMETRIC  LAMINATE  FAILURE  SURFACE  TSAI  UU  FSXY=-l./2* 

®T30  0/52  08 

SLAYER  NNL  =  8*TH  =  0.*90  .*30.*-3  0. »-30 . *30 . *90 • *0 • *PLNM=8*1 • * 
OT=0.»C-0.tNLDCirsi  * 

TSAI  UU  STRAIN  “0*5 

LAMINATE  INPLANE 

hybrid  STRNGTHPLTPLTEND 

/T\  I0/90/30/-30 )  SYM.  LAMINATE  FAIL*  ENVELOPES  TSAI  UU  FSXY=-l/2*  HYBRID* 

VZ/  SLAYERS  NNL=8*LMPI=2*l*4*5*2*l*TH=0.*90**30**2*-30.*30.»90.v0.t 

PLNM-8*!**  IUNIT=  1,  DT=0 . ,  C=0,  NLDCN=1  $ 

TSAI  UU  STRAIN  ‘°*5 

THEEND 


Output : 


(0/90/3C/-30)  SYMMETRIC  '  LAMINATE  FAILURE  SURFACE  TSAI  UU  FSXY=-l./2. 


NUMBER  OF  PLYS  =8 

ANGLES  OF  PL Y  OR IENTA T I  ON  IN  DEGREES*  BOTTOM  LAYER  =  1 
0.0  90.0  30.0-30.0-30.0  30.0  90.0  0.0 

NUMBER  OF  LAYERS  FOR  C  ORR  ESPOND I  N'G  PLY  ORIENTATION 
1.0  1.0  1.0  1.0  '  1*0  1.0  1.0  1.0 
MATERIAL  PROPERTY  IDENTIFICATION  NUMBER  FOR  CORR.  PLY 
11  1  1  1  1  1  1 
TEMPERATURE  DT  =  0.0000  M0ISTURE=  0.000  0 

SI  UNITS 


FAILURE  SURFACE  FOR  THIS  LAMINATE  HAS  BEEN  PLOTTED 


Note:  LMPI  in  input  data  card  number  11  represents  the  material 

property  identification  number  of  each  layer  in  the 
laminate. 
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I 


2 


(0/90/30/-30 )  SYM.  LAMINATE  FAIL.  ENVELOPES  TSAI  WU  FSXY=-l/2.  HYBRID. 


NUMBER  OF  PLYS  =  8 

ANGLES  OF  PLY  ORIENTATION  IN  DEGREES*  BOTTOM  LA  YER  =  1 
0.0  90.0  30.0-30.0-30.0  30.0  90.0  0.0 

NUMBER  OF  LAYERS  FOR  CORRESPONDING  PLY  ORIENTATION 

1.0  1.0  1.0  1.0  1.0  1.0  1.0  1.0 

MATERIAL  PROPERTY  IDENTIFICATION  NUMBER  FOR  CORR.  PLY 
1  1  5  5  5  5  1  1 

TEMPERATURE  DT=  0.0000  MOISTURE=  0.0000 
SI  UNITS 


FAILURE  SURFACE  FOR  THIS  LAMINATE  HAS  BEEN  PLOTTED 
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Figure  A-4:  Failure  Envelopes  for  <0/90/+30)s  -  Hybrid  Laminates 

0°  and  90°  Plies  T300/5208  Material  and  +30°  Layers 
Kevlar  49  Material. 


VIII.  Failure  surfaces  for  different  laminates  in  one  figure: 

In  this  illustration  the  failure  envelopes  of  0°  laminate  on  the 
basis  of  six  failure  theories  are  given  for  T300/5208  material, 
using  the  option  of  SUPERPOSE  and  MULTICURV. 


Input  data 


LAMINATE  INPLANE 

PURE  STRNGTHPLTPLT START 

0  DEGREE  LAMINATE  T300/5208  MATERIAL  TSAI  WU  FSXY=-.5 
T300/5208 

SLAYER  NNL  =  1  »TH=i*0.  *PLNM  =  l*l7,Df=0~.  »C  =  0.  *NLDCN  =  l  $ 

TSAI  WU  STRAIN  SUPERPOSE  -.5  MULTICURV 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOW 

•  0  DEGREE  LAMINATE  T300/5208  MATERIAL  CHAMIS 
T30O/52O3 

SLAYER  NNL=l*TH=l*0.*PLNK=l*l.tDT=0#*C=0#»NLDCN=l  $ 

CHAMIS  STRAIN  SUPERPOSE  -.7  MULTICURV 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOW 

0  DEGREE  LAMINATE  T300/5208  MATERIAL  HOFFMAN 
T300/52G8 

SLAYER  NNL=1*TH=1*0.*PLNM=1*1.*DT=0«*C=0.*NLDCN=1  $ 

HOFFMAN  STRAIN  SUPERPOSE  MULTICURV 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOW 

0  DEGREE  LAMINATE  T300/5208  MATERIAL  HILL  CRITERIA 
T 30 0/52 08 

SLAYER  NNL=1,TH=1*0.*PLNM=1*1.,DT=0.*C=0.»NLDCN=1  $ 

HILL  STRAIN  SUPERPOSE  MULTICURV 

LAMINATE  INPLANE 

PURE  '  STRNGTHPLTPLTFOLLOW 

0  DEGREE  LAMINATE  T300/5208  MATERIAL  MAXSTRAIN  CRITERIA 
T300/5208 

SLAYER  NNL=1,TH=1*0.,PLNM=1*1.,DT=0.,C=0.,NLDCN=1  $ 

MAXSTRAIN  STRAIN  SUPERPOSE  MULTICURV 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTEND 

0  DEGREE  LAMINATE  T300/5208  MATERIAL  MAXSTRESS  CRITERIA 
T300/5208 

SLAYER  NNL=l*TH=l*0.*PLNM=l*l.tDT=0.fC=0.tNLDCN=l  $ 
MAXSTRESS  STRAIN  SUPERPOSE  MULTICURV 

THE END 

Output:  Figure  A- 5 
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IX.  Failure  envelopes  for  (0/90/+30)  laminate  in  different 
spaces  on  the  basis  of  Tsai  Wu  failure  criteria,  F*  =  -.5. 

1.  qr  -  stress  space 

2.  qr  -  strain  space 

3.  principal  stress  space 


Input  Data: 


LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTSTART 

(0/90/30/-30 I  SYMMETRIC  LAMINATE  TSAI  WU-STRESS  QR. 
T300/5208 

SLAYER  NNL=8*TH=0. *9  0. *30. *-30. *-30. *30. *90. *0.* 
PLNM=8*1 .*DT=0.*C=0. *NLDCN=1  $ 

_  TSAI  WU  STRESS  QR  -0.5 

LAMINATE  INPLANE 

PURE  STRNGT  HP  LTP*LT  FOLLOW 

CO/90/30/-30)  SYMMETRIC  LAMINATE  TSAI  WU-STRAIN  QR. 
T300/5208 

SLAYER  NNL=8*TH=0.*90.*30.*-30.*-30.*30.*9Q.*0.* 
PLNM=8*1.,DT=0.,C=0.*NLDCN=1  $ 

_  TSAI  WU  STRAIN  QR  -0.5 

~  LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTEND 

CO/90/30/-30)  SYMM.  LAMINATE  TSAI  WU-  SiTRESS  SPACE 
T300/5208 

SLAYER  NNL=8,TH=0.,90. * 30 • ,-30 . * -30 • *  30. *90 . *0 . » 
PLNM=8*1.*DT=0.*C=0.,NLDCN=1  $ 

TSAI  WU  STRESS  -0.5 

THEEND 


Output:  Figures  A-6  through  A-8. 
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Figure  A -6:  Failure  Envelopes  for  (0/90/+30)s  Laminate  in 
Plane  of  Stress  Space. 


X 


.  Laminates  consisting  of  more  than  three  orientations: 

1.  Failure  envelopes  of  0°,  30°  and  60°  layers  in  a 
(0/+30/+45/+60/+75/90)  laminate. 

■i“”  1  o  - 

2.  failure  Envelopes  of  45° ,  75°  and  90°  layers  in  a 
(0/+30/+45/+60/+75/90)  laminate. 

3.  Failure  envelopes  of  0°,  30°,  45°,  60°,  75°  and  90° 
layer  of  the  aforementioned  laminate. 


Input  data: 


LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTSTART 

(0/30/60/45/75/90/-75/-45/-60/-30)  SYM  LAM*  FAILURE  SURF. 
T300/5208 

SLAYER  NNL=20»TH=0.*3G.*60.*45.*75.*9Q.*-75. *-45  .  * -60 . *-3 0 . * 

—  30#  *-60*  *  —  45*  *  —  75**90 • *75*  *45*  *60*  *3  0 • *0*  * 

PLNM=20*1.*DT=0.*C=0.*NLDCN=1  $ 

TSAI  UU  STRAIN  -*5 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOW 

<0/30/60/45/75 /90/-75/-45/-60/-30)  SYM  LAM.  FAILURE  SURF. 
T300/5208 

SLAYER  NNL=20,TH=45. *75. *90. *0.*30. *60.* 

-60.*-30.*-75.*-45.* -45 .*-75.*“30.*“60.*60.*30.*0.*90.*75.*45»* 
PLNM=20*1.*DT=0.*C=0.*NLDCN=1  $ 

TSAI  UU  STRAIN  -.5 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTEND 

C0/30/60/45/75/90/-75/-45/-60/-30)  SYM  LAM.  FAILURE  SURF. 

T 30 0/52 08 

SLAYER  NNL=20*TH=0.,30.*60.*45. *75. *90. *-75. *-45.»-60.*-30.* 

-  3  0  •  *  -  6  0  .  *  -  4  5  •  *  -  75  •  *  9  0  .  *  7  5  •  *  4  5  .  *  6  0  •  *  3  0  •  *  0  •  * 

PLNM=2  0  * 1 • *  DT-0 • *C  =  0 •*  NLDCN  =  1  $ 

TSAI  UU  STRAIN  SUPERPOSE  -.5 

THE END 


Output:  Figures  A-9  through  A-11  • 
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Figure  A-9:  Failure  Envelopes  of  0°,  30°  and  60°  Plies 
in  (0/90/+30/+45/+60/+75)  -  Laminate  in 

Principal  Strain  Space. 


no 


Hill 


VAWi 


!»i 


i 


Figure  A-ll: 


Failure  Envelopes  of  all  the  Fly  Orientations 
for  Laminate  Considered  in  Figs.  A- 9,  and  A-10 
Superposition  of  Figs.  A-9  and  A-10. 


XI.  Failure  envelopes  for  (0/90/+45I  -  laminate  made  of  T300/5208 

—  s 

graphite  epoxy  material  on  the  basis  of  six  failure  theories  in 
principal  stress  and  principal  strain  spaces. 


Input  data: 

LAMINATE  INPLANE 

PURL  STRNGTHPLTPLTSTART 

( 0/9C/45/-45 )  LAMINATE 

T 30 6/52 08 

SLAYER  NNL  =  8*TH=0.*9C.,45.*2*-45,.,4':;.*90.,0.,PLNM  =  8*1.,DT  =  0.*C  =  0., 
NLDCM=1  $ 

TSAI  MU  STRESS  -.5 

LAMINATE  INFLAME 

PURE  STRNGTHPLTPLTFOLLOU 

(C/90/45/-45)  LAMINATE 
T  3  0  0  /  5  2  0  8 

SLAYER  NNL  =  8  ,  TH  =  0.  ,90.  ,45.  ,  2*-45‘.  ,45.  ,90.  ,0.  ,PLNM  =  8*1  .  ,DT  =  0.  ,C  =  0.  » 
NLDCN=1  S 

TSAI  MU  STRAIN  -.5 

LAMINATE  INPLANE 

PURL  STRNGTHPLTPLTFOLLOU 

■(0/90/45/-45J  LAMINATE 
T3 00/52 08 

SLAYER  NNL=8,TH=0.»90.»45. , 2 * -45  . ,45 . ,90  . *  0 . * PLNM  =  8 * 1 .  ♦  DT= 0  .  ,  C  =  0  .  , 
NLDCN'-l  S 

CHAMIS  STRESS  -.7 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOU 

< 0/90/45/-45)  LAMINATE 
T300/520P 

SLAYER  NNL-8*TH=0.*9Q. *45 . , 2  * -4 5 . , 45 . *  90 . *  0 . , PLNM  =  8  *  1  .  , DT= 0  .  * G  =  0 . * 
NLDCN=1  S 

CHAMIS  STRAIN  -.7 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOU 

C 0/90/45/-45)  LAMINATE 
T 30 0/52 08 

SLAYER  NNL  =  8  *TH=0  . ,90. ,45 . ,2*-45. ,45. *90 . ,0 . ,PLNM  =  8*1 . *DT=0. *C  =  0 . * 
N  L  0  C  N  =  1  % 

HOFFMAN  STRESS  0. 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOU 

< 0/90/45/-45)  LAMINATE 

T3G0/5208 

SLAYER  NNL  =  8,TH=0.,90.,45.,2*-45.,4  5.,90.-,0.,PLNM  =  8*l.,pT  =  0.,C  =  0», 
NLUCN^l  S 

HOFFMAN  STRAIN  0. 

laminate  inplane 

PURE  STRNGTHPLTPLTFOLLOU 

( 0/90/45/-45)  LAMINATE 
T30C/520 8 
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SLAYER  NNL=P*TH=0.*90. *45. *2*-45.*4f .,90.*0. *PLNM=8*1.*DT=0.*C=0.» 
NL0CN=1  $ 

HILL  STRESS 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOW 

( 0/90/ A5/-45)  LAMINATE 
T3G0/52D« 

SLAYER  NNL=fi*TH=0. *9C.*45.*2*-45. ,45. *90 .*0.  *  PLNM  =  8 *  1  .  , DT  =  0 . * C  =  0 . * 
N  L  D  C  N  =  1  $ 

HILL  STRAIN 

LAMINATE  INPLANE 

PURE.  STRNGTHPLTPLTFOLLOW 

( 0/90/45/-45)  LAMINATE 
T  3  0  0  /  5  2  0  8 

SLAYER  NNL=R»TH=0.*9O.*45.*2*-45.*45.,9O. , 0 . *PLNM=8 * 1 . , DT= 0 . *C=0. * 
NLDCN-1  S 
MAXSTRAIN  STRESS 
LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOW 

<  0/9Cl/45/-45>  LAMINATE 
T300/5208 

SLAYER  NNL=fi*TH=0.*90.*45.*2*-45.*45.*90.*0.*PLNM=8*l.,DT=0.  *C=0.* 
NLDCN=1  S 
MAXSTRAIN  STRAIN 
LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOW 

(D/90/45/-45)  LAMINATE 
T300/52QR 

SLAYER  NNL=R*TH=0.,90.,45.*2+-45..45.,9O.*O.*PLNM-8*1.*DT=O.*C=O.* 
N  L  D  0  N  =  1  S 
MAXSTRESS  STRESS 
LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTEND 

( 0/90/45/-45)  LAMINATE 
T  3  0  0  /  e  2  0  8 

SLAYER  WNL  =  R  » TH=0 . *9G. *45. *2*-45. *45. *90.*  0. *PLNM  =  8*1 . *DT  =  0.  *C  =  0.  ♦ 
NLDCN= 1  S 
MAXSTRESS  STRAIN 
THEEND 

Output:  Figures  A  - 1 2  through  A-23. 
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Figure  A -12:  Failure  Envelopes  for  (0/90/+45)  -  Laminate  of 

■“  s 

T300/5208  Material  in  Stress  Space  on  the  Basis 
of  Tsai  Wu  Criterion. 
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Figure  A-14:  Failure  Envelopes  for  (Q/90/+45)  -  Laminate  of 

“  s 

T300/5208  Material  in  Stress  Space  on  the  Basis 
of  Chamis  Criterion. 
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Figure  A-15:  Failure  Envelopes  for  (0/90/+45)  -  Laminate  of 

"  s 

T300/5208  Material  in  Strain  Space  on  the  Basis 
of  Chamis  Criterion. 
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Figure  A-17:  Failure  Envelopes  for  (0/90/+45)  -  Laminate 

of  T300/5208  Material  in  Strain  Space  on  the 
Basis  of  Hoffman  Criterion. 
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Failure  Envelopes  for  (0/90/+45)  -  Laminate  of 

“  s 

T300/5208  Material  in  Strain  Space  on  the  Basis 
of  Hill  Criterion. 


Failure  Envelopes  for  (0/90/+45)  -  Laminate  of 

s 

T300/5208  Material  in  Strain  Space  on  the  Basis 
of  Max.  Stress  Criterion. 
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Figure  A-22:  Failure  Envelopes  for  (0/90/+45)  -  Laminate  of 

.  g 

T300/5208  Material  in  Stress  Space  on  the  Basis 
of  Max.  Strain  Criterion. 
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XII.  Change  of  Scale: 

The  scale  of  the  coordinate  axis  of  Figure,  (A-18)  has  been  reduced 
to  1/2  by  the  use  of  FCTRS. 


Input  data: 


LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTONE 

C0/90/45/-45)  LAMINATE 
T  300/52  08 

SLAYER  NNL=8*TH=0,#90.*45. t2*-45. t45. *90.*Q.*PLNM 
NLDCN=1  $ 

HILL  STRESS 

THEEND 


Output:  Figure  A-24. 


=8*1.*DT=0.*C=0., 

.5 
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Figure  A-24;  Failure  Envelopes  for  (0/90/+45)s  Laminate  of 

T300/5208  Material  in  Stress  Space  on  the 
Basis  of  Hill  Criterion.  Same  Figure  as  A-18 
with  Scale  Reduced  to  Half. 


Failure  surfaces  for  isotropic  materials: 

Isotropic  material  with  v  =  .5  and  X'  =  X,  Y'  =  Y,  in 
principal  strain  space  for  =  -0.5 

Isotropic  material  with  v  =  0,  (0.5)  and  X'  =  X,  Y'  =  Y,  in 

principal  stress  space  for  F*^  =  0,  (-.5). 

Isotropic  material  with  v  =  0  and  X'  =  X,  Y'  =  Y,  in 
principal  stress  space,  F*^  =  -.5. 

Isotropic  material  with  v  =  0,  (.5)  and  X'  =  2X,  Y'  =  2Y, 
in  principal  stress  space  for  F*^  =  -.5,  (0). 


Input  data: 

NEWMTRLS 


YLWMIRLS  SI 

£L  A MOAT  A  NNM=4tEX=4*69.*EY=4*69.»VX=0.,.5»0.,.5» 
ES=34.5t23.»34.5,23.»ALFX=4*0.*ALFY=4*0.tBTAX=4*0.< 
BTAY=4*0.tX=4*400.*XD=2*4C0. *2*800. »Y=4*400., 
YD=2*400. *2*800., 5=4*230. ,SK=. 000125  $ 

AM  7  m  a  t  r  INPLANE 


YD=2*400., 2*800. ,S=4*230.,SK=. 000125  $ 

LAMINATE  INPLANE 

PURE  S TRNGTHPLT PL T START 

V=.5,X=X0 
B4/5505 

SLAYER  NNL=lfTH=0.*PLNM=l.*DT=0.tC=0.*NLDCN=l  $ 
TSAI  WU  STRAIN  SUPERPOSE  -.5 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOW 

V=0,X=XD,FSXY=-.5 
T300/5208 

SLAYER  NNL=1,TH=0.*PLNM=1.,DT=0.,C=0.,NLDCN=1  $ 

TSAI  WU  STRESS  SUPERPOSE  -.5 

LAMINATE  INPLANE 

PURE  STRNGTHPLTPLTFOLLOW 

V=.5»X=XD,FSXY=*0 
B4/5505 

SLAYER  NNL=1*TH=0.*PLNM=1.*DT=0.,C=0.,NLDCN=1  S 
TSAI  WU  STRESS  SUPERPOSE  0.0 


MULTICURV 


AM  INATE 
JR  E 


STRESS 

INPLANE 


0.  25 


STRNGTHPLTPLTFOLLOW 


V=0,X=XD,FSXY=-.5 
5  0  C/52  Cb 

iL  AYER  NNL  =  1  *TH  =  0  .  ,PLNM=1  .  ,DT  =  0.  ,C=0.  ,NLDCN= 
J  A I  WU  STRESS  SUPERPOSE  -.5 

iM  INATE'  INPLANE 


0.25 


L. 

V=G,XD=2X,FSXY=-.5 
1 3  5  0 1 

-AYER  NNL=1 *TH=0. *PLNM=1 . , 
A  I  WU  STRESS 
1 1  NATE  INPLANE 
*  E 

/=.5,XD=2X*FSXY=  Q . 0 
3TCHPLY 

.AYER  NNL=1*TH=0. *PLNM=1 . , 
U  WU  STRESS 


STRNGTHPLTPLTFOLLOW 

DT=0.»C=0. »  N  L DC  N= 1  $ 

SUPERPOSE  -.5 

STRNGTHPLTPLTEND 

DT=0.*C=0. »NLDCN=1  S 
SUPERPOSE  0.0 


MULTICURV  0.5 


Output:  Figures  A- 2 6  through  A-29. 
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Figure  A-27:  Failure  Surfaces  for  an  Isotropic 
Material  in  Principal  Stress  Space 
=  0  (-.5)  and  F*  =  0,  (-.5) . 


V 
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Figure  A -2 9: 


Failure  Surfaces  for  an  Isotropic 
Material  with  v  =  0,  (.5),  X'  =  2x, 

Y1  =  2y  in  Principal  Stress  Space  for 
Fjy  =  0,  (-.5). 
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APPENDIX  B 
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CURRENT  AFWAL/MLBM  LAMINATE  PROGRAMS 

MAY  1983 


SYMMETRIC 

ONLY 

BALANCED 

ONLY 

IN-PLANE 

STIFFNESS 

FLEXURAL 

STIFFNESS 

LAMINATE 

STRENGTH 

THERMAL 

EFFECT 

HYGROSCOPIC 

EFFECT 

HARDCOPY 

MEDIUM 

STORAGE 

MEDIUM 

OPTIMIZATION 

(SIZING) 

SANDWICH 

CORE 

HYBRID 

MATERIALS 

FORTRAN 

(CDC) 

X 

X 

1,2 

X 

X 

4 

10 

X 

X 

Ti  CC-40 

X 

16 

16 

16 

16 

16 

14,15 

16 

T 1  -59 

W/  COMBO  CARDS 

X 

X 

X 

1 

5 

11 

X 

14 

APPLE  n* 

X 

X 

1.3 

X 

X 

6 

12 

X 

TIMEX  1000 
(+16K) 

X 

7 

13 

14 

TRS-80 

PC-1 

X 

X 

X 

1 

5 

13 

X 

SHARP  PC-1500 
(+8K) 

X 

X 

X 

1 

5,8 

13 

X 

EPSON  HX-20 

X 

X 

X 

1 

8,9 

13 

14 

X 

PANASONIC 

HHC  (4K) 

X 

X 

X 

1 

5 

X 

1. 

QUADRATIC  (TSAI -WU)  FAILURE 

7. 

TIMEX/SINCLAIR  COMPATIBLE 

CRITERION 

PRINTERS 

2. 

REPORT  IN  PROGRESS  INCLUDES 

8. 

ALSO  PLOTS  FAILURE  SURFACES 

ADDITIONAL  FAILURE  CRITERIA: 

9. 

BUILT-IN  PRINTER 

a.  CHAMIS 

b.  HOFFMAN 

10. 

PUNCHED  CARDS /MAG.  TAPE 

c.  HILL 

11. 

REQUIRES  CUSTOM  AFWAL/MLBM 

d.  MAX  STRESS 

ROM  MODULE  AND  MAG.  CARDS 

e.  MAX  STRAIN 

12. 

5  1/4"  DISK 

3. 

MAX  STRAIN  FAILURE  CRITERION 

13. 

AUDIO  CASSETTE  TAPE 

4. 

LINE  PRINTER  AND  CALCOMP 

FAILURE  SURFACE  PLOTS 

14. 

IN-PLANE  LOADS  ONLY 

5. 

PRINTING  CRADLE 

15. 

BALANCED  LAMINATES  ONLY 

6. 

80  COLUMN  PRINTER 

16. 

PROGRAMMING  IN, PROGRESS 

ALSO  AVAILABLE  ARE  DISKS  CONTAINING  PROGRAMS  TO  ANALYZE  BEAMS,  SHAFTS, 
PRESSURE  VESSELS,  MOISTURE  ABSORPTION  AND  DESORPTION,  CURING  OF  EPOXY-MATRIX 

COMPOSITES,  AND  LIFE  PREDICTION  (COMING  SOON). 
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